Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,636,286 G→T 6.1% H357Q (CAC→CAA der ← 50S ribosomal subunit stability factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,636,2860GT 6.1% 471.0 / 9.3 140H357Q (CAC→CAAder50S ribosomal subunit stability factor
Reads supporting (aligned to +/- strand):  ref base G (67/63);  new base T (6/4);  total (73/67)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

AGCTTCACACGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAAT  >  NC_000913/2636137‑2636432
                                                                                                                                                     |                                                                                                                                                  
aGCTTCACACGACGACCGCGTACCAGCTGTGTTTTGTTATCTTCAACATCCATCGTGATGATGCGCGTTAGAATATATTTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCg                                                                                                                                                    <  2:1230196/150‑1 (MQ=255)
aGCTTCACACGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCg                                                                                                                                                    <  2:1280734/150‑1 (MQ=255)
 gCTTCACACGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCCAACAAGTTACCAACACCACTGCCGt                                                                                                                                                   >  2:2138308/1‑150 (MQ=255)
   ttCACACGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGc                                                                                                                                                 >  2:982821/1‑150 (MQ=255)
   ttCACACGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGc                                                                                                                                                 <  1:659993/150‑1 (MQ=255)
   ttCACACGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGc                                                                                                                                                 <  2:1526502/150‑1 (MQ=255)
       cacGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAgg                                                                                                                                             <  1:84630/150‑1 (MQ=255)
         cgacgaCCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCa                                                                                                                                           >  1:226227/1‑150 (MQ=255)
          gacgacCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCac                                                                                                                                          >  1:1071894/1‑149 (MQ=255)
           acgacCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCaga                                                                                                                                         <  1:2429811/150‑1 (MQ=255)
            cgacCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCagag                                                                                                                                        >  1:2420540/1‑150 (MQ=255)
             gacCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAgaga                                                                                                                                       <  2:2090082/150‑1 (MQ=255)
             gacCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAgaga                                                                                                                                       <  2:7660/150‑1 (MQ=255)
              acCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCCGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTAAAACAAGTTACCAACACCACTGCCTTGAAAGCCAgagaa                                                                                                                                      >  2:2395483/1‑149 (MQ=255)
              acCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCTTGCACGGCagtggt                                                                                                                                      <  1:1306357/150‑5 (MQ=255)
               ccGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATa                                                                                                                                     >  2:1474196/1‑150 (MQ=255)
                cgcgTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATaa                                                                                                                                    <  2:5977/150‑1 (MQ=255)
                      ccAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGca                                                                                                                              >  1:1121864/1‑150 (MQ=255)
                      ccAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGca                                                                                                                              <  1:164490/150‑1 (MQ=255)
                      ccAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGca                                                                                                                              >  1:995663/1‑150 (MQ=255)
                        aGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGcaca                                                                                                                            >  1:1622380/1‑150 (MQ=255)
                         gcggcgGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCacac                                                                                                                           <  2:405660/150‑1 (MQ=255)
                            gcggTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCCAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAg                                                                                                                        <  1:2138308/150‑1 (MQ=255)
                            gcggTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAg                                                                                                                        <  2:1498114/150‑1 (MQ=255)
                             cggTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGc                                                                                                                       <  2:1464885/150‑1 (MQ=255)
                                 tGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCaaaa                                                                                                                   <  2:152723/150‑1 (MQ=255)
                                   gTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATc                                                                                                                 >  2:1807523/1‑150 (MQ=255)
                                     gATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGa                                                                                                               <  1:1901607/150‑1 (MQ=255)
                                      aTCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGAt                                                                                                              <  2:1829393/150‑1 (MQ=255)
                                        cTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATaa                                                                                                            >  2:222492/1‑150 (MQ=255)
                                             aCAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGccc                                                                                                       >  2:2150538/1‑150 (MQ=255)
                                              cAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGca                                                                                                                              <  2:320162/126‑1 (MQ=255)
                                              cAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGca                                                                                                                              >  1:320162/1‑126 (MQ=255)
                                              cAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCa                                                                                                      >  2:1186558/1‑150 (MQ=255)
                                                gCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGa                                                                                                    >  1:1275879/1‑150 (MQ=255)
                                                  cATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACg                                                                                                  <  2:2267102/150‑1 (MQ=255)
                                                    tagTCATGATGGGGGTGTTCATAGTGGTCCCCAGACGACGGGTGGAGCTGTCATACGCTTGACTTTCTGATTAAAACAAGTTACCTACACCACTGCCTTGAAAGGCAGAGATAAATTGCACACTATGAAATTCGTTAAAGACCAGACGGa                                                                                                <  2:514359/148‑1 (MQ=255)
                                                    tCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGa                                                                                                <  1:1474196/150‑1 (MQ=255)
                                                     cGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGaa                                                                                               <  1:1843002/150‑1 (MQ=255)
                                                      gTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAg                                                                                              <  1:481126/150‑1 (MQ=255)
                                                        cATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTc                                                                                            <  1:453724/150‑1 (MQ=255)
                                                            atgCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCTGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGc                                                                                        >  2:1272546/1‑150 (MQ=255)
                                                            atgCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGc                                                                                        >  2:2277155/1‑150 (MQ=255)
                                                               cgcgTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGtt                                                                                     >  1:1243244/1‑150 (MQ=255)
                                                                   tCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGtttctt                                                                                 >  1:234065/1‑150 (MQ=255)
                                                                     aGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGtttctttc                                                                               <  1:829231/150‑1 (MQ=255)
                                                                     aGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGtttctttc                                                                               >  1:1806317/1‑150 (MQ=255)
                                                                      gCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCa                                                                              >  1:609680/1‑150 (MQ=255)
                                                                      gCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCa                                                                              >  1:901899/1‑150 (MQ=255)
                                                                       cATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCAc                                                                             >  2:2133165/1‑150 (MQ=255)
                                                                       cATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCAc                                                                             >  2:156950/1‑150 (MQ=255)
                                                                        aTAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCAcc                                                                            <  1:185073/150‑1 (MQ=255)
                                                                         tAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCt                                                                           <  2:1889023/150‑1 (MQ=255)
                                                                             ggTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGctc                                                                       >  1:1808426/1‑150 (MQ=255)
                                                                             ggTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGctc                                                                       >  1:333615/1‑150 (MQ=255)
                                                                                ccccACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCttt                                                                    >  1:108796/1‑150 (MQ=255)
                                                                                ccccACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCttt                                                                    >  2:6318/1‑150 (MQ=255)
                                                                                 cccACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTc                                                                   <  1:394385/150‑1 (MQ=255)
                                                                                  ccACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTGCAGCGTTTCTTTCACCTGCTCTTTCa                                                                  >  2:870375/1‑150 (MQ=255)
                                                                                   cacaCGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCAc                                                                 <  1:659773/150‑1 (MQ=255)
                                                                                       cgacgGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTcc                                                             <  1:2111152/150‑1 (MQ=255)
                                                                                       cgacgGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTcc                                                             <  1:272631/150‑1 (MQ=255)
                                                                                               ggAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAg                                                     >  1:924795/1‑150 (MQ=255)
                                                                                               ggAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAg                                                     <  2:353959/150‑1 (MQ=255)
                                                                                                  gCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGcc                                                  <  2:473050/150‑1 (MQ=255)
                                                                                                     gTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCa                                                      <  2:1632210/143‑1 (MQ=255)
                                                                                                       cATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCCCACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGCCTCAGGCCATccc                                             >  2:2078735/1‑150 (MQ=255)
                                                                                                       cATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATccc                                             <  1:1186558/150‑1 (MQ=255)
                                                                                                        aTACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCa                                            >  1:19061/1‑150 (MQ=255)
                                                                                                        aTACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCa                                            >  1:1275900/1‑150 (MQ=255)
                                                                                                        aTACGCTTCACGTACTGATTCAAACAAGTTACCAACACAACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCa                                            >  1:280676/1‑150 (MQ=255)
                                                                                                           cGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTATCACCTGCTCTTTCACTTCCGGACTCAGGCCATCCCACtt                                         >  2:2231941/1‑150 (MQ=255)
                                                                                                               tCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTc                                     <  2:1383973/150‑1 (MQ=255)
                                                                                                                cACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTca                                    <  2:454659/150‑1 (MQ=255)
                                                                                                                cACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTca                                    <  2:484251/150‑1 (MQ=255)
                                                                                                                cACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTca                                    >  1:1105013/1‑150 (MQ=255)
                                                                                                                 aCGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACTGTAGTCGAGAGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCTGACTTATTcac                                   <  2:1286565/150‑1 (MQ=255)
                                                                                                                    tACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTcaccac                                <  1:1684171/150‑1 (MQ=255)
                                                                                                                        gATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATg                            >  1:811935/1‑150 (MQ=255)
                                                                                                                        gATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATg                            >  2:1067080/1‑150 (MQ=255)
                                                                                                                        gATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATg                            <  1:2387828/150‑1 (MQ=255)
                                                                                                                        gATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATg                            <  2:324747/150‑1 (MQ=255)
                                                                                                                           tCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACa                         <  2:272759/150‑1 (MQ=255)
                                                                                                                           tCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACa                         <  1:214049/150‑1 (MQ=255)
                                                                                                                             aaaCAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAg                       >  1:559980/1‑150 (MQ=255)
                                                                                                                              aaCAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGt                      >  2:130197/1‑150 (MQ=255)
                                                                                                                               aCAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTg                     <  2:1217468/150‑1 (MQ=255)
                                                                                                                                    ttAGCAACACCAATGCCTTGCAAGGCAGAGATAAATTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAgcgc                <  2:1828524/150‑1 (MQ=255)
                                                                                                                                    ttACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCGCCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAgcgc                <  2:1243244/150‑1 (MQ=255)
                                                                                                                                     tACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGcc               >  1:288147/1‑150 (MQ=255)
                                                                                                                                       ccAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCa             >  2:281715/1‑150 (MQ=255)
                                                                                                                                        cAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCAc            >  1:1871437/1‑150 (MQ=255)
                                                                                                                                         tttatctcTGCCTTGCACGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACt           >  2:240511/8‑150 (MQ=255)
                                                                                                                                           caccacTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTAt         >  1:2004577/1‑150 (MQ=255)
                                                                                                                                            atctcTGCCTTGCACGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTAtt        >  2:1912883/5‑150 (MQ=255)
                                                                                                                                            atctcTGCCTTGCACGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTAtt        >  2:13784/5‑150 (MQ=255)
                                                                                                                                            atctcTGCCTTGCACGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTAtt        >  2:908545/5‑150 (MQ=255)
                                                                                                                                             tctcTGCCTTGCACGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGCCTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTg       >  1:1686270/4‑150 (MQ=255)
                                                                                                                                               acTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTgag     <  1:130197/150‑1 (MQ=255)
                                                                                                                                                  gCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGGTTCTTTCACCTGATCTTTCACTTCCTGACTCAGGCCATCCCAATTATGCACCACAATGACAAGTGAGCGCCCACTATTgagagt  >  2:238087/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AGCTTCACACGACGACCGCGTACCAGCGGCGGTTGGTGATCTTCAACAGCCATCGTCATGATGCGCGTCAGCATAGAGGTCCCCACACGACGGGTGGAGCTGTCATACGCTTCACGTACTGATTCAAACAAGTTACCAACACCACTGCCGTGCAAGGCAGAGATAAAGTGCACACGAGCAAAATCGATAAAGCCCAGACGGAAGTCCAGCGTTTCTTTCACCTGCTCTTTCACTTCCTGACTCAGGCCATCCCACTTATTCACCACAATGACAAGTGAGCGCCCACTATTGAGAAT  >  NC_000913/2636137‑2636432

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