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.

TGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTC  >  NC_000913/2581617‑2581910
                                                                                                                                                  |                                                                                                                                                    
tGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGCTGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGACgc                                                                                                                                                   <  3:55402/149‑2 (MQ=255)
    gggAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                                                               <  3:88916/149‑1 (MQ=255)
    gggAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                                                               >  2:836757/1‑149 (MQ=255)
        gTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCa                                                                                                                                           <  1:103436/149‑1 (MQ=255)
           aaaGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGc                                                                                                                                        >  4:406255/1‑149 (MQ=255)
            aaGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGcc                                                                                                                                       >  4:98075/1‑149 (MQ=255)
              ggTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGACGCCTGGCCCacgttgt                                                                                                                                     <  4:1085220/149‑7 (MQ=255)
                tGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGtt                                                                                                                                   <  4:1009262/149‑1 (MQ=255)
                tGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGtt                                                                                                                                   >  1:108695/1‑149 (MQ=255)
                   cTGATCGTGATTGCCGGATGCTGATTGTCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATc                                                                                                                                <  4:1489498/149‑1 (MQ=255)
                    tGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCAGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCg                                                                                                                               >  2:1320462/1‑149 (MQ=255)
                      aTGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCa                                                                                                                             >  1:1672246/1‑149 (MQ=255)
                       tGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAg                                                                                                                            >  1:789687/1‑149 (MQ=255)
                           gattgccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTgg                                                                                                                        <  4:1463512/149‑1 (MQ=255)
                             ttgccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGt                                                                                                                      <  2:1079811/149‑1 (MQ=255)
                                 cggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTg                                                                                                                                                >  1:1566441/1‑119 (MQ=255)
                                   gatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcg                                                                                                                <  4:1240936/149‑1 (MQ=255)
                                       ctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGt                                                                                                            <  1:149783/149‑1 (MQ=255)
                                         gattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTa                                                                                                           >  4:245121/1‑149 (MQ=255)
                                         gattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTaa                                                                                                          >  2:1494908/1‑149 (MQ=255)
                                             gccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCa                                                                                                      <  4:1345775/149‑1 (MQ=255)
                                                 gatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAg                                                                                                  >  2:703193/1‑149 (MQ=255)
                                                      gACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAg                                                                                                  >  3:1512924/1‑144 (MQ=255)
                                                      gACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGaa                                                                                             <  1:934275/149‑1 (MQ=255)
                                                         gCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                                                               >  1:669104/1‑96 (MQ=255)
                                                         gCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                                                                                                               <  2:669104/96‑1 (MQ=255)
                                                           tGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGc                                                                                        <  2:100859/149‑1 (MQ=255)
                                                             aCGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTg                                                                                      >  3:697712/1‑149 (MQ=255)
                                                             aCGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTg                                                                                      <  3:374026/149‑1 (MQ=255)
                                                              cgcgTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCAAAAACGCTGa                                                                                     <  1:1494908/149‑1 (MQ=255)
                                                              cgcgTCTCATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGa                                                                                     <  3:764100/149‑1 (MQ=255)
                                                               gcgTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcgg                                                                                                               >  3:1143947/1‑122 (MQ=255)
                                                               gcgTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcgg                                                                                                               <  4:1143947/122‑1 (MQ=255)
                                                               gcgTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGAc                                                                                    >  2:1186869/1‑149 (MQ=255)
                                                                  tCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCTCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGt                                                                                 <  4:841235/149‑1 (MQ=255)
                                                                  tCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTgg                                                                                                                        >  1:531963/1‑110 (MQ=255)
                                                                   cTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTc                                                                                <  4:201783/149‑1 (MQ=255)
                                                                      aTCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCaaa                                                                             <  2:1640261/149‑1 (MQ=255)
                                                                       tCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAg                                                                            >  2:1048111/1‑149 (MQ=255)
                                                                       tCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAg                                                                            <  1:890796/149‑1 (MQ=255)
                                                                        ccGGCCTACAGCGTCTTATCCCGCCTACATGTCCCGCCATTGTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAgc                                                                           <  4:1624011/149‑1 (MQ=255)
                                                                         cGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAgcg                                                                          >  2:1602140/1‑149 (MQ=255)
                                                                         cGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAgcg                                                                          >  1:333375/1‑149 (MQ=255)
                                                                                   cGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGcc                                                                <  2:341876/149‑1 (MQ=255)
                                                                                    ttCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCt                                                               <  2:1492200/148‑1 (MQ=255)
                                                                                    gTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCt                                                               <  1:81130/149‑1 (MQ=255)
                                                                                     tCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTg                                                              >  3:553883/1‑149 (MQ=255)
                                                                                     tCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTg                                                              >  3:1394802/1‑149 (MQ=255)
                                                                                      cTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACTGTCAAAGCGCCAGCGCGCCTGa                                                             >  3:1058689/1‑149 (MQ=255)
                                                                                       ttATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGaa                                                            >  1:377426/1‑149 (MQ=255)
                                                                                          tCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCa                                                                                                   <  4:1579762/107‑1 (MQ=255)
                                                                                          tCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCa                                                                                                   >  3:1579762/1‑107 (MQ=255)
                                                                                          tCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACATGGGCCAGACGCCTGGCCCacgt                                                                                                                                        <  4:1545444/70‑4 (MQ=37)
                                                                                            cGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCca                                                       >  4:1585560/1‑149 (MQ=255)
                                                                                                cTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAATGCGCCAGCGCGCCTGAATGCCACAAt                                                   <  4:697712/149‑1 (MQ=255)
                                                                                                 tACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATc                                                  >  1:448612/1‑149 (MQ=255)
                                                                                                  aCATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCg                                                 >  3:1249375/1‑149 (MQ=255)
                                                                                                   cATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGc                                                <  3:318541/149‑1 (MQ=255)
                                                                                                   cATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGc                                                <  3:246767/149‑1 (MQ=255)
                                                                                                         ccGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATc                                          >  1:1187726/1‑149 (MQ=255)
                                                                                                         ccGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATc                                          >  2:1310989/1‑149 (MQ=255)
                                                                                                          cGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATcc                                         >  3:1020758/1‑149 (MQ=255)
                                                                                                          cGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATcc                                         >  3:714101/1‑149 (MQ=255)
                                                                                                          cGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATcc                                         >  3:1521756/1‑149 (MQ=255)
                                                                                                          cGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATcc                                         >  2:368897/1‑149 (MQ=255)
                                                                                                          cGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATcc                                         >  1:1296544/1‑149 (MQ=255)
                                                                                                           gCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCcg                                        >  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)
                                                                                                                                                  |                                                                                                                                                    
TGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGACCGTCAAAGCGCCAGCGCGCCTGAATGCCACAATCGCTTAATCCGCGCCCTTTCGCTAAGGTCACCAGTTCACGCGATTTCTC  >  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: