Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,581,762 G→A 7.3% A339V (GCC→GTC)  ypfG ← DUF1176 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,581,7620GA 7.3% 525.6 / 21.3 165A339V (GCC→GTC) ypfGDUF1176 family protein
Reads supporting (aligned to +/- strand):  ref base G (85/68);  new base A (6/6);  total (91/74)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 169 total aligned reads displayed.

TGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTC  >  NC_000913/2581617‑2581910
                                                                                                                                                  |                                                                                                                                                    
tGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGCTGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGACgc                                                                                                                                                   <  3:55402/149‑2 (MQ=255)
    gggAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                                                               <  3:88916/149‑1 (MQ=255)
    gggAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                                                               >  2:836757/1‑149 (MQ=255)
        gTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCa                                                                                                                                           <  1:103436/149‑1 (MQ=255)
           aaaGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGc                                                                                                                                        >  4:406255/1‑149 (MQ=255)
            aaGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGcc                                                                                                                                       >  4:98075/1‑149 (MQ=255)
              ggTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGACGCCTGGCCCacgttgt                                                                                                                                     <  4:1085220/149‑7 (MQ=255)
                tGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGtt                                                                                                                                   <  4:1009262/149‑1 (MQ=255)
                tGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGtt                                                                                                                                   >  1:108695/1‑149 (MQ=255)
                   cTGATCGTGATTGCCGGATGCTGATTGTCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATc                                                                                                                                <  4:1489498/149‑1 (MQ=255)
                    tGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCAGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCg                                                                                                                               >  2:1320462/1‑149 (MQ=255)
                      aTGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCa                                                                                                                             >  1:1672246/1‑149 (MQ=255)
                       tGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAg                                                                                                                            >  1:789687/1‑149 (MQ=255)
                           gattgccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTgg                                                                                                                        <  4:1463512/149‑1 (MQ=255)
                             ttgccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGt                                                                                                                      <  2:1079811/149‑1 (MQ=255)
                                 cggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTg                                                                                                                                                >  1:1566441/1‑119 (MQ=255)
                                   gatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcg                                                                                                                <  4:1240936/149‑1 (MQ=255)
                                       ctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGt                                                                                                            <  1:149783/149‑1 (MQ=255)
                                         gattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTa                                                                                                           >  4:245121/1‑149 (MQ=255)
                                         gattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTaa                                                                                                          >  2:1494908/1‑149 (MQ=255)
                                             gccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCa                                                                                                      <  4:1345775/149‑1 (MQ=255)
                                                 gatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAg                                                                                                  >  2:703193/1‑149 (MQ=255)
                                                      gACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAg                                                                                                  >  3:1512924/1‑144 (MQ=255)
                                                      gACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGaa                                                                                             <  1:934275/149‑1 (MQ=255)
                                                         gCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                                                               >  1:669104/1‑96 (MQ=255)
                                                         gCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                                                               <  2:669104/96‑1 (MQ=255)
                                                           tGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGc                                                                                        <  2:100859/149‑1 (MQ=255)
                                                             aCGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTg                                                                                      >  3:697712/1‑149 (MQ=255)
                                                             aCGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTg                                                                                      <  3:374026/149‑1 (MQ=255)
                                                              cgcgTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCAAAAACGCTGa                                                                                     <  1:1494908/149‑1 (MQ=255)
                                                              cgcgTCTCATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGa                                                                                     <  3:764100/149‑1 (MQ=255)
                                                               gcgTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcgg                                                                                                               >  3:1143947/1‑122 (MQ=255)
                                                               gcgTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcgg                                                                                                               <  4:1143947/122‑1 (MQ=255)
                                                               gcgTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGAc                                                                                    >  2:1186869/1‑149 (MQ=255)
                                                                  tCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCTCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGt                                                                                 <  4:841235/149‑1 (MQ=255)
                                                                  tCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTgg                                                                                                                        >  1:531963/1‑110 (MQ=255)
                                                                   cTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTc                                                                                <  4:201783/149‑1 (MQ=255)
                                                                      aTCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCaaa                                                                             <  2:1640261/149‑1 (MQ=255)
                                                                       tCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAg                                                                            >  2:1048111/1‑149 (MQ=255)
                                                                       tCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAg                                                                            <  1:890796/149‑1 (MQ=255)
                                                                        ccGGCCTACAGCGTCTTATCCCGCCTACATGTCCCGCCA‑TTGTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAgc                                                                           <  4:1624011/149‑1 (MQ=255)
                                                                         cGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAgcg                                                                          >  2:1602140/1‑149 (MQ=255)
                                                                         cGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAgcg                                                                          >  1:333375/1‑149 (MQ=255)
                                                                                   cGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGcc                                                                <  2:341876/149‑1 (MQ=255)
                                                                                    ttCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCt                                                               <  2:1492200/148‑1 (MQ=255)
                                                                                    gTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCt                                                               <  1:81130/149‑1 (MQ=255)
                                                                                     tCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTg                                                              >  3:553883/1‑149 (MQ=255)
                                                                                     tCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTg                                                              >  3:1394802/1‑149 (MQ=255)
                                                                                      cTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACTGTCAAAGCGCCAGCGCGCCTGa                                                             >  3:1058689/1‑149 (MQ=255)
                                                                                       ttATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGaa                                                            >  1:377426/1‑149 (MQ=255)
                                                                                          tCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCa                                                                                                   <  4:1579762/107‑1 (MQ=255)
                                                                                          tCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCa                                                                                                   >  3:1579762/1‑107 (MQ=255)
                                                                                          tCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACATGGGCCAGACGCCTGGCCCacgt                                                                                                                                        <  4:1545444/70‑4 (MQ=37)
                                                                                            cGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCca                                                       >  4:1585560/1‑149 (MQ=255)
                                                                                                cTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAATGCGCCAGCGCGCCTGAATGCCACAAt                                                   <  4:697712/149‑1 (MQ=255)
                                                                                                 tACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATc                                                  >  1:448612/1‑149 (MQ=255)
                                                                                                  aCATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCg                                                 >  3:1249375/1‑149 (MQ=255)
                                                                                                   cATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGc                                                <  3:318541/149‑1 (MQ=255)
                                                                                                   cATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGc                                                <  3:246767/149‑1 (MQ=255)
                                                                                                         ccGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATc                                          >  1:1187726/1‑149 (MQ=255)
                                                                                                         ccGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATc                                          >  2:1310989/1‑149 (MQ=255)
                                                                                                          cGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATcc                                         >  3:1020758/1‑149 (MQ=255)
                                                                                                          cGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATcc                                         >  3:714101/1‑149 (MQ=255)
                                                                                                          cGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATcc                                         >  3:1521756/1‑149 (MQ=255)
                                                                                                          cGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATcc                                         >  2:368897/1‑149 (MQ=255)
                                                                                                          cGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATcc                                         >  1:1296544/1‑149 (MQ=255)
                                                                                                           gCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCcg                                        >  1:1647659/1‑149 (MQ=255)
                                                                                                                ttttGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCt                                                                                       >  4:1431033/1‑98 (MQ=255)
                                                                                                                ttttGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCt                                                                                       <  3:1431033/98‑1 (MQ=255)
                                                                                                                ttttGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGcc                                    >  2:321752/1‑149 (MQ=255)
                                                                                                                ttttGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGcc                                    >  3:1445157/1‑149 (MQ=255)
                                                                                                                ttttGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGcc                                    >  1:13536/1‑149 (MQ=255)
                                                                                                                  ttGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCt                                  <  4:287658/149‑1 (MQ=255)
                                                                                                                     tttAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTc                               <  3:1613388/149‑1 (MQ=255)
                                                                                                                        aaCGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCa                                                                                                                                      >  4:1606252/1‑43 (MQ=255)
                                                                                                                        aaCGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCt                            >  1:1408464/1‑149 (MQ=255)
                                                                                                                        aaCGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCt                            >  3:162569/1‑149 (MQ=255)
                                                                                                                          cGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTaa                          <  1:1004874/149‑1 (MQ=255)
                                                                                                                             gTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGt                       <  4:931991/149‑1 (MQ=255)
                                                                                                                                aTCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCAc                    <  1:1186869/149‑1 (MQ=255)
                                                                                                                                 tCCACAACGTGGGCCAGGCGTCTGGTCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCAcc                   >  3:511474/1‑149 (MQ=255)
                                                                                                                                 tCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCAcc                   >  1:316364/1‑149 (MQ=255)
                                                                                                                                    acaACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGt                >  2:860665/1‑149 (MQ=255)
                                                                                                                                    acaACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGt                >  2:498820/1‑149 (MQ=255)
                                                                                                                                       gcatGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCa             >  4:74425/4‑149 (MQ=255)
                                                                                                                                        cGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAc            <  2:13536/149‑1 (MQ=255)
                                                                                                                                        cGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAc            <  4:656132/149‑1 (MQ=255)
                                                                                                                                         gTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAcg           <  4:374469/149‑1 (MQ=255)
                                                                                                                                         atGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAcg           >  2:1451720/2‑149 (MQ=255)
                                                                                                                                         atGGGCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAcg           >  4:504500/2‑149 (MQ=255)
                                                                                                                                           gggCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGt                                                                                 >  4:826762/1‑77 (MQ=255)
                                                                                                                                            ggCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGc                                                >  2:53408/1‑109 (MQ=255)
                                                                                                                                            ggCCAGACGCCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGa        <  1:1451720/149‑1 (MQ=255)
                                                                                                                                               cAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGAttt     >  1:1601393/1‑149 (MQ=255)
                                                                                                                                                aGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCgcgc                                     <  2:1684135/116‑1 (MQ=255)
                                                                                                                                                aGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTtc    >  2:288319/1‑149 (MQ=255)
                                                                                                                                                  gCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAc            >  4:1366504/1‑139 (MQ=255)
                                                                                                                                                  gCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTtctc  <  1:325554/149‑1 (MQ=255)
                                                                                                                                                  gCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTtctc  >  1:948559/1‑149 (MQ=255)
                                                                                                                                                  gCGTCTGGCCCAGGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCAc            <  3:1366504/139‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCA‑TTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTC  >  NC_000913/2581617‑2581910

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑