Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 236 C→G 8.3% T16R (ACA→AGA) ‡ thrL → thr operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132360CG 8.3% 325.5 / 18.9 108T16R (ACA→AGA) ‡thrLthr operon leader peptide
Reads supporting (aligned to +/- strand):  ref base C (47/51);  new base G (3/6);  total (51/57)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-02

Only 100 of 110 total aligned reads displayed.

GAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGT  >  NC_000913/93‑384
                                                                                                                                               |                                                                                                                                                        
gAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTaa                                                                                                                                                    <  2:1650679/150‑1 (MQ=255)
 aGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAAc                                                                                                                                                   <  1:469441/150‑1 (MQ=255)
  gTAAAATAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACAACCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTaaag                                                                                                                                                  >  1:50389/1‑148 (MQ=255)
   tAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTCCAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACgg                                                                                                                                                 <  1:1359835/150‑1 (MQ=255)
    aaaTTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCGTTACCTGTGGTAATGGt                                                                                                                                                <  1:1458839/150‑1 (MQ=255)
    aaaTTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGt                                                                                                                                                >  2:1346276/1‑150 (MQ=255)
    aaaTTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACAACCATTACCACCACCATCACCATTACCACAGGTAACGGt                                                                                                                                                >  2:1042639/1‑150 (MQ=255)
     aaTTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTg                                                                                                                                               <  2:449058/150‑1 (MQ=255)
      aTTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGc                                                                                                                                              <  2:1523711/150‑1 (MQ=255)
       ttAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCg                                                                                                                                             <  2:456338/150‑1 (MQ=255)
       ttAAAATTTTATTGACTTAGGTAACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACAAAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCg                                                                                                                                             <  2:325961/150‑1 (MQ=255)
        tAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCgg                                                                                                                                            >  1:1098441/1‑150 (MQ=255)
        aaagtatttAGTGACCTAAGTCAATAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCgg                                                                                                                                            >  2:823668/7‑150 (MQ=255)
          aaaTTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCGTTACCTGTGGTAATGGTgatggt                                                                                                                                          <  1:472417/150‑6 (MQ=255)
          aaaTTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCGTTACCTGTGGTAATGGTgatggt                                                                                                                                          <  1:151751/150‑6 (MQ=255)
           aaTTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCGTTACCTGTGGTAATGGTgatggtg                                                                                                                                         <  2:1154501/150‑7 (MQ=255)
            aTTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCGTTACCTGTGGTAATGGTgatggtgg                                                                                                                                        <  1:1438047/150‑8 (MQ=255)
               ttATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGAcg                                                                                                                                     <  2:352459/150‑1 (MQ=255)
                    gACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTAc                                                                                                                                >  1:435786/1‑150 (MQ=255)
                     aCTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGAAACGGTGCGGGCTGACGCGTACa                                                                                                                               <  2:386667/150‑1 (MQ=255)
                           gTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAac                                                                                                                         <  2:825606/150‑1 (MQ=255)
                           gTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAac                                                                                                                         >  1:1662501/1‑150 (MQ=255)
                           gTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAac                                                                                                                         <  1:1126133/150‑1 (MQ=255)
                           gTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAac                                                                                                                         >  2:178938/1‑150 (MQ=255)
                              aCTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAAcaca                                                                                                                      >  2:794818/1‑150 (MQ=255)
                                tAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGa                                                                                                                   <  1:1013869/150‑1 (MQ=255)
                                 aaaTACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGaa                                                                                                                  >  2:1359183/1‑150 (MQ=255)
                                     aCTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGaaaaaa                                                                                                              >  2:1467998/1‑150 (MQ=255)
                                     aCTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGaaaaaa                                                                                                              >  1:1178123/1‑150 (MQ=255)
                                      cTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAg                                                                                                             <  2:1178123/150‑1 (MQ=255)
                                       tttAACCAATATAGGCTTAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGc                                                                                                            <  2:160668/150‑1 (MQ=255)
                                        ttAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGcc                                                                                                           <  1:823668/150‑1 (MQ=255)
                                           aCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGc                                                                                                        >  2:1248978/1‑150 (MQ=255)
                                           aCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGc                                                                                                        <  2:1063421/150‑1 (MQ=255)
                                           aCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGc                                                                                                        >  1:146419/1‑150 (MQ=255)
                                             cAATATAGGCATAGCGCCCAGACAGATAAAAAATACCGAGTACACAACATCCATGAAACGAAATACCACCAACATTACCACAACCATCACCATTTCCACAGGGAACGGTGCGGCATGACGCGAACAGGAAGCACAGAACAAAACCCCCAc                                                                                                      >  1:1142372/1‑150 (MQ=255)
                                                  taGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGAc                                                                                                 >  1:993145/1‑150 (MQ=255)
                                                     gCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGt                                                                                             >  1:180572/1‑150 (MQ=255)
                                                        tAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCg                                                                                          >  2:864179/1‑150 (MQ=255)
                                                         aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCgg                                                                                         <  1:335805/150‑1 (MQ=255)
                                                          gcgcACAGACAGATAAAAATTACAGAGTACACAACATCGATGAAACGCATAAGCACCACCATTACCACCACCATCACCATTACCACAGGTATCGCGGGGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCggg                                                                                        <  2:1406653/150‑1 (MQ=255)
                                                          gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCggg                                                                                        <  2:1638288/150‑1 (MQ=255)
                                                           cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACTGTCAGGTGCGGGc                                                                                       <  2:1098441/150‑1 (MQ=255)
                                                           cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACTGTCAGGTGCGGGc                                                                                       <  1:1467998/150‑1 (MQ=255)
                                                            gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCt                                                                                     >  1:1422712/1‑150 (MQ=255)
                                                            gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACTGTCAGGTGCGGGCt                                                                                     <  1:1617580/150‑1 (MQ=255)
                                                             cacaGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCtt                                                                                    <  2:1037534/150‑1 (MQ=255)
                                                              acagacagATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACTGTCAGGTGCGGGCttt                                                                                   <  1:705800/150‑1 (MQ=255)
                                                                agacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttttt                                                                                <  1:1615715/150‑1 (MQ=255)
                                                                 gacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCCGACAGTGCGGGCtttttt                                                                               >  1:1468242/1‑150 (MQ=255)
                                                                  acagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGTTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttttttt                                                                              <  1:864179/150‑1 (MQ=255)
                                                                   cagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCGGACAGTGCGGGCtttttttt                                                                             >  1:186766/1‑150 (MQ=255)
                                                                   attaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACTGTCAGGTGCGGGCttttttct                                                                             <  2:1662501/147‑3 (MQ=39)
                                                                    agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttttttttt                                                                            <  2:1183107/150‑1 (MQ=255)
                                                                    agaTAAAAATTACAGAGGACTCAACATCCATGAACCGCATTAGCACCACCATAACCACCACCATCACCATAACCACAGGTAACGGTGCGGGCTGAAGCGTACAGGAAACACAGAAAAAAAGCCCGCACTGTCAGGTGCGGGCttttttct                                                                             <  2:278534/150‑3 (MQ=37)
                                                                     gaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTc                                                                           >  2:875165/1‑150 (MQ=255)
                                                                      aTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACTGTCAGGTGCGGGCttttttctgtg                                                                          <  2:149740/150‑6 (MQ=39)
                                                                       tAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCgc                                                                         >  1:941855/1‑149 (MQ=255)
                                                                        aaaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGATACACCGAAAAAAGCCCGCACTGTCAGGTGCGGGCttttttctgtgtt                                                                        <  1:271591/150‑8 (MQ=38)
                                                                        aaaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGAc                                                                        <  1:592041/150‑1 (MQ=255)
                                                                         aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATAACCACAGGTAACGGTGCGGGCTGACGGGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGAGCGGGCTTTTTTTTTCGagc                                                                       >  2:837152/1‑148 (MQ=255)
                                                                                cAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTa                                                                <  1:1346276/150‑1 (MQ=255)
                                                                                cAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTa                                                                <  1:1462508/150‑1 (MQ=255)
                                                                                     tACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCAACACCATTACCACAGGTAACGGTGCGGGCTGACGCGAACAGGAAACACAGAAAAAAGCCCGCACATGACAGTGCGGGCTTTTTTTTTCGACCAAAGATAACCAg                                                           >  1:88768/1‑150 (MQ=255)
                                                                                       cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGt                                                         >  1:1481330/1‑150 (MQ=255)
                                                                                            caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                    >  1:790084/1‑150 (MQ=255)
                                                                                            caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGGAACGAGGTaacaa                                                    >  2:852408/1‑150 (MQ=255)
                                                                                              aacATTAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACATAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAAcc                                                  <  2:1439670/148‑1 (MQ=255)
                                                                                               ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                 <  2:701122/150‑1 (MQ=255)
                                                                                                  tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCCACCAAAGGTACCGAGGTAACAACCATGc                                              >  1:313501/1‑150 (MQ=255)
                                                                                                     aaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                           <  1:178938/150‑1 (MQ=255)
                                                                                                       cGCATTAGCACCACCATTACCACCACCATCACCATTACCTCAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACTGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgtg                                         <  2:1408212/150‑1 (MQ=255)
                                                                                                            tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                    <  1:1665723/150‑1 (MQ=255)
                                                                                                             aGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAg                                   >  1:789747/1‑150 (MQ=255)
                                                                                                              gCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGt                                  >  2:1501786/1‑150 (MQ=255)
                                                                                                                accaccATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGCCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                >  2:1525902/1‑150 (MQ=255)
                                                                                                                 caaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGATACACTGAAAAAAGGCGGCACCTGACAGTGCGGGCTTTTTTTTTCGACCATAGGTAACGAGGTAACAAGCATGCGAGTGATGAAGTTcg                               <  2:939580/148‑1 (MQ=255)
                                                                                                                  caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                              <  1:1680191/150‑1 (MQ=255)
                                                                                                                    ccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACCACCATGCGAGTGGTGAAGTTcggcg                            >  2:1205127/1‑150 (MQ=255)
                                                                                                                       ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                         <  1:1359/150‑1 (MQ=255)
                                                                                                                       ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                         <  1:1359183/150‑1 (MQ=255)
                                                                                                                       ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                         <  2:568944/150‑1 (MQ=255)
                                                                                                                        tACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAc                        <  1:794818/150‑1 (MQ=255)
                                                                                                                        tACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAc                        <  2:1422712/150‑1 (MQ=255)
                                                                                                                          ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAt                      >  1:870123/1‑150 (MQ=255)
                                                                                                                          caaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCACGCACCTGACAGTGGGGGCTGTTTTTTTCGACCAAAGGTAAAGAGGTAACACCACTGCGAGTTTTGAAGTTAGGCGGTACAt                      >  2:1280236/3‑150 (MQ=255)
                                                                                                                           caccaccaTAACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACAAAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                     <  2:88768/150‑1 (MQ=255)
                                                                                                                             ccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAg                   >  1:931054/1‑150 (MQ=255)
                                                                                                                                 ccgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTATATTTCGACCAAAGGTAACGATGTAACAACCATGCGAGTGTTGAAGTTAGGCGGTACAACAGTGGc               >  2:589140/4‑150 (MQ=255)
                                                                                                                                 catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc               >  1:1573245/1‑150 (MQ=255)
                                                                                                                                  cgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa              >  2:1243236/3‑150 (MQ=255)
                                                                                                                                    caccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTGTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaa             >  1:1109522/1‑150 (MQ=255)
                                                                                                                                       caatACAACAGGTAACGGTCCGGGCTGACGGGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCATTGGCAAATg          <  2:878638/147‑1 (MQ=255)
                                                                                                                                        atTACCACAGGTAACGGTGCGGGCTGACGGGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa        <  1:1609719/150‑1 (MQ=255)
                                                                                                                                           aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa     >  1:762714/1‑150 (MQ=255)
                                                                                                                                           aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAAAAAACATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa     >  1:562923/1‑150 (MQ=255)
                                                                                                                                           aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGAAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATTAAGaa     >  1:762716/1‑150 (MQ=255)
                                                                                                                                           aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTGTTTTTTTCGACCAAAGGTAACGAGTGTACAACCATGCGAGTGTTGAAGTTGGGGGGTACATAAGTGGCAAATGCAGaa     >  2:1089475/1‑150 (MQ=255)
                                                                                                                                            ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc    >  1:1507716/1‑150 (MQ=255)
                                                                                                                                              acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATACAGAACGt  <  2:1483218/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                        
GAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGT  >  NC_000913/93‑384

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