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.

GAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGT  >  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)
                                tAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑a                                                                                                                   <  1:1013869/150‑1 (MQ=255)
                                 aaaTACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑aa                                                                                                                  >  2:1359183/1‑150 (MQ=255)
                                     aCTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑aaaaaa                                                                                                              >  2:1467998/1‑150 (MQ=255)
                                     aCTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑aaaaaa                                                                                                              >  1:1178123/1‑150 (MQ=255)
                                      cTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAg                                                                                                             <  2:1178123/150‑1 (MQ=255)
                                       tttAACCAATATAGGCTTAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGc                                                                                                            <  2:160668/150‑1 (MQ=255)
                                        ttAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGcc                                                                                                           <  1:823668/150‑1 (MQ=255)
                                           aCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGc                                                                                                        >  2:1248978/1‑150 (MQ=255)
                                           aCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGc                                                                                                        <  2:1063421/150‑1 (MQ=255)
                                           aCCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGc                                                                                                        >  1:146419/1‑150 (MQ=255)
                                             cAATATAGGCATAGCGCCCAGACAGATAAAAAATACCGAGTACACAACATCCATGAAACGAAATACCACCAACATTACCACAACCATCACCATTTCCACAGGGAACGGTGCGGCATGACGCGAACAGGAAGCACAG‑AACAAAACCCCCAc                                                                                                      >  1:1142372/1‑150 (MQ=255)
                                                  taGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGAc                                                                                                 >  1:993145/1‑150 (MQ=255)
                                                     gCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑Gt                                                                                             >  1:180572/1‑150 (MQ=255)
                                                        tAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCg                                                                                          >  2:864179/1‑150 (MQ=255)
                                                         aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                         <  1:335805/150‑1 (MQ=255)
                                                          gcgcACAGACAGATAAAAATTACAGAGTACACAACATCGATGAAACGCATAAGCACCACCATTACCACCACCATCACCATTACCACAGGTATCGCGGGGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCggg                                                                                        <  2:1406653/150‑1 (MQ=255)
                                                          gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCggg                                                                                        <  2:1638288/150‑1 (MQ=255)
                                                           cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCA‑CTGTCAGGTGCGGGc                                                                                       <  2:1098441/150‑1 (MQ=255)
                                                           cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCA‑CTGTCAGGTGCGGGc                                                                                       <  1:1467998/150‑1 (MQ=255)
                                                            gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑t                                                                                     >  1:1422712/1‑150 (MQ=255)
                                                            gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑t                                                                                     <  1:1617580/150‑1 (MQ=255)
                                                             cacaGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑tt                                                                                    <  2:1037534/150‑1 (MQ=255)
                                                              acagacagATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑ttt                                                                                   <  1:705800/150‑1 (MQ=255)
                                                                agacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttt‑tt                                                                                <  1:1615715/150‑1 (MQ=255)
                                                                 gacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCCGACA‑GTGCGGGC‑ttt‑ttt                                                                               >  1:1468242/1‑150 (MQ=255)
                                                                  acagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGTTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttt‑tttt                                                                              <  1:864179/150‑1 (MQ=255)
                                                                   cagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCGGACA‑GTGCGGGC‑ttt‑ttttt                                                                             >  1:186766/1‑150 (MQ=255)
                                                                   attaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑ttt‑tttct                                                                             <  2:1662501/147‑3 (MQ=39)
                                                                    agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttt‑tttttt                                                                            <  2:1183107/150‑1 (MQ=255)
                                                                    agaTAAAAATTACAGAGGACTCAACATCCATGAACCGCATTAGCACCACCATAACCACCACCATCACCATAACCACAGGTAACGGTGCGGGCTGAAGCGTACAGGAAACACAGAAAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑ttt‑tttct                                                                             <  2:278534/150‑3 (MQ=37)
                                                                     gaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTc                                                                           >  2:875165/1‑150 (MQ=255)
                                                                      aTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑ttt‑tttctgtg                                                                          <  2:149740/150‑6 (MQ=39)
                                                                       tAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCgc                                                                         >  1:941855/1‑149 (MQ=255)
                                                                        aaaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGATACACCG‑AAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑ttt‑tttctgtgtt                                                                        <  1:271591/150‑8 (MQ=38)
                                                                        aaaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGAc                                                                        <  1:592041/150‑1 (MQ=255)
                                                                         aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATAACCACAGGTAACGGTGCGGGCTGACGGGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GAGCGGGC‑TTT‑TTTTTTCGagc                                                                       >  2:837152/1‑148 (MQ=255)
                                                                                cAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTa                                                                <  1:1346276/150‑1 (MQ=255)
                                                                                cAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTa                                                                <  1:1462508/150‑1 (MQ=255)
                                                                                     tACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCAACACCATTACCACAGGTAACGGTGCGGGCTGACGCGAACAGGAAACACAG‑AAAAAAGCCCGCACATGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGATAACCAg                                                           >  1:88768/1‑150 (MQ=255)
                                                                                       cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGt                                                         >  1:1481330/1‑150 (MQ=255)
                                                                                            caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTaacaa                                                    >  1:790084/1‑150 (MQ=255)
                                                                                            caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGGAACGAGGTaacaa                                                    >  2:852408/1‑150 (MQ=255)
                                                                                              aacATTAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAT‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAAcc                                                  <  2:1439670/148‑1 (MQ=255)
                                                                                               ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                 <  2:701122/150‑1 (MQ=255)
                                                                                                  tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCCACCAAAGGTACCGAGGTAACAACCATGc                                              >  1:313501/1‑150 (MQ=255)
                                                                                                     aaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                           <  1:178938/150‑1 (MQ=255)
                                                                                                       cGCATTAGCACCACCATTACCACCACCATCACCATTACCTCAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACTG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgtg                                         <  2:1408212/150‑1 (MQ=255)
                                                                                                            tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                    <  1:1665723/150‑1 (MQ=255)
                                                                                                             aGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAg                                   >  1:789747/1‑150 (MQ=255)
                                                                                                              gCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGt                                  >  2:1501786/1‑150 (MQ=255)
                                                                                                                accaccATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGCCGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                >  2:1525902/1‑150 (MQ=255)
                                                                                                                 caaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGATACACTG‑AAAAAAGGCGGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCATAGGTAACGAGGTAACAAGCATGCGAGTGATGAAGTTcg                               <  2:939580/148‑1 (MQ=255)
                                                                                                                  caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                              <  1:1680191/150‑1 (MQ=255)
                                                                                                                    ccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACCACCATGCGAGTGGTGAAGTTcggcg                            >  2:1205127/1‑150 (MQ=255)
                                                                                                                       ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                         <  1:1359/150‑1 (MQ=255)
                                                                                                                       ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                         <  1:1359183/150‑1 (MQ=255)
                                                                                                                       ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                         <  2:568944/150‑1 (MQ=255)
                                                                                                                        tACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAc                        <  1:794818/150‑1 (MQ=255)
                                                                                                                        tACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAc                        <  2:1422712/150‑1 (MQ=255)
                                                                                                                          ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAt                      >  1:870123/1‑150 (MQ=255)
                                                                                                                          caaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCACGCACCTGACA‑GTGGGGGC‑TGT‑TTTTTTCGACCAAAGGTAAAGAGGTAACACCACTGCGAGTTTTGAAGTTAGGCGGTACAt                      >  2:1280236/3‑150 (MQ=255)
                                                                                                                           caccaccaTAACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACAAAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                     <  2:88768/150‑1 (MQ=255)
                                                                                                                             ccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAg                   >  1:931054/1‑150 (MQ=255)
                                                                                                                                 ccgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑ATATTTCGACCAAAGGTAACGATGTAACAACCATGCGAGTGTTGAAGTTAGGCGGTACAACAGTGGc               >  2:589140/4‑150 (MQ=255)
                                                                                                                                 catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc               >  1:1573245/1‑150 (MQ=255)
                                                                                                                                  cgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa              >  2:1243236/3‑150 (MQ=255)
                                                                                                                                    caccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTGTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaa             >  1:1109522/1‑150 (MQ=255)
                                                                                                                                       caatACAACAGGTAACGGTCCGGGCTGACGGGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCATTGGCAAATg          <  2:878638/147‑1 (MQ=255)
                                                                                                                                        atTACCACAGGTAACGGTGCGGGCTGACGGGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa        <  1:1609719/150‑1 (MQ=255)
                                                                                                                                           aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa     >  1:762714/1‑150 (MQ=255)
                                                                                                                                           aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAAAAAACATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa     >  1:562923/1‑150 (MQ=255)
                                                                                                                                           aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGAAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATTAAGaa     >  1:762716/1‑150 (MQ=255)
                                                                                                                                           aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TGT‑TTTTTTCGACCAAAGGTAACGAGTGTACAACCATGCGAGTGTTGAAGTTGGGGGGTACATAAGTGGCAAATGCAGaa     >  2:1089475/1‑150 (MQ=255)
                                                                                                                                            ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc    >  1:1507716/1‑150 (MQ=255)
                                                                                                                                              acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATACAGAACGt  <  2:1483218/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                        
GAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTT‑TTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGT  >  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: ‑