Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,633,090 A→T 7.8% K216M (AAG→ATG) ‡ putA → proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,633,0900AT 7.8% 593.9 / 31.1 192K216M (AAG→ATG) ‡putAproline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (90/86);  new base T (7/8);  total (98/94)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

Only 100 of 197 total aligned reads displayed.

GGCCGAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCATTGCCGAAG  >  NC_002947/5632941‑5633239
                                                                                                                                                     |                                                                                                                                                     
 gCCGAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGACGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAg                                                                                                                                                      <  2:1784807/150‑1 (MQ=255)
  ccGAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGc                                                                                                                                                     <  2:5447910/150‑1 (MQ=255)
   cGAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCt                                                                                                                                                    >  2:4673970/1‑150 (MQ=255)
    gAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTg                                                                                                                                                   <  1:3365299/150‑1 (MQ=255)
      aGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGt                                                                                                                                                 <  1:3967655/150‑1 (MQ=255)
      aGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGt                                                                                                                                                 <  1:4072161/150‑1 (MQ=255)
      aGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGt                                                                                                                                                 <  2:4733033/150‑1 (MQ=255)
      aGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGt                                                                                                                                                 <  1:4072547/150‑1 (MQ=255)
      aGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGt                                                                                                                                                 <  1:2281745/150‑1 (MQ=255)
       gCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTc                                                                                                                                                >  1:2855407/1‑150 (MQ=255)
        cccTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCa                                                                                                                                               <  1:3842334/150‑1 (MQ=255)
        cccTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCa                                                                                                                                               <  2:1256804/150‑1 (MQ=255)
        cccTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCa                                                                                                                                               <  2:2187453/150‑1 (MQ=255)
          ctgctgCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGc                                                                                                                                             <  1:3458352/150‑1 (MQ=255)
          ctgctgCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCAGCTTGCCGGTCAGc                                                                                                                                             <  2:935329/150‑1 (MQ=255)
           tgctgcGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCa                                                                                                                                            <  2:2502500/150‑1 (MQ=255)
           tgctgcGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCAGCTTGCCGGTCAGCa                                                                                                                                            <  1:1927085/150‑1 (MQ=255)
           tgctgcGTATCCCCGACAAAGGGACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCAGCTTGCCGGTCAGCa                                                                                                                                            <  2:1828144/150‑1 (MQ=255)
            gctgcGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGGTGCTGACCAGCTTGCCGGTCAGCag                                                                                                                                           <  1:311438/150‑2 (MQ=255)
            gctgcGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGGTGCTGACCAGCTTGCCGGTCAGCag                                                                                                                                           <  1:4528002/150‑2 (MQ=255)
            gctgcGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCAc                                                                                                                                           <  1:5277936/150‑1 (MQ=255)
            gctgcGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCAGCTTGCCGGTCAGCag                                                                                                                                           <  1:2252329/150‑2 (MQ=255)
              tgcGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCAccc                                                                                                                                         <  1:121971/150‑1 (MQ=255)
                cGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCAGCTTGCCGGTCAGCag                                                                                                                                           >  2:4369126/1‑145 (MQ=255)
                 gTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGACAAGCTGGTCAGCACCCATa                                                                                                                                      >  1:1410247/1‑150 (MQ=255)
                    tCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACg                                                                                                                                   >  1:4614385/1‑150 (MQ=255)
                     ccccGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGa                                                                                                                                  >  2:636041/1‑150 (MQ=255)
                       ccGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGaaa                                                                                                                                >  2:122550/1‑150 (MQ=255)
                       ccGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGaaa                                                                                                                                >  1:2357519/1‑150 (MQ=255)
                           cAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCgg                                                                                                                            <  1:4637025/150‑1 (MQ=255)
                            aaaGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGc                                                                                                                           >  2:4855479/1‑150 (MQ=255)
                            aaaGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGc                                                                                                                           >  1:5108137/1‑150 (MQ=255)
                               ggTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTc                                                                                                                        <  1:125626/150‑1 (MQ=255)
                                   cccGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCAcct                                                                                                                    <  1:1148297/150‑1 (MQ=255)
                                    ccGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCAcctc                                                                                                                   >  2:841187/1‑150 (MQ=255)
                                     cGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCAcctcc                                                                                                                  >  2:2628191/1‑150 (MQ=255)
                                             gcTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCAc                                                                                                          <  2:4853044/150‑1 (MQ=255)
                                              cTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCAcc                                                                                                         <  2:2238185/150‑1 (MQ=255)
                                               tGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCAccc                                                                                                        <  1:2519207/150‑1 (MQ=255)
                                                gATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCg                                                                                                       <  1:2486146/150‑1 (MQ=255)
                                                gATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCg                                                                                                       <  2:3152294/150‑1 (MQ=255)
                                                    cgcgACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGcatc                                                                                                   <  1:622763/150‑1 (MQ=255)
                                                     gcgACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGcatca                                                                                                  >  1:3781223/1‑150 (MQ=255)
                                                     gcgACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGcatca                                                                                                  >  2:5033280/1‑150 (MQ=255)
                                                     gcgACAAGATCAGTACCGGCAACTGGCAGCCGCACCGGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAAc                                                                                                                                    <  2:3332303/116‑1 (MQ=255)
                                                            gATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCaa                                                                                           <  1:1021287/150‑1 (MQ=255)
                                                            gATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCaa                                                                                           >  1:1476239/1‑150 (MQ=255)
                                                                 gTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAgcg                                                                                      <  1:4313464/150‑1 (MQ=255)
                                                                  tACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAgcgg                                                                                     <  2:45856/150‑1 (MQ=255)
                                                                    ccGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAgcggcg                                                                                   >  2:4494937/1‑150 (MQ=255)
                                                                    ccGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAgcggcg                                                                                   >  1:3323054/1‑150 (MQ=255)
                                                                     cGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGa                                                                                  >  1:1999289/1‑150 (MQ=255)
                                                                      cggcaCTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAg                                                                                 >  2:1661289/5‑150 (MQ=255)
                                                                            tGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCgatg                                                                           >  2:1983327/1‑150 (MQ=255)
                                                                              gCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCgatgat                                                                         <  1:4756716/150‑1 (MQ=255)
                                                                              gCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCgatgat                                                                         >  1:2885745/1‑150 (MQ=255)
                                                                                  ccGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGc                                                                     >  2:4646792/1‑150 (MQ=255)
                                                                                   cGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCaa                                                                                           >  1:3279660/1‑127 (MQ=255)
                                                                                   cGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCaa                                                                                           <  2:3279660/127‑1 (MQ=255)
                                                                                   cGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCa                                                                    >  1:5222946/1‑150 (MQ=255)
                                                                                      aCCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGc                                                                                                                           <  1:3201508/92‑1 (MQ=255)
                                                                                      aCCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCAGCGGCACGAGCGGCGAGCCGATGATCCGCAAgg                                                                 >  2:1418719/1‑150 (MQ=255)
                                                                                         tGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCg                                                              >  1:3565879/1‑150 (MQ=255)
                                                                                           ggCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTc                                                            <  1:3055303/150‑1 (MQ=255)
                                                                                             cAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGa                                                          >  2:5195345/1‑150 (MQ=255)
                                                                                             cAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGa                                                          <  2:1069517/150‑1 (MQ=255)
                                                                                              aaCAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGAc                                                         >  2:2664354/1‑150 (MQ=255)
                                                                                                  gCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGGTCCGCAAGGGCGTCGACATgg                                                     <  2:2357519/150‑1 (MQ=255)
                                                                                                  gCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATgg                                                     <  2:4048970/150‑1 (MQ=255)
                                                                                                      gTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCGGCGGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCCCCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCCAGGGCGTCGACATGGCCAt                                                 <  1:43753/150‑1 (MQ=255)
                                                                                                        cACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATgc                                               <  2:3204234/150‑1 (MQ=255)
                                                                                                           tGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCatcatc                                                                                                >  2:3058529/1‑98 (MQ=255)
                                                                                                           tGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGcc                                            >  1:351451/1‑150 (MQ=255)
                                                                                                           tGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGcc                                            <  2:1715869/150‑1 (MQ=255)
                                                                                                             ttCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTg                                          >  1:4171922/1‑150 (MQ=255)
                                                                                                                gTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATg                                       >  2:2297970/1‑150 (MQ=255)
                                                                                                                gTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATg                                       >  2:2845234/1‑150 (MQ=255)
                                                                                                                gTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGACTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATg                                       >  2:295415/1‑150 (MQ=255)
                                                                                                                 tCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCTGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATgg                                      <  2:736494/150‑1 (MQ=255)
                                                                                                                 tCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATgg                                      >  2:2755049/1‑150 (MQ=255)
                                                                                                                 tCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATgg                                      <  2:4546082/150‑1 (MQ=255)
                                                                                                                   aaCGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGt                                    <  2:4232833/150‑1 (MQ=255)
                                                                                                                       ccgccACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGc                                >  1:5331978/1‑150 (MQ=255)
                                                                                                                        cgccACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCa                               <  2:3249432/150‑1 (MQ=255)
                                                                                                                         gccACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCAc                                                                                                                      >  1:4681034/1‑62 (MQ=255)
                                                                                                                           cACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGtt                            >  2:3883201/1‑150 (MQ=255)
                                                                                                                               tGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTc                        <  1:5020627/150‑1 (MQ=255)
                                                                                                                               tGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTc                        <  2:3861321/150‑1 (MQ=255)
                                                                                                                                    atgggtgctgACCAGCTTGCCGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCgg                   >  1:4331571/6‑150 (MQ=255)
                                                                                                                                      gggtgctgACCAGCTTGCCGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCg                 >  1:1007696/4‑150 (MQ=255)
                                                                                                                                       gctgctgACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCga                <  2:5115690/150‑1 (MQ=255)
                                                                                                                                       gctgctgACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCga                >  2:2113080/1‑150 (MQ=255)
                                                                                                                                        ctgctgACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCgag               >  2:1587618/1‑150 (MQ=255)
                                                                                                                                         tgctgACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCgaga              >  1:154005/1‑150 (MQ=255)
                                                                                                                                         tgctgACCAGCTTGCCGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCgaga              >  1:1719322/1‑150 (MQ=255)
                                                                                                                                            tgACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCa           >  1:663431/1‑150 (MQ=255)
                                                                                                                                            tgACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCACAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCa           >  2:4030601/1‑150 (MQ=255)
                                                                                                                                             gACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCAt          >  2:1281475/1‑150 (MQ=255)
                                                                                                                                                cGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCATTGc       >  1:582154/1‑150 (MQ=255)
                                                                                                                                                     aGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCATTGCCGAAg  <  2:606605/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCCGAAGCCCTGCTGCGTATCCCCGACAAAGGTACCCGTGATGCGCTGATCCGCGACAAGATCAGTACCGGCAACTGGCAGCCGCACCTGGGCAACAGCCCGTCACTGTTCGTCAACGCCGCCACCTGGGGGCTGCTGCTGACCGGCAAGCTGGTCAGCACCCATAACGAAACCGGCCTCACCTCCTCGCTCACCCGCATCATCGGCAAGAGCGGCGAGCCGATGATCCGCAAGGGCGTCGACATGGCCATGCGCCTGATGGGTGAGCAGTTCGTCACCGGCGAGACCATTGCCGAAG  >  NC_002947/5632941‑5633239

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