Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,720,240 A→T 6.4% N122I (AAT→ATT)  slyB → outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,720,2400AT 6.4% 552.4 / 21.1 157N122I (AAT→ATT) slyBouter membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (64/83);  new base T (5/5);  total (69/88)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

Only 100 of 162 total aligned reads displayed.

CTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAT  >  NC_000913/1720095‑1720388
                                                                                                                                                 |                                                                                                                                                    
 tGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGCCGATGGTAATAc                                                                                                                                                   <  2:490908/148‑1 (MQ=255)
 tGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATAcc                                                                                                                                                  <  4:1444735/149‑1 (MQ=255)
    ggAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACcatc                                                                                                                                               >  1:1632785/1‑149 (MQ=255)
       aTACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTATTACcatcgtc                                                                                                                                            <  4:761025/149‑4 (MQ=255)
       aTACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTATTACcatcgtc                                                                                                                                            <  3:1161383/149‑4 (MQ=255)
        tACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACtat                                                                                                                                                >  3:1450456/1‑144 (MQ=255)
        tACTGTTGGTGGCGGAACCGGGCCTTCTCTGGCTACTGCAGCAGGCGCTGTTGTAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACcat                                                                                                                                                <  4:1450456/144‑1 (MQ=255)
         aCTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCAtggt                                                                                                                                          <  1:1109980/149‑1 (MQ=255)
          cTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCAtggtg                                                                                                                                         <  4:1539338/149‑1 (MQ=255)
               ggtggCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACa                                                                                                                                    >  4:1652436/1‑149 (MQ=255)
                     ggAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACa                                                                                                                              <  1:1275752/149‑1 (MQ=255)
                      gAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACaa                                                                                                                             <  3:699479/149‑1 (MQ=255)
                      gAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACaa                                                                                                                             <  2:1632785/149‑1 (MQ=255)
                      gAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACaa                                                                                                                             <  1:132544/149‑1 (MQ=255)
                       aaCCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGGTTGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAg                                                                                                                            <  3:183663/149‑1 (MQ=255)
                       aaCCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAg                                                                                                                            <  2:1426891/149‑1 (MQ=255)
                       aaCCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAg                                                                                                                            <  2:1098131/149‑1 (MQ=255)
                       aaCCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAg                                                                                                                            <  1:1218637/149‑1 (MQ=255)
                        aCCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAgg                                                                                                                           >  2:1559289/1‑149 (MQ=255)
                        aCCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAgg                                                                                                                           <  4:1005342/149‑1 (MQ=255)
                             gCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAaca                                                                                                                      >  2:540007/1‑149 (MQ=255)
                               gTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACaa                                                                                                                             >  4:1371389/1‑140 (MQ=255)
                                ttCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTc                                                                                                                   <  4:919714/149‑1 (MQ=255)
                                  ctctGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGt                                                                                                                 <  2:1030775/149‑1 (MQ=255)
                                     tGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTtc                                                                                                              >  3:1342487/1‑149 (MQ=255)
                                       gCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTtctc                                                                                                            <  2:883720/149‑1 (MQ=255)
                                         tACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTtctctc                                                                                                          >  4:1164878/1‑149 (MQ=255)
                                            tGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCgg                                                                                                       <  2:461994/149‑1 (MQ=255)
                                             gcagcaGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCggg                                                                                                      >  4:1440135/1‑149 (MQ=255)
                                                   ggCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACg                                                                                                >  2:1373263/1‑149 (MQ=255)
                                                   ggCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACg                                                                                                >  2:1077297/1‑149 (MQ=255)
                                                    gcgcTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCaa                                                                                                                        <  1:1076840/124‑1 (MQ=255)
                                                      gcTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCAtggtgg                                                                                                                                        >  3:861148/1‑106 (MQ=255)
                                                       cTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTc                                                                                            <  2:402570/149‑1 (MQ=255)
                                                          ttGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTa                                                                                         <  3:1652436/149‑1 (MQ=255)
                                                                  ggCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAg                                                                                 <  1:426859/149‑1 (MQ=255)
                                                                      tAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTTCAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAAt                                                                             <  3:1139639/149‑1 (MQ=255)
                                                                      tAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAAt                                                                             <  1:1514753/149‑1 (MQ=255)
                                                                       aGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATg                                                                            <  1:467785/149‑1 (MQ=255)
                                                                          tGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCa                                                                         >  2:696337/1‑149 (MQ=255)
                                                                            gTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGt                                                                       >  4:1344547/1‑149 (MQ=255)
                                                                               aGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAg                                                                    <  1:1330391/149‑1 (MQ=255)
                                                                                 ggCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGt                                                                  <  2:1084905/149‑1 (MQ=255)
                                                                                 ggCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGt                                                                  >  2:277996/1‑149 (MQ=255)
                                                                                       cAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGt                                                                       <  2:74409/138‑1 (MQ=255)
                                                                                       cAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGt                                                                       >  1:74409/1‑138 (MQ=255)
                                                                                           gTGCAATGAACAAAACGCAGGGTATCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTtct                                                        <  1:1641333/149‑1 (MQ=255)
                                                                                              cAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCcg                                                     >  2:1012507/1‑149 (MQ=255)
                                                                                                aTGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAg                                                                                                                            >  3:1188026/1‑76 (MQ=255)
                                                                                                 tGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCgcgc                                                  >  4:564876/1‑149 (MQ=255)
                                                                                                   aaCAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGc                                                                                    >  4:48989/1‑113 (MQ=255)
                                                                                                   aaCAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTa                                                <  1:696337/149‑1 (MQ=255)
                                                                                                   aaCAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTa                                                >  3:757897/1‑149 (MQ=255)
                                                                                                     cAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTc                                                                                            >  2:1234650/1‑103 (MQ=255)
                                                                                                     cAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTc                                                                                            <  1:1234650/103‑1 (MQ=255)
                                                                                                     cAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTaaa                                              >  4:1031702/1‑149 (MQ=255)
                                                                                                      aaaaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTtctc                                                                                                            <  3:105449/86‑1 (MQ=255)
                                                                                                      aaaaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTtctc                                                                                                            >  4:105449/1‑86 (MQ=255)
                                                                                                      aaaaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAAt                                             >  1:495742/1‑149 (MQ=255)
                                                                                                      aaaaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAAt                                             >  1:91216/1‑149 (MQ=255)
                                                                                                      aaaaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAAt                                             >  1:863091/1‑149 (MQ=255)
                                                                                                       aaaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATa                                            >  4:621094/1‑149 (MQ=255)
                                                                                                       aaaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATa                                            >  1:959683/1‑149 (MQ=255)
                                                                                                        aaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATaa                                           <  1:1077297/149‑1 (MQ=255)
                                                                                                        aaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATaa                                           <  3:1164878/149‑1 (MQ=255)
                                                                                                        aaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATaa                                           <  3:1028618/149‑1 (MQ=255)
                                                                                                        aaCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATaa                                           <  1:1447389/149‑1 (MQ=255)
                                                                                                         aCGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAg                                                                    >  3:516448/1‑123 (MQ=255)
                                                                                                          cGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGt                                         <  3:503082/149‑1 (MQ=255)
                                                                                                            cAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGttt                                       <  1:443830/149‑1 (MQ=255)
                                                                                                              gggTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGc                                     <  3:384297/149‑1 (MQ=255)
                                                                                                              gggTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGc                                     >  4:1212079/1‑149 (MQ=255)
                                                                                                                gtgtCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAAt                                             >  3:1627296/1‑139 (MQ=255)
                                                                                                                    cGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGt                                                            >  2:1202808/1‑120 (MQ=255)
                                                                                                                     gAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGt                              <  3:1323771/149‑1 (MQ=255)
                                                                                                                      aGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGGAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTc                             >  1:1499094/1‑149 (MQ=255)
                                                                                                                       gCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCa                            <  1:636804/149‑1 (MQ=255)
                                                                                                                        cTGGAAATTCGTAAAGACGATGGTAATACTATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTACGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGAGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAg                           >  4:184665/1‑149 (MQ=255)
                                                                                                                        cTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAg                           <  2:1499094/149‑1 (MQ=255)
                                                                                                                         tGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAgg                          >  1:277677/1‑149 (MQ=255)
                                                                                                                              aTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAgtaacct                     <  4:846345/149‑7 (MQ=255)
                                                                                                                                 cGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGAc                  <  2:122558/149‑1 (MQ=255)
                                                                                                                                   tAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACgtgt                                                                                             <  3:823077/72‑1 (MQ=255)
                                                                                                                                   tAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACgtgt                                                                                             >  4:823077/1‑72 (MQ=255)
                                                                                                                                      ccatgATGGTATTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCcgc                                                    <  2:1683373/106‑1 (MQ=255)
                                                                                                                                      ccatgATGGTATTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCcgc                                                    >  1:1683373/5‑110 (MQ=255)
                                                                                                                                      aGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGTAACCTGACCacac             <  1:1148016/149‑1 (MQ=255)
                                                                                                                                      aGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCacac             <  3:414133/149‑1 (MQ=255)
                                                                                                                                       catgATGGTATTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGt                                                                  >  4:333932/4‑95 (MQ=255)
                                                                                                                                       catgATGGTATTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGt                                                                  <  3:333932/92‑1 (MQ=255)
                                                                                                                                          gATGGTATTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAAGt                                                                  >  1:477350/1‑92 (MQ=255)
                                                                                                                                          gATGGTATTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAAGt                                                                  <  2:477350/92‑1 (MQ=255)
                                                                                                                                          gATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCGGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGtt                        <  2:749075/134‑1 (MQ=255)
                                                                                                                                          gATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCGGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGtt                        >  1:749075/1‑134 (MQ=255)
                                                                                                                                          gATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGccc         >  2:662552/1‑149 (MQ=255)
                                                                                                                                               tAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTAc                      <  3:1472649/131‑1 (MQ=255)
                                                                                                                                               tAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAAGGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTAc                      >  4:1472649/1‑131 (MQ=255)
                                                                                                                                                aaTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGt                         <  2:86998/127‑1 (MQ=255)
                                                                                                                                                aaTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCa   >  4:460077/1‑149 (MQ=255)
                                                                                                                                                 aTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAt  <  3:1435486/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACTGCAGCAGGCGCTGTTGCAGGTGGCGTAGCTGGTCAGGGCGTACAGAGTGCAATGAACAAAACGCAGGGTGTCGAGCTGGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAT  >  NC_000913/1720095‑1720388

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