Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,121,962 C→A 9.9% G109V (GGT→GTT)  hslV ← peptidase component of the HslVU protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,121,9620CA 9.9% 285.6 / 14.4 91G109V (GGT→GTT) hslVpeptidase component of the HslVU protease
Reads supporting (aligned to +/- strand):  ref base C (39/43);  new base A (7/2);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

TCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGTTCAAACA  >  NC_000913/4121815‑4122109
                                                                                                                                                   |                                                                                                                                                   
tCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCg                                                                                                                                                   >  1:1632879/1‑150 (MQ=255)
 cGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCgg                                                                                                                                                  <  1:136701/150‑1 (MQ=255)
  gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGt                                                                                                                                                 >  1:307439/1‑150 (MQ=255)
       cAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatga                                                                                                                                            >  1:1150081/1‑150 (MQ=255)
       cAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatga                                                                                                                                            >  2:242001/1‑150 (MQ=255)
         atatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGGTAACGGtgacgtgg                                                                                                                                          <  1:302272/150‑6 (MQ=255)
          tatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATaa                                                                                                                                         >  1:1282926/1‑150 (MQ=255)
           atCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGGTAACGGtgacgtggtg                                                                                                                                        <  2:359038/150‑8 (MQ=255)
            tCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGc                                                                                                                                       <  1:1397751/150‑1 (MQ=255)
             ccAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCg                                                                                                                                      >  1:765914/1‑150 (MQ=255)
                      tttCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGttt                                                                                                                             <  1:259793/150‑1 (MQ=255)
                       ttCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGTCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTc                                                                                                                            <  1:1224596/150‑1 (MQ=255)
                        tCAGCAATTTCACGGCCGCTAAGTGCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAACATCGTGTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGCTTCa                                                                                                                           <  2:1162630/150‑1 (MQ=255)
                          aGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATc                                                                                                                         >  2:707500/1‑150 (MQ=255)
                          aGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATc                                                                                                                         >  1:323671/1‑150 (MQ=255)
                              aTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCcgcg                                                                                                                     <  1:1430223/150‑1 (MQ=255)
                              aTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCcgcg                                                                                                                     <  1:138756/150‑1 (MQ=255)
                                ttCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCcgcggc                                                                                                                   >  1:1258426/1‑148 (MQ=255)
                                ttCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGAc                                                                                                                   <  1:1544568/150‑1 (MQ=255)
                                ttCACGGGCGCTAAGTTCAGTGTTTTCTAAAA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTg                                                                                                                 >  2:1015966/1‑150 (MQ=255)
                                  cACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTg                                                                                                                 <  2:307439/150‑1 (MQ=255)
                                  cACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTg                                                                                                                 <  2:1150081/150‑1 (MQ=255)
                                       gcgcTAATTTCAGTGTTTTCTAACAGGGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCcagc                                                                                                            <  2:443306/150‑1 (MQ=255)
                                          cTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGt                                                                                                         >  1:1321844/1‑150 (MQ=255)
                                              gTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCtt                                                                                                     <  1:1519512/150‑1 (MQ=255)
                                              gTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCtt                                                                                                     >  2:1266465/1‑150 (MQ=255)
                                              gTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCtt                                                                                                     <  1:242001/150‑1 (MQ=255)
                                                tCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCa                                                                                                   <  1:1474808/150‑1 (MQ=255)
                                                    tgtTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGtt                                                                                               <  1:930688/150‑1 (MQ=255)
                                                      ttttCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTg                                                                                             <  2:18509/150‑1 (MQ=255)
                                                       tttCTATCAGCGCGCGGGGCGCATCCTGGTTGTATGGGCCGGCGGAGCCGATAGGTATAAGATCGTTTTCTGGGTGCACTACGTCTCCGTTCCCGGTGATGATATGGGTTGCAGTTTCTTCGGCGAGTGCCAGCAGTGCTTCAAGTTTgc                                                                                            <  2:311755/150‑1 (MQ=255)
                                                       tttCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTgc                                                                                            >  2:990287/1‑150 (MQ=255)
                                                       tttCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTgc                                                                                            >  2:1267391/1‑150 (MQ=255)
                                                       tttCTAACAGCGCCCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCAGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTATCATCCGCGACTGCCAGCAGTGCTTCAAGTTTgc                                                                                            >  2:1017672/1‑150 (MQ=255)
                                                         tCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTgcgc                                                                                          <  1:1625288/150‑1 (MQ=255)
                                                           tAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCgcag                                                                                        >  1:1217270/1‑150 (MQ=255)
                                                           tAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCgcag                                                                                        >  1:426757/1‑150 (MQ=255)
                                                            aaCAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCgcagc                                                                                       <  1:928516/150‑1 (MQ=255)
                                                              cAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCAt                                                                                     >  2:808459/1‑150 (MQ=255)
                                                               aGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATg                                                                                    >  1:415628/1‑150 (MQ=255)
                                                                gcgcgcgcgcCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGc                                                                                   <  2:415628/150‑1 (MQ=255)
                                                                 cgcgcgcgcCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCg                                                                                  >  2:568350/1‑150 (MQ=255)
                                                                       cgccgcAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGt                                                                            >  1:717854/1‑150 (MQ=255)
                                                                       cgccgcAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGGCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGt                                                                            <  2:1511674/150‑1 (MQ=255)
                                                                        gccgcAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTa                                                                           <  2:667598/150‑1 (MQ=255)
                                                                           gcAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGc                                                                        <  1:5388/150‑1 (MQ=255)
                                                                              gCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAg                                                                     <  1:1015966/150‑1 (MQ=255)
                                                                              gCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAg                                                                     >  2:1083451/1‑150 (MQ=255)
                                                                                   ggCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACTGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCttt                                                                <  1:638835/150‑1 (MQ=255)
                                                                                     cGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGCCTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTg                                                              >  2:519779/1‑150 (MQ=255)
                                                                                          gggCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGc                                                         <  1:1254799/150‑1 (MQ=255)
                                                                                           ggCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCt                                                        <  2:765914/150‑1 (MQ=255)
                                                                                           ggCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCt                                                        >  2:1169077/1‑150 (MQ=255)
                                                                                            gccgccGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTc                                                       >  2:847919/1‑150 (MQ=255)
                                                                                            gccgccGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTc                                                       <  2:1321844/150‑1 (MQ=255)
                                                                                                  ggAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAAcggc                                                 >  2:717985/1‑150 (MQ=255)
                                                                                                    aGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAAcggcgg                                               >  1:1423735/1‑150 (MQ=255)
                                                                                                    aGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAAcggcgg                                               <  2:583270/150‑1 (MQ=255)
                                                                                                    aGCCGATAGCAATAAGAACGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAAcggcgg                                               <  2:600751/150‑1 (MQ=255)
                                                                                                      ccGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCt                                             >  2:248163/1‑150 (MQ=255)
                                                                                                      ccGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCt                                             >  1:1300523/1‑150 (MQ=255)
                                                                                                         aTAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTg                                          <  2:560049/150‑1 (MQ=255)
                                                                                                            gCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGAcc                                       >  2:644403/1‑150 (MQ=255)
                                                                                                             cAATAAGATCTTTTTCTGGCTCCACCACGTCACCCTTACAGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCa                                      >  2:88922/1‑150 (MQ=255)
                                                                                                              aaTAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAg                                     >  1:394966/1‑150 (MQ=255)
                                                                                                                     tCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGccc                              <  2:1278709/150‑1 (MQ=255)
                                                                                                                     tCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGccc                              <  2:717854/150‑1 (MQ=255)
                                                                                                                        ttttCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGa                           <  2:596538/150‑1 (MQ=255)
                                                                                                                         tttCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGGCTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGAt                          <  2:1258426/150‑1 (MQ=255)
                                                                                                                         tttCTGGCTGCACCACGTCACCGTTACCCGTGATGATAAGCGATGCAGTTGCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGAt                          <  2:781675/150‑1 (MQ=255)
                                                                                                                            cTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCa                       <  2:323671/150‑1 (MQ=255)
                                                                                                                              ggCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCAtt                     >  1:660954/1‑150 (MQ=255)
                                                                                                                              ggCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCAtt                     >  2:632489/1‑150 (MQ=255)
                                                                                                                               gCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCAttt                    >  1:745820/1‑150 (MQ=255)
                                                                                                                                cTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTc                   <  1:1169077/150‑1 (MQ=255)
                                                                                                                                  gCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCa                 >  2:520639/1‑150 (MQ=255)
                                                                                                                                  gCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCa                 >  2:923154/1‑150 (MQ=255)
                                                                                                                                   caccacGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAg                <  1:1083451/150‑1 (MQ=255)
                                                                                                                                    ttatcatCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGt               >  2:910469/7‑150 (MQ=255)
                                                                                                                                     tatcatCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGtt              >  1:1549482/6‑150 (MQ=255)
                                                                                                                                     tatcatCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGtt              >  1:187324/6‑150 (MQ=255)
                                                                                                                                      cacGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGttt             >  1:1216765/1‑150 (MQ=255)
                                                                                                                                      atcatCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGttt             >  2:474077/5‑150 (MQ=255)
                                                                                                                                      atcatCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGttt             >  1:141447/5‑150 (MQ=255)
                                                                                                                                      atcatCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGttt             >  1:583688/5‑150 (MQ=255)
                                                                                                                                       tcatCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTa            >  2:532192/4‑150 (MQ=255)
                                                                                                                                       acGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTa            <  1:632489/150‑1 (MQ=255)
                                                                                                                                           ctgcGTTACCGTTGATGATAATCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACTCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATTGCCCTGATGCATTTCCAGTTTACGtt        <  1:605680/147‑1 (MQ=255)
                                                                                                                                            aCCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGTTc       >  2:596916/1‑150 (MQ=255)
                                                                                                                                            aCCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGAATTCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGTTc       <  2:1323600/150‑1 (MQ=255)
                                                                                                                                                 tACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGTTCAAACa  <  2:660954/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGTTCAAACA  >  NC_000913/4121815‑4122109

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