Predicted mutation
evidence seq id position mutation annotation gene description
RA AM406671 2,044,751 G→A A37V (GCA→GTA)  ftsA ← cell division protein ftsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM4066712,044,7510GA99.6% 3896.6 / NA 983A37V (GCA→GTA) ftsAcell division protein ftsA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (549/428);  minor base C (3/1);  total (553/430)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 983 total aligned reads displayed.

TTTCGAGGGAATTTGCTGGAATGCTCACGTTAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTGCGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTCCCTCACACATCACGATTCAT  >  AM406671/2044605‑2044896
                                                                                                                                                  |                                                                                                                                                 
  tCGAGGGAATTTGCTGGAATGCTCACGTTAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGtt                                                                                                                                               <  2:9042005/149‑1 (MQ=255)
   cGAGGGAATTTGCTGGAATGCTCACGTTAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGttt                                                                                                                                              >  2:2151270/1‑149 (MQ=255)
     aGGGAATTTGCTGGAATGCTCACGTTAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTcc                                                                                                                                            >  1:6852363/1‑149 (MQ=255)
        gAATTTGCTGGAATGCTCACGTTAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGac                                                                                                                                         >  1:5953690/1‑149 (MQ=255)
        gAATTTGCTGGAATGCTCACGTTAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGac                                                                                                                                         >  2:7149800/1‑149 (MQ=255)
         aaTTTGCTGGAATGCTCACGTTAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGaca                                                                                                                                        >  1:3675441/1‑149 (MQ=255)
              gCTGGAATGCTCACGTTAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAAt                                                                                                                                   >  2:4243555/1‑149 (MQ=255)
                 ggAATGCTCACGTTAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATtat                                                                                                                                >  1:4074817/1‑149 (MQ=255)
                  gAATGCTCACGTTAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATtata                                                                                                                               >  1:3395512/1‑149 (MQ=255)
                    aTGCTCACGTTAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCAATTTTTAAGGCCATTGGATTTTACGTTTCCGACACCAATTATAtt                                                                                                                             >  2:2468469/1‑149 (MQ=255)
                       cTCACGTTAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCAt                                                                                                                          >  2:8448917/1‑149 (MQ=255)
                        tCACGTTAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCAtt                                                                                                                         <  1:7477511/149‑1 (MQ=255)
                          aCGTTAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTc                                                                                                                       >  1:5852008/1‑149 (MQ=255)
                           cGTTAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCa                                                                                                                      >  1:4621937/1‑149 (MQ=255)
                            gTTAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCAc                                                                                                                     <  2:1483985/149‑1 (MQ=255)
                             ttAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCAcc                                                                                                                    >  2:3061791/1‑149 (MQ=255)
                             ttAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCAcc                                                                                                                    >  1:1203221/1‑149 (MQ=255)
                              tAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCt                                                                                                                   >  1:9160599/1‑149 (MQ=255)
                                aTTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTga                                                                                                                 >  1:7750529/1‑149 (MQ=255)
                                aTTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTga                                                                                                                 >  2:4399851/1‑149 (MQ=255)
                                aTTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTga                                                                                                                 >  1:857263/1‑149 (MQ=255)
                                aTTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTga                                                                                                                 >  2:2734627/1‑149 (MQ=255)
                                 ttCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTgag                                                                                                                <  1:6741039/149‑1 (MQ=255)
                                    gATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACg                                                                                                             >  1:4243907/1‑149 (MQ=255)
                                    gATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACg                                                                                                             >  2:7814160/1‑149 (MQ=255)
                                       cAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTAt                                                                                                          >  1:3033981/1‑149 (MQ=255)
                                         aTCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTTATGCTTGTGCTCCCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTTCGTTTCCGACACCAATTATATTCATTTCACCTGAACCGTATTc                                                                                                        >  2:1268644/1‑149 (MQ=255)
                                         aTCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTc                                                                                                        >  1:1716345/1‑149 (MQ=255)
                                          tCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCa                                                                                                       <  1:959593/149‑1 (MQ=255)
                                           cGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAg                                                                                                      <  1:4825159/149‑1 (MQ=255)
                                           cGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAg                                                                                                      >  1:618644/1‑149 (MQ=255)
                                             tGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGca                                                                                                    <  2:2874731/149‑1 (MQ=255)
                                              gATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGcaa                                                                                                   >  2:3967774/1‑149 (MQ=255)
                                               aTTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGAGGTATTCAGcaac                                                                                                  >  2:5456945/1‑149 (MQ=255)
                                                ttCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGcaaca                                                                                                 >  1:4367941/1‑149 (MQ=255)
                                                ttCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGcaaca                                                                                                 >  1:350128/1‑149 (MQ=255)
                                                 tCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCATAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGcaacaa                                                                                                <  1:886437/149‑1 (MQ=255)
                                                    ggCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCa                                                                                             >  1:1894383/1‑149 (MQ=255)
                                                     gcgcGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCAc                                                                                            >  1:3799202/1‑149 (MQ=255)
                                                      cgcgTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCAcc                                                                                           >  2:9154923/1‑149 (MQ=255)
                                                        cgTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCtt                                                                                         <  1:6374274/149‑1 (MQ=255)
                                                         gTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCttt                                                                                        >  1:1022705/1‑149 (MQ=255)
                                                           tcttcAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTAt                                                                                      >  2:6771042/1‑149 (MQ=255)
                                                           tcttcAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTAt                                                                                      >  2:3287407/1‑149 (MQ=255)
                                                           tcttcAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTAt                                                                                      <  2:1129280/149‑1 (MQ=255)
                                                              tcAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATtgt                                                                                   <  2:2013359/149‑1 (MQ=255)
                                                                agcagcATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTgtgt                                                                                 <  2:8897962/149‑1 (MQ=255)
                                                                agcagcATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTgtgt                                                                                 >  2:4119867/1‑149 (MQ=255)
                                                                   agcaTTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTg                                                                              <  1:8119834/149‑1 (MQ=255)
                                                                    gcaTTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGt                                                                             >  1:8976598/1‑149 (MQ=255)
                                                                      aTTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTcc                                                                           <  1:7308561/149‑1 (MQ=255)
                                                                      aTTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTcc                                                                           >  1:3917818/1‑149 (MQ=255)
                                                                       tttACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTccc                                                                          >  2:7939200/1‑149 (MQ=255)
                                                                         tACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGa                                                                        >  2:7260738/1‑149 (MQ=255)
                                                                          aCAGCTTTTCGTAATACTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGat                                                                       <  2:6072649/149‑1 (MQ=255)
                                                                           cAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGata                                                                      >  1:2275391/1‑149 (MQ=255)
                                                                            aGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGatat                                                                     <  1:1581258/149‑1 (MQ=255)
                                                                              cTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCt                                                                   <  2:8162629/149‑1 (MQ=255)
                                                                                  tCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGc                                                               >  1:8387443/1‑149 (MQ=255)
                                                                                  tCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGc                                                               >  1:7147574/1‑149 (MQ=255)
                                                                                   cGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGtt                                                                                <  1:7608112/131‑1 (MQ=255)
                                                                                        tGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGtat                                                         >  1:791487/1‑149 (MQ=255)
                                                                                         gCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGCTTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACCCCTGtata                                                        >  1:545866/1‑149 (MQ=255)
                                                                                             gtgCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGt                                                    <  2:3804803/149‑1 (MQ=255)
                                                                                                 tACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCAc                                                <  1:6976947/149‑1 (MQ=255)
                                                                                                 tACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGGCCCGATATCTACGCCTGTATAAAGTCCAc                                                <  2:6466974/149‑1 (MQ=255)
                                                                                                  aCCTTCTCAATATCTACAATCATTCCATTTTTAAGGCCATCCGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTCTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACt                                               <  2:476500/149‑1 (MQ=255)
                                                                                                   ccTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACtt                                              <  2:761309/149‑1 (MQ=255)
                                                                                                   ccTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACtt                                              >  1:564149/1‑149 (MQ=255)
                                                                                                    cTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACttt                                             >  1:2785744/1‑149 (MQ=255)
                                                                                                        tcAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGc                                         >  2:2330012/1‑149 (MQ=255)
                                                                                                        tcAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGc                                         >  2:4944564/1‑149 (MQ=255)
                                                                                                         cAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTCTCCGACAACAATTATATTCATTTCACCTGCGACGTATTCAGGGACAAGCACCTTTATTGGGTGTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTggac                                        >  2:6628924/1‑147 (MQ=255)
                                                                                                          aaTATCTACAATAATTCCATTTTTTTTTTTTTTTTATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCTCTTTTGGCCa                                       >  2:3648363/1‑149 (MQ=255)
                                                                                                          aaTATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCa                                       >  2:4617849/1‑149 (MQ=255)
                                                                                                           atatCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCAt                                      >  1:1088012/1‑149 (MQ=255)
                                                                                                              tCTACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATaaa                                   >  2:6389878/1‑149 (MQ=255)
                                                                                                                tACAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAg                                 >  2:1845414/1‑149 (MQ=255)
                                                                                                                  cAATAATTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCg                               >  2:743803/1‑149 (MQ=255)
                                                                                                                    ataatTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGAc                             <  1:244749/149‑1 (MQ=255)
                                                                                                                    ataatTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGAc                             >  2:3405873/1‑149 (MQ=255)
                                                                                                                       atTCCATTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACATTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTcc                          <  2:5325467/149‑1 (MQ=255)
                                                                                                                        ttCCATTTTTAAGGCCATCGGATTTAACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCt                         >  2:5056296/1‑149 (MQ=255)
                                                                                                                            aTTTTTAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTcc                     <  1:5209206/149‑1 (MQ=255)
                                                                                                                            aTTTTTAACGCCATCTGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCATTTTTTTTTTTAATGATACGGCGACCCCCGAGATCTACACCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTcc                     <  2:5849181/149‑1 (MQ=255)
                                                                                                                             tttttAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTccc                    >  2:1716924/1‑149 (MQ=255)
                                                                                                                              ttttAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTCCCt                   >  2:85153/1‑149 (MQ=255)
                                                                                                                                ttAAGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTCCCTca                 >  1:3244848/1‑149 (MQ=255)
                                                                                                                                  aaGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTCCCTcaca               >  1:6089515/1‑149 (MQ=255)
                                                                                                                                  aaGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTCCCTcaca               >  1:8676270/1‑149 (MQ=255)
                                                                                                                                   aGGCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTCCCTcacac              >  1:499283/1‑149 (MQ=255)
                                                                                                                                    ggCCATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTCCCTcacaca             <  2:6924901/149‑1 (MQ=255)
                                                                                                                                      ccATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTCCCTCACACATc           <  1:2602822/149‑1 (MQ=255)
                                                                                                                                      ccATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTCCCTCACACATc           <  1:6481540/149‑1 (MQ=255)
                                                                                                                                       cATCGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTCCCTCACACATCa          >  2:6227499/1‑149 (MQ=255)
                                                                                                                                          cGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTCCCTCACACATCACGa       >  1:8787098/1‑149 (MQ=255)
                                                                                                                                          cGGATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTCCCTCACACATCACGa       >  2:568128/1‑149 (MQ=255)
                                                                                                                                            gATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTCCCTCACACATCACGAtt     >  1:4160217/1‑149 (MQ=255)
                                                                                                                                            gATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTCCCTCACACATCACGAtt     <  1:1391752/149‑1 (MQ=255)
                                                                                                                                            gATTTTACGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTCCCTCACACATCACGAtt     <  2:8672907/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTTCGAGGGAATTTGCTGGAATGCTCACGTTAATTCGATCAATCGTGATTCCGGCGCGTTCTTCAGCAGCATTTACAGCTTTTCGTAATGCTTGTGCTACCTTCTCAATATCTACAATAATTCCATTTTTAAGGCCATCGGATTTTGCGTTTCCGACACCAATTATATTCATTTCACCTGAGACGTATTCAGCAACAAGCACCTTTATTGTGTTTGTCCCGATATCTACGCCTGTATAAAGTCCACTTTTGGCCATAAAAGCGACTCCTTTCCCTCACACATCACGATTCAT  >  AM406671/2044605‑2044896

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