Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,985,112 C→G 9.5% intergenic (‑43/+27) araG ← / ← araF L‑arabinose ABC transporter ATPase/L‑arabinose ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,985,1120CG 9.5% 220.8 / 11.4 72intergenic (‑43/+27)araG/araFL‑arabinose ABC transporter ATPase/L‑arabinose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (29/35);  new base G (4/3);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

TGAACCTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTAAAG  >  NC_000913/1984963‑1985261
                                                                                                                                                     |                                                                                                                                                      
tGAACCTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGcatc                                                                                                                                                        >  1:1380172/1‑150 (MQ=255)
tGAACCTGACCGGCATAGCAGTCAAAACTAATATCCGGCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGcatc                                                                                                                                                        <  2:841690/150‑1 (MQ=255)
tGAACCTGACCGGCATAGCAGTCAAAACTAAACTCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACAGATGCCGCGAAAGGATAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGcatc                                                                                                                                                        <  1:906692/150‑1 (MQ=255)
   aCCTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATg                                                                                                                                                     <  2:430684/150‑1 (MQ=255)
    ccTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACCCCCCCTCTGCATCATGc                                                                                                                                                    >  1:1487117/1‑150 (MQ=255)
    ccTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATgt                                                                                                                                                    >  1:776477/1‑149 (MQ=255)
    ccTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGc                                                                                                                                                    <  1:1606776/150‑1 (MQ=255)
    ccTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGc                                                                                                                                                    <  2:487859/150‑1 (MQ=255)
    ccTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGc                                                                                                                                                    <  2:1106184/150‑1 (MQ=255)
     cTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCa                                                                                                                                                   >  2:988196/1‑150 (MQ=255)
     cTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTGTTACCGCTGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCa                                                                                                                                                   <  2:1032264/150‑1 (MQ=255)
      tGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATACCGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCag                                                                                                                                                  >  1:218828/1‑150 (MQ=255)
       gACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCaga                                                                                                                                                 >  2:1353895/1‑150 (MQ=255)
        aCCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCagag                                                                                                                                                <  1:725891/150‑1 (MQ=255)
        aCCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCagag                                                                                                                                                <  2:1610132/150‑1 (MQ=255)
          cGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggg                                                                                                                                              <  1:988196/150‑1 (MQ=255)
          cGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGGTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGGCGTTCACACCCCCTCTGCATCATGCAGAggg                                                                                                                                              <  1:1464527/150‑1 (MQ=255)
            gCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTGGCATAACGTCTCCGTGACTCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggggg                                                                                                                                            <  2:1498814/150‑2 (MQ=255)
             cATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAgggggt                                                                                                                                           <  2:1531005/150‑3 (MQ=255)
                aGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAAttt                                                                                                                                        <  2:953293/150‑1 (MQ=255)
                      cAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCgg                                                                                                                                  >  2:569953/1‑150 (MQ=255)
                       aaaaCTAATATCCGTCCGCGCCTTACCGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGGAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGc                                                                                                                                 <  1:1480684/150‑1 (MQ=255)
                              atctCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAttac                                                                                                                          <  2:1189262/147‑1 (MQ=255)
                              atatCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAttac                                                                                                                          >  2:1349026/1‑150 (MQ=255)
                                  ccGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGCAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCACCTATGCATCATGCAGAGGGGAATTTTTCCGGGATAttacttac                                                                                                                      >  1:848777/1‑150 (MQ=255)
                                  ccGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAttacttac                                                                                                                      >  2:556402/1‑150 (MQ=255)
                                    gTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCg                                                                                                                    >  1:95838/1‑150 (MQ=255)
                                        gcgcCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTa                                                                                                                <  1:402736/150‑1 (MQ=255)
                                         cgcCTTAACGCCGGGAAACGTTTTACCGAAGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTATTCCTGCAAATTACTTACCGCCTaa                                                                                                               >  2:1413067/1‑150 (MQ=255)
                                               aaCGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑ttt                                                                                                        <  1:447940/150‑1 (MQ=255)
                                               aaCGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑ttt                                                                                                        >  2:1425403/1‑150 (MQ=255)
                                                aCGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑tttt                                                                                                       >  1:742050/1‑150 (MQ=255)
                                                    cGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTtct                                                                                                   >  1:394755/1‑150 (MQ=255)
                                                       gAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGCTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCa                                                                                                <  1:1349290/150‑1 (MQ=255)
                                                                tACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTcttctt                                                                                       <  1:1349026/150‑1 (MQ=255)
                                                                 aCCGCTGCCGCGAAATGATAGTTTCGGGGTAGACTTTTGGATAACGTCTCAGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGGAGAGGGGATTTTTTCCGGCAAATTAGTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCttt                                                                                      <  1:1291541/150‑1 (MQ=255)
                                                                  ccGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTcttctt                                                                                       <  2:1423506/148‑1 (MQ=255)
                                                                  ccGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTcttctt                                                                                       >  1:1423506/1‑148 (MQ=255)
                                                                   cGATGCCGGGTAATGAGAGATACGGGGTAGACTGTTGCATATATTCTCAGTGAATCACTGGTCGTTCACACCCCCTCTGCATAATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTaa                                                                                    <  1:1703325/150‑1 (MQ=255)
                                                                    gATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACACCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTaaa                                                                                   <  2:697558/150‑1 (MQ=255)
                                                                     aTGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAg                                                                                  <  1:1620813/150‑1 (MQ=255)
                                                                        ccGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTg                                                                               >  1:686826/1‑150 (MQ=255)
                                                                          gcgAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTc                                                                             <  2:1487117/150‑1 (MQ=255)
                                                                           cgAAATGAGAGATACGGGGGAGACTGTTGCATAACGGCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCa                                                                            <  2:1488750/150‑1 (MQ=255)
                                                                                tGAGAGATACGGGGTAGGCTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCGCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAGGTTGTCACGCGt                                                                       <  1:1606115/150‑1 (MQ=255)
                                                                                tGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGt                                                                       <  1:1666755/150‑1 (MQ=255)
                                                                                    agaTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATc                                                                   <  2:884887/150‑1 (MQ=255)
                                                                                     gaTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCa                                                                  >  2:577143/1‑150 (MQ=255)
                                                                                             gTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TGTATTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACg                                                          >  1:30697/1‑150 (MQ=255)
                                                                                               aGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTc                                                        <  2:780838/150‑1 (MQ=255)
                                                                                                 aCTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCgg                                                      <  1:983534/150‑1 (MQ=255)
                                                                                                 aCTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCgg                                                      >  1:1121925/1‑150 (MQ=255)
                                                                                                  cTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGt                                                     >  1:356216/1‑150 (MQ=255)
                                                                                                     ttGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAAc                                                  <  1:556402/150‑1 (MQ=255)
                                                                                                        cATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTc                                               >  2:1233266/1‑150 (MQ=255)
                                                                                                        cATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTc                                               >  1:1145000/1‑150 (MQ=255)
                                                                                                          tAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCg                                              >  2:492354/1‑150 (MQ=255)
                                                                                                           aaCGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGTAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGt                                            <  2:776477/150‑1 (MQ=255)
                                                                                                                  ccGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAAtttt                                     >  2:1696540/1‑150 (MQ=255)
                                                                                                                   cGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTg                                    >  1:522548/1‑150 (MQ=255)
                                                                                                                    gTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTgg                                   >  2:1011054/1‑150 (MQ=255)
                                                                                                                     tGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGc                                  >  2:1470621/1‑150 (MQ=255)
                                                                                                                     tGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGc                                  <  2:95838/150‑1 (MQ=255)
                                                                                                                     tGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGc                                  <  1:577143/150‑1 (MQ=255)
                                                                                                                       aaTCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAATTTTTGGCgg                                >  2:761496/1‑150 (MQ=255)
                                                                                                                        aTCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGt                               <  1:1284027/150‑1 (MQ=255)
                                                                                                                        aTCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGt                               <  1:1277209/150‑1 (MQ=255)
                                                                                                                              ggTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCTGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACg                         <  1:977430/150‑1 (MQ=255)
                                                                                                                                  gTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCtt                     <  1:1233266/150‑1 (MQ=255)
                                                                                                                                  gTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TATTTTCTCCAGTTCTTCTTTAAAGTTGACACTCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGGTAAACGTCtt                     >  2:120919/1‑150 (MQ=255)
                                                                                                                                    tCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCtttt                   <  2:1121925/150‑1 (MQ=255)
                                                                                                                                         cccccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAc              >  1:1097715/1‑150 (MQ=255)
                                                                                                                                          ccccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAcc             >  2:319089/1‑150 (MQ=255)
                                                                                                                                           cccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAccc            >  1:1260849/1‑150 (MQ=255)
                                                                                                                                           cccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAccc            >  2:812534/1‑150 (MQ=255)
                                                                                                                                                   atcatGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTaa    <  1:1007693/150‑1 (MQ=255)
                                                                                                                                                     catGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TATTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTAAAg  >  2:267690/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TGAACCTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACC‑TTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTAAAG  >  NC_000913/1984963‑1985261

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