Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 73,319 T→C P401P (CCA→CCG)  thiP ← thiamine ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091373,3190TC97.4% 119.0 / ‑6.7 38P401P (CCA→CCG) thiPthiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (20/17);  minor base G (0/1);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGA‑CGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCAC  >  NC_000913/73073‑73529
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                   
gAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGA‑CGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCa                                                                                                                                                                                                                 <  1:111855/251‑1 (MQ=255)
         caccacGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGC‑CGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGccaatgg                                                                                                                                                                                                        >  2:24006/1‑246 (MQ=255)
          ccccccccAAAATCACCCACCGACAGCACGCCGGCCAAAGCCAGCCCCTGCGCCAGTGGAACGTT‑CATGGCCCTCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCTGCAGGccaatggt                                                                                                                                                                                                       <  2:105180/245‑7 (MQ=255)
          accacGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGTA‑CGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGccaatggt                                                                                                                                                                                                       <  1:33855/251‑7 (MQ=255)
              cGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGA‑CGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGccaatggggttg                                                                                                                                                                                                   >  2:52532/1‑241 (MQ=255)
                       cACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGA‑CGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCATCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGt                                                                                                                                                                                         <  2:22919/251‑1 (MQ=255)
                              cGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGA‑CGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTaaaa                                                                                                                                                                                     >  1:106517/1‑248 (MQ=255)
                                  aGCACGCAGGCACAAGCCAGAGCCTGCGCCCGTGGA‑CGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑AGTTGTTGAGCAGTAAAAAGAAGcc                                                                                                                                                                              >  2:166586/1‑251 (MQ=255)
                                    cACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGA‑CGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCgg                                                                                                                                                                            <  1:21783/251‑1 (MQ=255)
                                       gcaggCAAAAGCCAGCGCCTGCGCCAGTGGA‑CGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGGCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGTGCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAg                                                                                                                                                                         <  2:254273/247‑1 (MQ=255)
                                         aggCAAAAGCCAGCGCCTGCGCCAGTGGA‑CGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGcc                                                                                                                                                                       >  1:186676/3‑251 (MQ=255)
                                           gCAAAAGCCAGCGCCTGCGCCAGTGGA‑CGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAg                                                                                                                                                                     <  2:53527/251‑1 (MQ=255)
                                             aaaagCCAGCGCCTGCGCCAGTGGA‑CGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAAATG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGc                                                                                                                                                                    >  1:211683/5‑250 (MQ=255)
                                             aaaagCCAGCGCCTGCGCCAGTGGA‑CGTTTC‑GGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGcac                                                                                                                                                                  <  1:119802/247‑1 (MQ=255)
                                              aaagCCAGCGCCTGCGCCAGTGGA‑CGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGcac                                                                                                                                                                  >  2:38646/4‑251 (MQ=255)
                                              aaagCCACCGCCTGCGCCAGAGGA‑CGCTTCAGGGCGCGCCGCTCGACCACTTTTAAGCCAGACCAGCCTTCAATGCCCAGCGACTGCCACAACATACTGTCGCGGTCGGTGATATCGCGCATCGGGTTTTCCCGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGCGCCGTAAAAATCAGCCGGTAACCAGcac                                                                                                                                                                  >  2:170855/4‑251 (MQ=255)
                                                  ccAGCGCCTGCGCCAGTGGA‑CGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATc                                                                                                                                                              >  1:94499/1‑251 (MQ=255)
                                                        ccTGCGCCAGTGGA‑AGTTTCAGGGGGAGCAGCTCGACCCCTTTTAAGCGAGACACGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTAATTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATAG‑GGTTGTTGAGCAGTAAAAAGACGCCGGTAACCCGCCCACTCCCCGGc                                                                                                                                                        >  2:259086/1‑251 (MQ=255)
                                                           gcgcCAGTGGA‑CGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCa                                                                                                                                                       <  1:158149/249‑1 (MQ=255)
                                                                   ggA‑CGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTACGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAt                                                                                                                                             >  1:44298/1‑251 (MQ=255)
                                                                                 cgcgcAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGGAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTc                                                                                                                                >  1:192086/1‑251 (MQ=255)
                                                                                          cgaccacTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCa                                                                                                                         <  1:20842/247‑1 (MQ=255)
                                                                                           gaccacTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGcg                                                                                                                      <  1:132009/250‑2 (MQ=255)
                                                                                                  tttAAGCGCGCCCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTACAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTATCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGAc                                                                                                                <  2:203310/250‑1 (MQ=255)
                                                                                                    tAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGAccc                                                                                                              <  1:226151/250‑1 (MQ=255)
                                                                                                    aaaGCGCGACCAGCCTTCAAGGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCg                                                                                                             <  2:44298/250‑1 (MQ=255)
                                                                                                                                 aCTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTcg                                                                                >  1:145998/1‑251 (MQ=255)
                                                                                                                                    gACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTcgcg                                                                              >  2:30629/1‑250 (MQ=255)
                                                                                                                                                        gtgttaatatGGCGCTTCGGGTTTTCCAGTACTTTCAGCGAAAAAGGGTTCGCCATTACCGCATTGGTGAAAA‑GCACAATGCGGTCAGAAGATGGCGGAAGGCCAATGT‑TGTTGTTGAGAAGTAAAAAGAAGCCGGTAGCCAGCACAACCCCCGGCAGGGCGAGAAGCAACATGCCGCTCACCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGcatcat                                                         <  1:257208/245‑1 (MQ=255)
                                                                                                                                                                     gcATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAAAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACg                                            >  1:48206/1‑251 (MQ=255)
                                                                                                                                                                      catcggtttgtcgattattttCAGCGCATAAGGGTTCGCCTTTACCTCATTG‑TGAAAAATCAAAATCCCGTCAGCAGATGGGGGCGGGCC‑ATGGGTGTTGTTGAGCAGTAAAAAAAGGCCGGTGGCCAGCAAAATCCCGGGCATGGCGAGGATAAACGGCCCGCTCACCCCCAGCGCCTGACCCGCCAGCATTTTCGGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc                                           <  2:167515/233‑1 (MQ=255)
                                                                                                                                                                         cGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGcac                                         >  1:209332/1‑250 (MQ=255)
                                                                                                                                                                                                 aTAAGGGATCGCCATTAACGCATTGGTGAAAG‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATAc                >  2:117618/1‑251 (MQ=255)
                                                                                                                                                                                                 aTAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATAc                >  1:113269/1‑251 (MQ=255)
                                                                                                                                                                                                    ggggATCCCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCACCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCCCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCa             <  1:88349/250‑1 (MQ=255)
                                                                                                                                                                                                          cgcattaaacgcttgggtaaaaa‑TCACAATCCCTTAAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCACCAAAACCCCCGGCAGGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAgg       <  1:179583/232‑1 (MQ=255)
                                                                                                                                                                                                          cGCCATTAACGCATTGGTGAAAA‑TCACCATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGCTCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAgg       >  2:30140/1‑251 (MQ=255)
                                                                                                                                                                                                               ttAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCcac  >  2:36085/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                   
GAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGA‑CGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAA‑TCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAG‑TGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCAC  >  NC_000913/73073‑73529

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑