Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,679,639 A→C 7.1% Q28P (CAA→CCA)  ydgI → putative arginine/ornithine antiporter transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,679,6390AC 7.1% 267.5 / 8.1 93Q28P (CAA→CCA) ydgIputative arginine/ornithine antiporter transporter
Reads supporting (aligned to +/- strand):  ref base A (43/43);  new base C (5/2);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

ACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGC  >  NC_000913/1679491‑1679780
                                                                                                                                                    |                                                                                                                                             
aCTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCaa                                                                                                                                              <  1:1425665/150‑1 (MQ=255)
 cTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGCCTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCaaa                                                                                                                                             >  2:1047318/1‑150 (MQ=255)
  tCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCaaaa                                                                                                                                            >  2:1586854/1‑150 (MQ=255)
      tttttATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGGTGATATTAAGCTCAATGATGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATg                                                                                                                                        >  2:36578/1‑150 (MQ=255)
        tttATGAGTTTTGCCCTAATATTTATTATTATTATCACGTTTTAATCACTGGATATCGAAGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAAAGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATggc                                                                                                                                      <  2:1018804/150‑1 (MQ=255)
          tATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATggcgg                                                                                                                                    >  1:536254/1‑150 (MQ=255)
           ttGCTTTTTGCCCTTATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCTATGGAAAAGAAACTGGGACTGAGCGGACTTACGGCTCTGGTATTAAGCTCTATGCTGGGCGCGTTTTTTTTCAGTCTGCCGCAAAATATggcggc                                                                                                                                   <  1:1687398/149‑1 (MQ=255)
               gTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCCATATTTTGCGGCAgac                                                                                                                               <  1:1796677/150‑3 (MQ=255)
                ttttGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTg                                                                                                                              <  1:13482/150‑1 (MQ=255)
                 tttGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCCATATTTTGCGGCAgactg                                                                                                                             <  2:1591606/150‑5 (MQ=255)
                  ttGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGcc                                                                                                                            <  2:1280606/150‑1 (MQ=255)
                    gCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAg                                                                                                                          >  1:695198/1‑150 (MQ=255)
                    gCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAg                                                                                                                          <  1:1381336/150‑1 (MQ=255)
                     cccTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGc                                                                                                                         <  1:703812/150‑1 (MQ=255)
                           tatttattCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGgca                                                                                                                   >  1:850292/1‑150 (MQ=255)
                           tatttattCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGAGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAATATATGGCGGCAGTTGCCAGCCCGgcc                                                                                                                   >  1:1256955/1‑149 (MQ=255)
                            atttattCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGgcag                                                                                                                  >  2:446595/1‑150 (MQ=255)
                            atttattCTTCTTATCACGTTTAAATCACTGGATATCGATCGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGgcag                                                                                                                  <  1:378933/150‑1 (MQ=255)
                             tttattCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGgcagc                                                                                                                 <  1:695327/150‑1 (MQ=255)
                                  tcttctTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTGAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGc                                                                                                            >  2:1110714/1‑150 (MQ=255)
                                  tcttctTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGc                                                                                                            >  1:856129/1‑150 (MQ=255)
                                   cttcttATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCt                                                                                                           >  2:342923/1‑150 (MQ=255)
                                    ttcttATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTc                                                                                                          <  2:1104093/150‑1 (MQ=255)
                                         aTCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCgg                                                                                                     >  2:1652207/1‑150 (MQ=255)
                                         aTCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCgg                                                                                                     >  2:894774/1‑150 (MQ=255)
                                         aTCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCgg                                                                                                     >  1:219728/1‑150 (MQ=255)
                                           cACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCt                                                                                                   >  2:1104820/1‑150 (MQ=255)
                                            aCGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTg                                                                                                  <  2:1351812/150‑1 (MQ=255)
                                             cGTTTTAATCACTGGATATCGATGGATAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTgg                                                                                                 <  2:1512037/150‑1 (MQ=255)
                                               ttttAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTgggg                                                                                               <  1:217548/150‑1 (MQ=255)
                                                 ttAATCACTGGATATCGATGGATAAGAACCTGGGACTGCGCGCACTCACCGCGCTGGGATTAAGCTCAATGGTGGGCGCGGGTGTTTTCACTCTGCCGCAAAATATGGCGGCAGTTTCCAGCCCGGCCGCAGGGCTAAGCGGCGGGGGTa                                                                                             >  2:1151713/1‑150 (MQ=255)
                                                   aaTCACTGGATATCGATGGAAAAGAAACTGGGACTGCGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTAtt                                                                                           >  2:1787224/1‑150 (MQ=255)
                                                    aTCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTa                                                                                          >  1:1544782/1‑150 (MQ=255)
                                                        cTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTgg                                                                                      <  1:1635096/150‑1 (MQ=255)
                                                          ggATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGgcg                                                                                    <  2:446676/150‑1 (MQ=255)
                                                          ggATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGgcg                                                                                    >  1:1546562/1‑150 (MQ=255)
                                                              atCGATGGACAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTgg                                                                                <  1:104742/150‑1 (MQ=255)
                                                              atCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTgg                                                                                <  1:530222/150‑1 (MQ=255)
                                                                 gATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCAGCGGCTGGGGTATTACTGGCGCTGGCAt                                                                             >  2:251764/1‑150 (MQ=255)
                                                                 gATGGAAAAGAAACTGGGACTGAGCGAACACACCGCGCTGCTATAAAGCTCAAAGCTGGGCGCGGGTGTTTTCAGTCTCCCGCAAAATAAGGCGGCACTTGCCAGCCCGGCAGCACTGCTCATCGCCTGGGGTATTACTGGCGCTGTCaa                                                                             >  2:1804404/1‑149 (MQ=255)
                                                                  aTGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCAtt                                                                            >  1:1703235/1‑150 (MQ=255)
                                                                     gAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATtata                                                                         >  1:541997/1‑150 (MQ=255)
                                                                     gAAAAGAAACTGGGACTGAGCGCACTCACCGCGCAGGTATTAAGCTCTATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCATCCCGGCAGCCCTGCTCATCGTCTGGGGTATTACTGGCGCTGGCAGTTTa                                                                         >  2:1640121/1‑150 (MQ=255)
                                                                      aacaGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTAt                                                                        <  1:383259/147‑1 (MQ=255)
                                                                       aaaGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGACAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTAtt                                                                       <  2:568622/150‑1 (MQ=255)
                                                                          gAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCt                                                                    >  1:780268/1‑150 (MQ=255)
                                                                            aaCTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTgg                                                                  >  1:1436646/1‑150 (MQ=255)
                                                                              cTGGGACTGAGCGCACTCACCTCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGcc                                                                <  2:1446369/150‑1 (MQ=255)
                                                                              cTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGcc                                                                <  2:490548/150‑1 (MQ=255)
                                                                              cTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGcc                                                                <  2:1533528/150‑1 (MQ=255)
                                                                               tGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCTCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGACt                                                               >  1:1510713/1‑150 (MQ=255)
                                                                                  gACTGCGCGCACTCACCGCGCTGGGATTAAGCTCAATGCTGGCCGCGGGTGGTTTCAGGCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGAATTACTGGCGCTGGCATTTTATTGCTGGCCTTTg                                                            <  2:277100/150‑1 (MQ=255)
                                                                                  gACTGAGCGCACTCACCGCGCTTGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTg                                                            <  1:573331/150‑1 (MQ=255)
                                                                                  gACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTg                                                            >  1:1587389/1‑150 (MQ=255)
                                                                                    cTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGcc                                                          <  1:1367652/150‑1 (MQ=255)
                                                                                    cTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGcc                                                          <  1:1257721/150‑1 (MQ=255)
                                                                                      gAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCAt                                                        >  2:307076/1‑150 (MQ=255)
                                                                                      gAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGGGTTTTCAGTCTGACGCAAAATATGTCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCAGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCAt                                                        >  2:193336/1‑150 (MQ=255)
                                                                                         cgcACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGGGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCt                                                     <  2:1423654/150‑1 (MQ=255)
                                                                                            aCTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTACCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGAt                                                  <  2:563695/150‑1 (MQ=255)
                                                                                             cTCACCTCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGTGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATc                                                 <  1:1529387/150‑1 (MQ=255)
                                                                                             cTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCCCTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATc                                                 >  1:1145031/1‑150 (MQ=255)
                                                                                               cACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCt                                               >  2:539201/1‑150 (MQ=255)
                                                                                                 ccGCGCTGGTATAAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCa                                             <  2:1659946/150‑1 (MQ=255)
                                                                                                      cTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACgcgc                                        >  1:518384/1‑150 (MQ=255)
                                                                                                        ggTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATACGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCAt                                      <  2:246527/150‑1 (MQ=255)
                                                                                                          tATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCCCGCGCATTc                                    >  1:848155/1‑150 (MQ=255)
                                                                                                          tATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTc                                    <  2:1205503/150‑1 (MQ=255)
                                                                                                          tATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTc                                    >  1:413702/1‑150 (MQ=255)
                                                                                                           aTTAAGCTCAAGGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCg                                   <  2:1677108/150‑1 (MQ=255)
                                                                                                               aGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTccc                               >  2:1488280/1‑150 (MQ=255)
                                                                                                               aGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTccc                               >  1:1661835/1‑150 (MQ=255)
                                                                                                                gCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCg                              >  2:680620/1‑150 (MQ=255)
                                                                                                                gCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCg                              >  1:1577635/1‑150 (MQ=255)
                                                                                                                gCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGAATTTTATTGCTGGCCATTGCCATGCTGATCCTCACGCGCATTCGTcccc                              >  2:1513914/1‑149 (MQ=255)
                                                                                                                  tCACTGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGaa                            <  2:1122762/150‑1 (MQ=255)
                                                                                                                  tCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGGTCCTCACGCGCATTCGTCCCGaa                            >  1:1318966/1‑150 (MQ=255)
                                                                                                                   cAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAAc                           <  1:770219/150‑1 (MQ=255)
                                                                                                                     aTGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACtt                         >  1:654583/1‑150 (MQ=255)
                                                                                                                       gCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGa                       <  2:1861986/150‑1 (MQ=255)
                                                                                                                        cTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAc                      <  1:1277293/150‑1 (MQ=255)
                                                                                                                        cTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAc                      <  1:219435/150‑1 (MQ=255)
                                                                                                                        cTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAc                      <  1:1844211/150‑1 (MQ=255)
                                                                                                                             cgcgGGTGTTTTCAGTCAGCCGCAAAATATGGCGGCTGTTGCCAGCCCGGCAGCACTGCTGATCGGCTGTGGTATTACTGGCGCTGGCATTTTTTTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAcggcgg                 <  1:1565276/150‑1 (MQ=255)
                                                                                                                              gcgGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGt                <  2:1770252/150‑1 (MQ=255)
                                                                                                                                   tgtTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCtt           >  1:510133/1‑150 (MQ=255)
                                                                                                                                      tttCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTAc        <  1:1652207/150‑1 (MQ=255)
                                                                                                                                       ttCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTAcc       <  2:850292/150‑1 (MQ=255)
                                                                                                                                       ttCAGGCTGCCGAAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTAcc       <  1:1682199/150‑1 (MQ=255)
                                                                                                                                        tCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCCCTGCTCATCGGCGGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTCCCATGCTGATCCTCACGCGCCTTCGTCCCGAACTTGACGGTGGTATCTTTACCt      >  2:949565/1‑150 (MQ=255)
                                                                                                                                        ggcaacTGCCGCCATATTTTGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCt      >  1:999022/6‑150 (MQ=255)
                                                                                                                                          caacTGCCGCCATATTTTGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTAt    >  2:1002691/4‑150 (MQ=255)
                                                                                                                                          aGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGTTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTAt    >  1:874642/1‑150 (MQ=255)
                                                                                                                                           aacTGCCGCCATATTTTGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATg   >  1:497816/3‑150 (MQ=255)
                                                                                                                                            acTGCCGCCATATTTTGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGc  >  1:715053/2‑150 (MQ=255)
                                                                                                                                            acTGCCGCCATATTTTGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGAACCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGc  >  2:1366852/2‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
ACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGC  >  NC_000913/1679491‑1679780

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: