Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 235 A→T 9.0% T16S (ACA→TCA) ‡ thrL → thr operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132350AT 9.0% 260.1 / 14.0 108T16S (ACA→TCA) ‡thrLthr operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (48/50);  new base T (3/7);  total (51/57)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-02

Only 100 of 111 total aligned reads displayed.

CTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGT  >  NC_000913/86‑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:1462508/1‑150 (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)
                 aaaTTTTATTGACTTAGGTCACTAAATACTTTACCCAATATAGGCATAGCGCACAGACAGATAAAAATAACAGAGTACACAACATCAAAGAAACGCATTAGCAACACCAATACCACCACCATAACAATTACCAAAGGTAACGGTGAGGGc                                                                                                                                         >  1:93108/1‑150 (MQ=255)
                 aaaTTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCGTTACCTGTGGTAATGGTgatggt                                                                                                                                         <  1:472417/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:1126133/150‑1 (MQ=255)
                                  gTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAac                                                                                                                        >  2:178938/1‑150 (MQ=255)
                                  gTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAac                                                                                                                        >  1:1662501/1‑150 (MQ=255)
                                       tAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGa                                                                                                                   <  1:1013869/150‑1 (MQ=255)
                                        aaaTACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGaa                                                                                                                  >  2:1359183/1‑150 (MQ=255)
                                            aCTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGaaaaaa                                                                                                              >  1:1178123/1‑150 (MQ=255)
                                            aCTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGaaaaaa                                                                                                              >  2:1467998/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                                                                                                        >  1:146419/1‑150 (MQ=255)
                                                  aCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGc                                                                                                        >  2:1248978/1‑150 (MQ=255)
                                                    cAATATAGGCATAGCGCCCAGACAGATAAAAAATACCGAGTACACAACATCCATGAAACGAAATACCACCAACATTACCACAACCATCACCATTTCCACAGGGAACGGTGCGGCATGACGCGAACAGGAAGCACAGAACAAAACCCCCAc                                                                                                      >  1:1142372/1‑150 (MQ=255)
                                                         taGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGAc                                                                                                 >  1:993145/1‑150 (MQ=255)
                                                               tAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCg                                                                                          >  2:864179/1‑150 (MQ=255)
                                                                aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                         <  1:335805/150‑1 (MQ=255)
                                                                 gcgcACAGACAGATAAAAATTACAGAGTACACAACATCGATGAAACGCATAAGCACCACCATTACCACCACCATCACCATTACCACAGGTATCGCGGGGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCggg                                                                                        <  2:1406653/150‑1 (MQ=255)
                                                                 gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCggg                                                                                        <  2:1638288/150‑1 (MQ=255)
                                                                  cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGc                                                                                       <  2:1098441/150‑1 (MQ=255)
                                                                  cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGc                                                                                       <  1:1467998/150‑1 (MQ=255)
                                                                   gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑t                                                                                     >  1:1422712/1‑150 (MQ=255)
                                                                   gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑t                                                                                     <  1:1617580/150‑1 (MQ=255)
                                                                    cacaGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑tt                                                                                    <  2:1037534/150‑1 (MQ=255)
                                                                     acagacagATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑ttt                                                                                   <  1:705800/150‑1 (MQ=255)
                                                                       agacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttt‑tt                                                                                <  1:1615715/150‑1 (MQ=255)
                                                                        gacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCCGACA‑GTGCGGGC‑ttt‑ttt                                                                               >  1:1468242/1‑150 (MQ=255)
                                                                         acagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGTTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttt‑tttt                                                                              <  1:864179/150‑1 (MQ=255)
                                                                          cagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttt‑ttttt                                                                             >  2:1665723/1‑150 (MQ=255)
                                                                          cagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCGGACA‑GTGCGGGC‑ttt‑ttttt                                                                             >  1:186766/1‑150 (MQ=255)
                                                                          attaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑ttt‑tttct                                                                             <  2:1662501/147‑3 (MQ=39)
                                                                           agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttt‑tttttt                                                                            <  2:1183107/150‑1 (MQ=255)
                                                                            gaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTc                                                                           >  2:875165/1‑150 (MQ=255)
                                                                              tAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCgc                                                                         >  1:941855/1‑149 (MQ=255)
                                                                               aaaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGATACACCGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑ttt‑tttctgtgtt                                                                        <  1:271591/150‑8 (MQ=38)
                                                                               aaaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGAc                                                                        <  1:592041/150‑1 (MQ=255)
                                                                                aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATAACCACAGGTAACGGTGCGGGCTGACGGGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GAGCGGGC‑TTT‑TTTTTTCGagc                                                                       >  2:837152/1‑148 (MQ=255)
                                                                                       cAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTa                                                                <  1:1346276/150‑1 (MQ=255)
                                                                                       cAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTa                                                                <  1:1462508/150‑1 (MQ=255)
                                                                                            tACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCAACACCATTACCACAGGTAACGGTGCGGGCTGACGCGAACAGGAAACACAGAAAAAAGCCCGCACATGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGATAACCAg                                                           >  1:88768/1‑150 (MQ=255)
                                                                                              cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGt                                                         >  1:1481330/1‑150 (MQ=255)
                                                                                                   caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTaacaa                                                    >  1:790084/1‑150 (MQ=255)
                                                                                                   caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGGAACGAGGTaacaa                                                    >  2:852408/1‑150 (MQ=255)
                                                                                                     aacATTAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACATAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAAcc                                                  <  2:1439670/148‑1 (MQ=255)
                                                                                                      ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                 <  2:707146/150‑1 (MQ=255)
                                                                                                      ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                 <  2:701122/150‑1 (MQ=255)
                                                                                                         tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCCACCAAAGGTACCGAGGTAACAACCATGc                                              >  1:313501/1‑150 (MQ=255)
                                                                                                            aaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                           <  1:178938/150‑1 (MQ=255)
                                                                                                              cGCATTAGCACCACCATTACCACCACCATCACCATTACCTCAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACTGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgtg                                         <  2:1408212/150‑1 (MQ=255)
                                                                                                                   tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                    <  1:1665723/150‑1 (MQ=255)
                                                                                                                    aGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAg                                   >  1:789747/1‑150 (MQ=255)
                                                                                                                     gCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGt                                  >  2:1501786/1‑150 (MQ=255)
                                                                                                                       accaccATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGCCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                >  2:1525902/1‑150 (MQ=255)
                                                                                                                        ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                               >  2:1192090/1‑150 (MQ=255)
                                                                                                                        caaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGATACACTGAAAAAAGGCGGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCATAGGTAACGAGGTAACAAGCATGCGAGTGATGAAGTTcg                               <  2:939580/148‑1 (MQ=255)
                                                                                                                         caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                              <  1:1680191/150‑1 (MQ=255)
                                                                                                                           ccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACCACCATGCGAGTGGTGAAGTTcggcg                            >  2:1205127/1‑150 (MQ=255)
                                                                                                                              ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                         <  1:1359/150‑1 (MQ=255)
                                                                                                                              ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                         <  2:568944/150‑1 (MQ=255)
                                                                                                                              ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                         <  1:1359183/150‑1 (MQ=255)
                                                                                                                               tACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAc                        <  2:1422712/150‑1 (MQ=255)
                                                                                                                               tACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAc                        <  1:794818/150‑1 (MQ=255)
                                                                                                                                 ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAt                      >  1:870123/1‑150 (MQ=255)
                                                                                                                                  caccaccaTAACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACAAAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                     <  2:88768/150‑1 (MQ=255)
                                                                                                                                    ccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAg                   >  1:931054/1‑150 (MQ=255)
                                                                                                                                        ccgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑ATATTTCGACCAAAGGTAACGATGTAACAACCATGCGAGTGTTGAAGTTAGGCGGTACAACAGTGGc               >  2:589140/4‑150 (MQ=255)
                                                                                                                                        catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc               >  1:1573245/1‑150 (MQ=255)
                                                                                                                                         cgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa              >  2:1243236/3‑150 (MQ=255)
                                                                                                                                           caccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTGTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaa             >  1:1109522/1‑150 (MQ=255)
                                                                                                                                            accatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑GTT‑TTTTTTCGCCCAAAGGTAACGAGGTTACAACCATGCGAGTGTTTAAGTACGGCGGTACATCAGTGGCAAAt           >  1:1528644/1‑150 (MQ=255)
                                                                                                                                              caatACAACAGGTAACGGTCCGGGCTGACGGGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCATTGGCAAATg          <  2:878638/147‑1 (MQ=255)
                                                                                                                                               atTACCACAGGTAACGGTGCGGGCTGACGGGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa        <  1:1609719/150‑1 (MQ=255)
                                                                                                                                                  aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa     >  1:762714/1‑150 (MQ=255)
                                                                                                                                                  aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGAAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATTAAGaa     >  1:762716/1‑150 (MQ=255)
                                                                                                                                                  aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TGT‑TTTTTTCGACCAAAGGTAACGAGTGTACAACCATGCGAGTGTTGAAGTTGGGGGGTACATAAGTGGCAAATGCAGaa     >  2:1089475/1‑150 (MQ=255)
                                                                                                                                                   ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc    >  1:1507716/1‑150 (MQ=255)
                                                                                                                                                   ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc    >  1:257779/1‑150 (MQ=255)
                                                                                                                                                     acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATACAGAACGt  <  2:1483218/150‑1 (MQ=255)
                                                                                                                                                     acaGGTAACGGTGCGGGATGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTTACAACCATGCGAGTGTTTAAGTTGGGCGGTACATAAGTGCCAAAAGCAGAACGt  >  1:756242/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                        
CTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGT  >  NC_000913/86‑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: ‑