Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,915,455 G→C 7.4% A27G (GCG→GGG)  proQ ← RNA chaperone ProQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,915,4550GC 7.4% 283.8 / 9.3 92A27G (GCG→GGG) proQRNA chaperone ProQ
Reads supporting (aligned to +/- strand):  ref base G (38/44);  new base C (2/5);  total (41/51)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

CGGGTTTAACACCGTAAAGATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTA  >  NC_000913/1915306‑1915604
                                                                                                                                                     |                                                                                                                                                     
cGGGTTTAACACCGTAAAGATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCg                                                                                                                                                       <  2:1025165/150‑1 (MQ=255)
 gggTTTAACACCGTAAAGATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGc                                                                                                                                                      <  2:409195/150‑1 (MQ=255)
  ggTTTAACACCGTAAAGATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCa                                                                                                                                                     <  2:574025/150‑1 (MQ=255)
   gTTTAACACCGTAAAGATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACCCGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCCCCTTCCGCAc                                                                                                                                                    >  2:818492/1‑150 (MQ=255)
      tAACACCGTAAAGATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCCCACTgt                                                                                                                                                 >  2:468483/1‑149 (MQ=255)
       aaCACCGTAAAGATACCGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCCGCAACACGATCGCCCAAATCCGGAAAAATACCGATTTTCCGCGGCGGGGCTTCACCTTCCGCCCGGaa                                                                                                                                                >  2:868705/1‑150 (MQ=255)
       aaCACCGTAAAGATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGaa                                                                                                                                                <  1:616276/150‑1 (MQ=255)
       aaCACCGTAAAGATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCCCACTgtt                                                                                                                                                >  2:1449993/1‑148 (MQ=255)
        acacCGTAAAGATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGaaa                                                                                                                                               >  1:373124/1‑150 (MQ=255)
        acacCGTAAAGATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCACACTGaaa                                                                                                                                               >  1:1310445/1‑150 (MQ=255)
         cacCGTAAAGATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAAc                                                                                                                                              <  2:40427/150‑1 (MQ=255)
         cacCGTAAAGATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAAc                                                                                                                                              <  2:1616136/150‑1 (MQ=255)
              tAAAGATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTgg                                                                                                                                         <  1:1607756/150‑1 (MQ=255)
              tAAAGATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTgg                                                                                                                                         <  1:53527/150‑1 (MQ=255)
               aaaGATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTggg                                                                                                                                        <  1:656181/150‑1 (MQ=255)
                 aGATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGa                                                                                                                                      <  2:590319/150‑1 (MQ=255)
                  gATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCCCACTGTTTCAGTGCGGaa                                                                                                                                     <  1:1449993/150‑1 (MQ=255)
                  gATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCCCACTGTTTCAGTGCGGaa                                                                                                                                     <  1:468483/150‑1 (MQ=255)
                   aTAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCCCACTGTTTCAGTGCGGaag                                                                                                                                    <  2:224883/150‑2 (MQ=255)
                   aTAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCCCACTGTTTCAGTGCGGaag                                                                                                                                    <  2:1353467/150‑2 (MQ=255)
                    tAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGaaaa                                                                                                                                   >  2:1118128/1‑150 (MQ=255)
                    tAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGaaaa                                                                                                                                   >  2:666976/1‑150 (MQ=255)
                    tAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCCCACTGTTTCAGTGCGGaagg                                                                                                                                   <  1:83103/150‑3 (MQ=255)
                    tAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCACCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGaaaa                                                                                                                                   >  2:602763/1‑150 (MQ=255)
                      aCGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACg                                                                                                                                 <  1:623564/150‑1 (MQ=255)
                       cGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGt                                                                                                                                >  1:548744/1‑150 (MQ=255)
                          cAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCg                                                                                                                             <  1:1273782/150‑1 (MQ=255)
                           aGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCgg                                                                                                                            <  1:116173/150‑1 (MQ=255)
                           aGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTGCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCgg                                                                                                                            <  1:238302/150‑1 (MQ=255)
                             cTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGcc                                                                                                                          <  1:784967/150‑1 (MQ=255)
                                gAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCGCTGAAACAGTGGGGAAAACGTTCGGCCAGa                                                                                                                       >  1:554461/1‑150 (MQ=255)
                                gAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTCAAACAGTGGGGAAAACGTTCGGCAAGa                                                                                                                       >  1:844500/1‑150 (MQ=255)
                                    tgtAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAAcg                                                                                                                   <  1:297057/150‑1 (MQ=255)
                                      tAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAAcgcg                                                                                                                 <  2:468463/150‑1 (MQ=255)
                                            cGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTAct                                                                                                           >  1:361284/1‑150 (MQ=255)
                                             gTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTActt                                                                                                          >  1:188461/1‑150 (MQ=255)
                                               aaaGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGTGGAAAACGTTCGGCCAGAAACGCGATTActtct                                                                                                        >  2:1241119/1‑150 (MQ=255)
                                               aaaGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTActtct                                                                                                        >  2:338486/1‑150 (MQ=255)
                                                aaGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTActtctt                                                                                                       <  1:331610/150‑1 (MQ=255)
                                                  gCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGCCCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTa                                                                                                     >  1:526714/1‑150 (MQ=255)
                                                    ggATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACt                                                                                                   <  2:51740/150‑1 (MQ=255)
                                                      aTCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGc                                                                                                 >  1:1399059/1‑150 (MQ=255)
                                                          cAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTAtt                                                                                             >  2:472316/1‑150 (MQ=255)
                                                                       cTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGa                                                                                >  2:374580/1‑150 (MQ=255)
                                                                         cAGGTTCATTTGCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGAtt                                                                              >  1:448755/1‑150 (MQ=255)
                                                                         cAGGTTCATTTCCCCAGCAACACGATCGCCCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGAtt                                                                              >  2:747931/1‑150 (MQ=255)
                                                                           ggTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGAtttt                                                                            >  1:99307/1‑150 (MQ=255)
                                                                           ggTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGAtttt                                                                            >  1:454607/1‑150 (MQ=255)
                                                                            gTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTc                                                                           <  2:735058/150‑1 (MQ=255)
                                                                            gTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTc                                                                           >  1:301693/1‑150 (MQ=255)
                                                                            gTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACCTAAGTTGaatagc                                                                           >  1:1262479/1‑145 (MQ=255)
                                                                             ttCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACCTAAGTTGaatagca                                                                          <  2:1262479/150‑7 (MQ=255)
                                                                              tCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGGGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCa                                                                         <  1:811699/150‑1 (MQ=255)
                                                                               cATTTACCCAGCAACACGATGGACCAAATACTGAAAAATACCGATTTTAAGCGGACGCGATTACCCTTACGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTAATTCTTTACTGCTATTAAACTTAGGTTGATTTTCCAt                                                                        >  2:901589/1‑150 (MQ=255)
                                                                                     cccAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAAtt                                                                  >  1:1361074/1‑150 (MQ=255)
                                                                                         gCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCt                                                              <  2:735374/150‑1 (MQ=255)
                                                                                         gCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTAATCAACTTAGGTTGATTTTCCATGAAATTTCCt                                                              <  2:741172/150‑1 (MQ=255)
                                                                                           aaCACGATAGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACACTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGa                                                            >  1:49435/1‑150 (MQ=255)
                                                                                            acacGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGAt                                                           <  2:426947/150‑1 (MQ=255)
                                                                                               cGATCGACCAAATCCTGAAAAATACCGATTTTGAGCGGAAGCGCTTCACCTTCCGGACTGAAACTGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTAc                                                        <  2:572086/150‑1 (MQ=255)
                                                                                                    gACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACgg                                                   >  1:676882/1‑150 (MQ=255)
                                                                                                    gACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCACACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACgg                                                   <  2:1310445/150‑1 (MQ=255)
                                                                                                       cAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACg                                                <  2:309233/150‑1 (MQ=255)
                                                                                                              tGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCcaa                                         <  2:933698/150‑1 (MQ=255)
                                                                                                                  aaaTACCGATTTTCAGCGGACGCGCTTCACCTTCCTCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAg                                     <  1:788902/150‑1 (MQ=255)
                                                                                                                    aTACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAgcg                                   >  2:1389048/1‑150 (MQ=255)
                                                                                                                      aCCGATTTTCAGCGGACGCGCTTCCCCTTCCGCACTGAAACAGTGGCCAAAACGTTCGGCCAGAAACGCGATTACTTGTTTACTGCAATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCa                                 >  2:397114/1‑150 (MQ=255)
                                                                                                                      aCCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCa                                 >  2:391235/1‑150 (MQ=255)
                                                                                                                       ccGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAg                                <  1:666976/150‑1 (MQ=255)
                                                                                                                       ccGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAg                                >  1:875317/1‑150 (MQ=255)
                                                                                                                          aTTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCa                             <  1:1118128/150‑1 (MQ=255)
                                                                                                                          aTTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCa                             <  2:301693/150‑1 (MQ=255)
                                                                                                                            tttCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATg                           <  1:598581/150‑1 (MQ=255)
                                                                                                                              tCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGaa                         <  1:338486/150‑1 (MQ=255)
                                                                                                                               cAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAAc                        >  2:1403928/1‑150 (MQ=255)
                                                                                                                                 gCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAg                      >  2:870234/1‑150 (MQ=255)
                                                                                                                                 gCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAg                      >  2:864682/1‑150 (MQ=255)
                                                                                                                                  cGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAgg                     >  1:432091/1‑150 (MQ=255)
                                                                                                                                      cgcgCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTc                 <  1:514639/150‑1 (MQ=255)
                                                                                                                                        cgcTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCAt               <  1:1166065/150‑1 (MQ=255)
                                                                                                                                        cgcTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCAt               <  2:52903/150‑1 (MQ=255)
                                                                                                                                         gcTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCAtt              >  2:1305899/1‑150 (MQ=255)
                                                                                                                                            tCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGCCCAGAAACGCGATTTCTTCTTTACTGCTATTCTATTTATTTTGATTTTGCATGAAATTTCGTGATTACATCGGATGTATCCATCTAGAGCAGTCTTTATTAGGAGTCAttttt           <  2:844500/150‑4 (MQ=255)
                                                                                                                                             cACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATaa          <  1:642818/150‑1 (MQ=255)
                                                                                                                                              aCCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAAt         <  1:828418/150‑1 (MQ=255)
                                                                                                                                                 ttCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGAc      >  2:23068/1‑150 (MQ=255)
                                                                                                                                                  tCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACg     <  2:1531617/150‑1 (MQ=255)
                                                                                                                                                  tCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACg     >  1:1155897/1‑150 (MQ=255)
                                                                                                                                                  tCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACg     <  2:237480/150‑1 (MQ=255)
                                                                                                                                                   ccGCACTGAAACAGTGGGGAAAACGTGCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGc    <  1:868705/150‑1 (MQ=255)
                                                                                                                                                    cGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGcc   <  1:1389048/150‑2 (MQ=255)
                                                                                                                                                    cGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGcc   <  1:1305899/150‑2 (MQ=255)
                                                                                                                                                    cGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGcc   <  1:1284008/150‑2 (MQ=255)
                                                                                                                                                     gCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTa  >  2:164031/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGGTTTAACACCGTAAAGATAACGCCAGCTCGAAGTGTAGAGACGTAAAGCGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTA  >  NC_000913/1915306‑1915604

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