Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,075,518 A→C 8.9% L1122R (CTC→CGC)  putA ← fused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,075,5180AC 8.9% 161.2 / 2.7 67L1122R (CTC→CGC) putAfused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (34/26);  new base C (4/3);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTA  >  NC_000913/1075372‑1075661
                                                                                                                                                  |                                                                                                                                               
ttCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGCCGCGCGAGc                                                                                                                                               >  1:416928/1‑149 (MQ=255)
      cctgcGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACt                                                                                                                                         >  2:367684/1‑149 (MQ=255)
       ctgcGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCACTGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTg                                                                                                                                        >  1:320802/1‑149 (MQ=255)
        tgcGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCTCGCGCGTCaggat                                                                                                                                       <  2:651784/149‑5 (MQ=255)
        gtcgcCAGCGCGCCATATTGCGTACATACCGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTATCTGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCTCGCGCGTCaggat                                                                                                                                       <  2:487994/147‑5 (MQ=255)
        cgcgcCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCGCGACGATTTGCTGCCCATTCCCGCAGTGCAGTTAGCGGCTGAGTCACTGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCTCGCGCGTCaggat                                                                                                                                       <  1:777589/148‑5 (MQ=255)
                    ccATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCTTCGACCGGATACTTTGCATCCTGACGCGCGCGCGTCACTGCAAGCGCACTTTc                                                                                                                           >  2:834462/1‑149 (MQ=255)
                       tatTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCCTTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGTCCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCgg                                                                                                                        <  2:102331/149‑1 (MQ=255)
                        atTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCTATGCGGCTTTCAACTGCGCAGCGTCCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCggg                                                                                                                       <  1:474480/149‑1 (MQ=255)
                        atTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCggg                                                                                                                       <  1:823351/149‑1 (MQ=255)
                              tACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCTCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTc                                                                                                                 >  2:118190/1‑149 (MQ=255)
                               aCATAACGCCTGCAATTCTGGACGCTTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCg                                                                                                                <  2:678968/149‑1 (MQ=255)
                                cATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGc                                                                                                               >  1:392753/1‑149 (MQ=255)
                                  tAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGAGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCca                                                                                                             <  2:533304/149‑1 (MQ=255)
                                   aaCGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCcag                                                                                                            <  1:224304/149‑1 (MQ=255)
                                         tGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGCCCGGATACTTTGCATCCTGCCGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACg                                                                                                      >  1:275249/1‑149 (MQ=255)
                                              ttCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTCCGCATCGCCCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGAACTTACCGGGCGATTCGCCAGCAGACGGTaga                                                                                                 >  1:653620/1‑149 (MQ=255)
                                                 tGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGAt                                                                                              <  2:276606/149‑1 (MQ=255)
                                                   gACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGTATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAACGCACTTTCCGGGCGATTCGACAGCAGCAGGTAGAGATag                                                                                            >  1:668224/1‑149 (MQ=255)
                                                   gACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATag                                                                                            <  1:367684/149‑1 (MQ=255)
                                                      gATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGc                                                                                         <  2:23876/149‑1 (MQ=255)
                                                      gATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGGCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGc                                                                                         <  1:137438/149‑1 (MQ=255)
                                                          tgctgcCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGAc                                                                                     <  2:392753/149‑1 (MQ=255)
                                                             tgcCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCCTCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGTGCGATTCGCCAGCACACGGAAGAGATAGAGCGGAACGc                                                                                  >  1:319449/1‑149 (MQ=255)
                                                             tgcCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGc                                                                                  <  2:394508/149‑1 (MQ=255)
                                                             tgcCCATTCCCGCAGGGCAGTTCGCGGCTGAGTCAATGCGGCTTTCAACTGCTCATCGACCGGAGACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTGTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGc                                                                                  <  1:62421/149‑1 (MQ=255)
                                                                     cccGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTc                                                                          <  2:547156/149‑1 (MQ=255)
                                                                      ccGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCg                                                                         >  2:885047/1‑149 (MQ=255)
                                                                      ccGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACATTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCg                                                                         >  2:135387/1‑149 (MQ=255)
                                                                       cGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCTGATACTTTGCATCCTGACTCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTATCgg                                                                        >  1:915533/1‑149 (MQ=255)
                                                                       cGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCgg                                                                        >  1:284876/1‑149 (MQ=255)
                                                                         cAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCTAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGcc                                                                      <  2:294532/149‑1 (MQ=255)
                                                                         cAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGAATACTGACGCGCGAGCGCCACTGACAGCGCACTTTCCGGGCGATTCGCCAGCAGTCGGTAGAGAAAGAGAGGACCGCCTGCTTTCGGcc                                                                      >  1:926910/1‑149 (MQ=255)
                                                                           gTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGCCCGCCTGCTTTCGGCCCg                                                                    >  2:952604/1‑149 (MQ=255)
                                                                             gCATTTAGCGGCTGAGTCAATGCGGCGTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGt                                                                  <  1:71609/149‑1 (MQ=255)
                                                                               aTTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGCCCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTAc                                                                >  2:831898/1‑149 (MQ=255)
                                                                               aTTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACTAGCGAGCGTCACTGCCAGAGCACTTTCCTGGCGATTCGCCAGAAGACGGTAGAGATATAGCGGACCGCCTGCTATCGGCCCGGTAc                                                                >  2:246849/1‑149 (MQ=255)
                                                                                   aGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGCCGCGCGAGCCTCACTGGCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCTCCTGCTTTCGGCCCGGTAACGGa                                                            >  2:226510/1‑149 (MQ=255)
                                                                                    gCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGGAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGAc                                                           >  1:532940/1‑149 (MQ=255)
                                                                                         tGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCCGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAAccc                                                      >  2:742125/1‑149 (MQ=255)
                                                                                            gTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGGCTGCTTTCGGCCCCGGACCGGACAGCCCTTc                                                   >  1:818654/1‑149 (MQ=255)
                                                                                               aaTGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGCCGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTcgcc                                                >  1:606332/1‑149 (MQ=255)
                                                                                               aaTGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTcgcc                                                >  2:270462/1‑149 (MQ=255)
                                                                                               aaTGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGAAAGAGCGGCCCGCCTGCTGTCGGCCCGGTCCCGGACAACCCTTcgcc                                                >  2:687624/1‑149 (MQ=255)
                                                                                                    ggCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGCTAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGCCCCCCCTTCGccgccgc                                           >  1:264179/1‑148 (MQ=255)
                                                                                                       tttCAACTGCGCATCGACCGGATAGTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGGCAACCCTTCGCCGCCGAAcc                                        >  2:4982/1‑148 (MQ=255)
                                                                                                          cAACTGCGCATCGACCGGATCCTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGCGCGGACCGCCTGCGTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCt                                     <  2:277590/149‑1 (MQ=255)
                                                                                                           aaCTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTg                                    <  2:301045/149‑1 (MQ=255)
                                                                                                           aaCTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTg                                    >  1:981712/1‑149 (MQ=255)
                                                                                                           aaCTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTATCGGCCCGGAACCGGAACACCCTTCGCCGCCCAACGGCTg                                    >  2:276781/1‑149 (MQ=255)
                                                                                                           aaCTGAGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCCGCAGACGGTAGAGATAGAGAGGACCGCCTGCTTTCGGCCCGGTACCGTACAACCCTTCGCCGCCGACCGGcgg                                    >  1:522066/1‑147 (MQ=255)
                                                                                                                 gcATCGACCGTATACTTTGCATCCTGTCGCGCGAGCGTGACTGCCAGCGCACTTTCCGGGCGATTCTCCAGCAGCCGGTAGAGATAGAGCGGACCGCCTGCGTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACacc                              <  1:799967/149‑1 (MQ=255)
                                                                                                                 gcATCGACCGGATACTTTGCATCCTGACGCGCGCGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGCCCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGCACGGCTGctcccc                              >  2:994175/1‑144 (MQ=255)
                                                                                                                    tCGACCGGATACTTTGCAGCATGACGCGCGAGCGCCACTTCCAGCGCAGTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCGGCACaccaac                           <  2:926910/149‑1 (MQ=255)
                                                                                                                      gACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACaccaacca                         <  2:743235/149‑1 (MQ=255)
                                                                                                                           gATCCTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTgcg                    <  2:1030712/149‑1 (MQ=255)
                                                                                                                              aCTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGccc                 <  1:118190/149‑1 (MQ=255)
                                                                                                                               cTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACCCCAACCACTGCGcccc                >  1:374872/1‑148 (MQ=255)
                                                                                                                                 ttGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCccacc              >  1:823309/1‑149 (MQ=255)
                                                                                                                                   gCATCCTGACGCGCGCGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCCCACCACCCACTGCGCccacgct            >  2:575100/1‑146 (MQ=255)
                                                                                                                                   gCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTACCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCat            <  2:1017862/149‑1 (MQ=255)
                                                                                                                                     gcagtGACGCTCGCGCGTCCCTGCCAGCGCACTTTCCGGGCGATTCGCCCGCAGACGGTAGAGATAGAGCGGCCCGCCTGCTTTCGGCCCCGTACCGGCCAACCCTTCGCCGCCGAACCGCTGCACACCAACCACTGCGCCCACCatat          >  2:109771/5‑149 (MQ=255)
                                                                                                                                       ccTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGCGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCGACCATATTa        >  2:820177/1‑149 (MQ=255)
                                                                                                                                         tGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACg      <  2:416928/149‑1 (MQ=255)
                                                                                                                                          gACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACgg     >  2:786021/1‑149 (MQ=255)
                                                                                                                                           aCGCGCGAGCGTCACTGCCAGCGCACTTTCAGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGTTTTCTGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGt    <  2:319449/149‑1 (MQ=255)
                                                                                                                                             gcgcgAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTa  >  1:507277/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTGACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTA  >  NC_000913/1075372‑1075661

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: ‑