Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,984,895 G→T 5.1% P233H (CCT→CAT)  yqeF ← short chain acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,8950GT 5.1% 270.5 / 2.3 79P233H (CCT→CAT) yqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (41/34);  new base T (2/2);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

CGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTC  >  NC_000913/2984764‑2985033
                                                                                                                                   |                                                                                                                                          
cgcgAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTaaa                                                                                                                                    <  1:1029256/140‑1 (MQ=255)
cgcgAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTaaa                                                                                                                                    <  2:880749/140‑1 (MQ=255)
cgcgAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTaaa                                                                                                                                    <  2:135188/140‑1 (MQ=255)
cgcgAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTaaa                                                                                                                                    <  2:135148/140‑1 (MQ=255)
cgcgAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTaaa                                                                                                                                    <  2:130761/140‑1 (MQ=255)
cgcgAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTaaa                                                                                                                                    <  1:1832539/140‑1 (MQ=255)
cgcgAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTaaa                                                                                                                                    <  1:1299239/140‑1 (MQ=255)
      tGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCt                                                                                                                              <  1:1197277/140‑1 (MQ=255)
        cGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTaa                                                                                                                            <  1:1243232/140‑1 (MQ=255)
             cAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCtt                                                                                                                       <  2:1569827/140‑1 (MQ=255)
                   gggTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTgcgc                                                                                                                 <  2:742281/140‑1 (MQ=255)
                     gTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTg                                                                                                               >  1:1311985/1‑140 (MQ=255)
                     gTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTg                                                                                                               >  1:1366984/1‑140 (MQ=255)
                     gTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTg                                                                                                               >  2:1902171/1‑140 (MQ=255)
                         aTTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGt                                                                                                           <  1:1138125/140‑1 (MQ=255)
                         aTTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGt                                                                                                           <  2:790876/140‑1 (MQ=255)
                         aTTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGt                                                                                                           <  2:1879064/140‑1 (MQ=255)
                         aTTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGt                                                                                                           <  1:560080/140‑1 (MQ=255)
                         aTTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGt                                                                                                           <  1:184884/140‑1 (MQ=255)
                          ttCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTc                                                                                                          >  2:1382580/1‑140 (MQ=255)
                          ttCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTc                                                                                                          >  1:230817/1‑140 (MQ=255)
                          ttCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTc                                                                                                          >  2:1170334/1‑140 (MQ=255)
                               cGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGc                                                                                                     >  1:433444/1‑140 (MQ=255)
                               cGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGc                                                                                                     >  2:1202487/1‑140 (MQ=255)
                               cGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGc                                                                                                     >  2:1711510/1‑140 (MQ=255)
                               cGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTACTGCGCTGGCGTCAGTGc                                                                                                     >  1:433428/1‑140 (MQ=255)
                               cGCTCGTGCTGTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCAACGTGTATGGATGATGCAATACCCGCTGTCGCAGAACCGAGGCTAACAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCAGGCGGCAGTGc                                                                                                     >  2:476962/1‑140 (MQ=255)
                                 cTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGt                                                                                                   <  2:123195/140‑1 (MQ=255)
                                          tGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAATGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCa                                                                                          >  2:1471017/1‑140 (MQ=255)
                                          tGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAATGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCa                                                                                          >  2:1471027/1‑140 (MQ=255)
                                             cTTCGCTCATCAACATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCg                                                                                       <  2:1902613/140‑1 (MQ=255)
                                                cGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTa                                                                                    >  2:1276737/1‑140 (MQ=255)
                                                cGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTa                                                                                    >  2:1276909/1‑140 (MQ=255)
                                                cGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTa                                                                                    >  1:1809964/1‑140 (MQ=255)
                                                cGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAAAGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTa                                                                                    >  1:1705675/1‑140 (MQ=255)
                                                  cTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATc                                                                                  <  2:648500/140‑1 (MQ=255)
                                                    catcatcatTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAATGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATcaa                                                                                <  1:1471027/140‑1 (MQ=255)
                                                    catcatcatTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAATGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATcaa                                                                                <  1:1471017/140‑1 (MQ=255)
                                                    catcatcatTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATcaa                                                                                <  1:637777/140‑1 (MQ=255)
                                                      tcatcatTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATcaaca                                                                              >  1:1557940/1‑140 (MQ=255)
                                                      tcatcatTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATcaaca                                                                              >  1:1557375/1‑140 (MQ=255)
                                                               ccGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTc                                                                     >  2:704384/1‑140 (MQ=255)
                                                                cGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCt                                                                    <  2:1788212/140‑1 (MQ=255)
                                                                 gCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTg                                                                   <  1:716941/140‑1 (MQ=255)
                                                                  cAGCTGCGCCATCGGTTATGGAGGATGCAGTACCCGCTGTCACAGAACCGAGACGATCAACTGAAGGTGTTAAACGGGCCTAGCCTGCTGCGCAGCCGGCAGTGCGTGGCTGTGCATCGGGATCTACAGCCAACCGCTGc                                                                  >  2:1673208/1‑140 (MQ=255)
                                                                       gcgcCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGtt                                                             >  1:1427250/1‑140 (MQ=255)
                                                                          ccATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACt                                                          <  2:631642/140‑1 (MQ=255)
                                                                                tttATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCTTCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGt                                                    <  2:1100287/140‑1 (MQ=255)
                                                                                tttATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGt                                                    <  2:232948/140‑1 (MQ=255)
                                                                                tttATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGt                                                    <  2:1100410/140‑1 (MQ=255)
                                                                                tttATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGt                                                    <  2:1100300/140‑1 (MQ=255)
                                                                                tttATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGt                                                    <  1:1309710/140‑1 (MQ=255)
                                                                                       atgatgCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTAcc                                             >  1:340514/1‑140 (MQ=255)
                                                                                          atgCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCggg                                          >  1:1074908/1‑140 (MQ=255)
                                                                                          atgCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCggg                                          >  2:125729/1‑140 (MQ=255)
                                                                                          atgCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCggg                                          >  2:1523172/1‑140 (MQ=255)
                                                                                          atgCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCggg                                          >  2:1758092/1‑140 (MQ=255)
                                                                                          atgCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGACCGTTACTTTGGGTCATTACCggg                                          >  2:1229016/1‑140 (MQ=255)
                                                                                           tgCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGa                                         <  1:1168379/140‑1 (MQ=255)
                                                                                                 aCCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtct                                   >  2:397694/1‑140 (MQ=255)
                                                                                                 aCCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtct                                   <  2:365912/140‑1 (MQ=255)
                                                                                                 aCCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtct                                   <  2:1454135/140‑1 (MQ=255)
                                                                                                 aCCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtct                                   >  2:1511564/1‑140 (MQ=255)
                                                                                                 aCCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGAtct                                   <  2:1451310/140‑1 (MQ=255)
                                                                                                 aCCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCATTCTGCGCTGGCGACAGTTCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGACCGTTACTTTGTGTCATAACCTGGAAGAtgt                                   >  2:1918427/1‑138 (MQ=255)
                                                                                                     gCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATc                               >  2:1323343/1‑140 (MQ=255)
                                                                                                     gCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATc                               >  1:1425863/1‑140 (MQ=255)
                                                                                                       tGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTaa                                                                                                                            <  1:1119214/45‑1 (MQ=255)
                                                                                                       tGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTaa                                                                                                                            >  2:1119214/1‑45 (MQ=255)
                                                                                                       tGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCtt                             >  1:1787766/1‑140 (MQ=255)
                                                                                                       tGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCtt                             >  2:671569/1‑140 (MQ=255)
                                                                                                                 ccGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTcc                   >  2:656494/1‑140 (MQ=255)
                                                                                                                    agaCTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGc                >  1:1012974/1‑140 (MQ=255)
                                                                                                                         aTCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCaa           >  2:1115738/1‑140 (MQ=255)
                                                                                                                         aTCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCaa           >  2:745377/1‑140 (MQ=255)
                                                                                                                         aTCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCaa           >  2:1090430/1‑140 (MQ=255)
                                                                                                                                  aGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTc  >  2:574561/1‑140 (MQ=255)
                                                                                                                                  aGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTc  >  2:1633088/1‑140 (MQ=255)
                                                                                                                                  aGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTc  <  2:234002/140‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                          
CGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCGGTATCAACAACCAACGTCTGCCCGTTACTTTGGGTCATTACCGGGACGATCTCATCTTTAAATCGTCCGGCGTCAATCGCCGCTC  >  NC_000913/2984764‑2985033

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: