Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,900,324 G→T 6.1% intergenic (‑51/‑268) ygcW ← / → yqcE putative SDR family oxidoreductase/putative MFS transporter, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,900,3240GT 6.1% 409.0 / 14.6 148intergenic (‑51/‑268)ygcW/yqcEputative SDR family oxidoreductase/putative MFS transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (63/76);  new base T (4/5);  total (67/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.07e-02

Only 100 of 151 total aligned reads displayed.

TGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACACAA  >  NC_000913/2900175‑2900473
                                                                                                                                                       |                                                                                                                                                     
tGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCcg                                                                                                                                                       >  2:1327205/1‑150 (MQ=255)
   cTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCga                                                                                                                                                    <  1:566715/150‑1 (MQ=255)
   cTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCga                                                                                                                                                    >  1:1743099/1‑150 (MQ=255)
   cTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCga                                                                                                                                                    <  1:1361462/150‑1 (MQ=255)
    tGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCgag                                                                                                                                                   >  2:323468/1‑150 (MQ=255)
    tGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCgag                                                                                                                                                   <  1:1500552/150‑1 (MQ=255)
    tGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCgag                                                                                                                                                   >  1:815946/1‑150 (MQ=255)
      gCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAGTCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCgaggg                                                                                                                                                 >  1:1903866/1‑148 (MQ=255)
      gCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAAAGAAAGGTAACTTCCCCCGCgag                                                                                                                                                   <  1:947835/150‑1 (MQ=255)
       ccTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAgg                                                                                                                                                >  2:272388/1‑150 (MQ=255)
       ccTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGCGAGAAATAATCCATTGAGAACGCATTGAGAGTTTCGTTTGACATAGTTAGCTCCAATTTCCGAAACGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAgg                                                                                                                                                <  2:71239/150‑1 (MQ=255)
         tAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGa                                                                                                                                              <  2:1012834/150‑1 (MQ=255)
          aaaCCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGaa                                                                                                                                             <  1:274660/150‑1 (MQ=255)
          aaaCCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGaa                                                                                                                                             <  1:1210452/150‑1 (MQ=255)
           aaCCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGt                                                                                                                                           >  1:897432/1‑150 (MQ=255)
                cTATTCCCACCTGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGGAAGGTAACTTCCCCCGCGAGAGGGAAGTTAcc                                                                                                                                       <  1:1669472/150‑1 (MQ=255)
                cTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACTCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCTCTCTCGGGGGAAGTTAcc                                                                                                                                       <  1:1438833/150‑1 (MQ=255)
                cTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCTCTCGCGGGGGAAGTTAcc                                                                                                                                       <  1:1577835/150‑1 (MQ=255)
                 taatCCCACAGTTAACAATTGCGGTTTTACCTTTCATGGAGAAAAACACCATTTAGAACGCATTGATAGATTCGCTTGACATAGTTAGCTCCAATTTCCGTAAGGA‑‑AAAAAAAAAGAAAGGTAACTTCCCTCTCGCTGGGGAAGTTACCt                                                                                                                                      <  1:990242/147‑1 (MQ=255)
                 tATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCTCTCGCGGGGGAAGTTACCt                                                                                                                                      <  1:1722654/150‑1 (MQ=255)
                      ccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTtctc                                                                                                                                 <  1:691999/150‑1 (MQ=255)
                         ccGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCAtt                                                                                                                              <  2:338400/150‑1 (MQ=255)
                          cGGTAACAATTGCGGTTTTAACTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTa                                                                                                                             <  2:185503/150‑1 (MQ=255)
                            gTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGAAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTAtc                                                                                                                           <  2:1399388/150‑1 (MQ=255)
                               aCAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATGGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAACAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATctct                                                                                                                        <  1:1511498/150‑1 (MQ=255)
                                cAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTg                                                                                                                       >  1:606439/1‑150 (MQ=255)
                                 aaTTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTgg                                                                                                                      >  1:417551/1‑150 (MQ=255)
                                   ttGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTa                                                                                                                    >  1:1793781/1‑150 (MQ=255)
                                    tGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATATATGGTaa                                                                                                                   >  1:830310/1‑150 (MQ=255)
                                       ggTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATATCTGGTAACTg                                                                                                                >  1:183733/1‑150 (MQ=255)
                                         ttttACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTgg                                                                                                               <  2:695076/150‑1 (MQ=255)
                                          tttACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGgag                                                                                                             >  2:465271/1‑150 (MQ=255)
                                               cTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAAt                                                                                                        >  2:22571/1‑150 (MQ=255)
                                                 ttCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAAttt                                                                                                      >  1:1981757/1‑150 (MQ=255)
                                                 ttCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCCAGTTACCAGAGATAAtga                                                                                                      <  2:1235462/150‑3 (MQ=255)
                                                   cAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAAttttt                                                                                                    <  2:525551/150‑1 (MQ=255)
                                                     gtgagaAAAAATACATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGTATGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAg                                                                                                  <  1:43722/148‑1 (MQ=255)
                                                     gggAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAg                                                                                                  <  1:1732987/150‑1 (MQ=255)
                                                     gggAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCCAGTTACCAGAGATAAtgagaaa                                                                                                  <  1:185904/150‑7 (MQ=255)
                                                       gagaAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAAATTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCa                                                                                                <  1:633331/150‑1 (MQ=255)
                                                             aaaTCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGc                                                                                         <  2:897432/150‑1 (MQ=255)
                                                                 ccATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATa                                                                                      <  2:1514752/150‑1 (MQ=255)
                                                                 ccATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATa                                                                                      <  2:1163066/150‑1 (MQ=255)
                                                                  cATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATaa                                                                                     >  2:1392322/1‑150 (MQ=255)
                                                                   aTTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTATCATTATATATGGTAACTGGAGATAATTTTTAGCATAAGGGCATAAc                                                                                    >  1:1243442/1‑150 (MQ=255)
                                                                       agaACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCt                                                                                <  1:323468/150‑1 (MQ=255)
                                                                           cGCATTTACAGATTAGATTGACATACTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAAtg                                                                            <  2:1332578/150‑1 (MQ=255)
                                                                           cGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAAtg                                                                            <  2:860308/150‑1 (MQ=255)
                                                                           cGCATTGAGAGATTCGATTGACATAGTTACCTCCAATTTCCGAAACGC‑‑AAAAAAAAAGAAAGGGCACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATACGGGCATAACTGCTAAtg                                                                            <  2:281230/150‑1 (MQ=255)
                                                                              aTTGAGAGATTCGATTGACATAGTTAGATCCAATTTCCGTAAGGC‑‑AAATAAAAAGAAAGGTAACTTCACCCGCTAGAGGGAAGTTACCTTTCTCATTATATATGGTTACTGGAGATAATTTTTAGAAGAAGGGCATAAATGATAATTTGa                                                                         >  1:129205/1‑150 (MQ=255)
                                                                               ttGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGAt                                                                        <  2:209468/150‑1 (MQ=255)
                                                                               ttGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGAt                                                                        <  2:1743099/150‑1 (MQ=255)
                                                                                tGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATc                                                                       <  2:770245/150‑1 (MQ=255)
                                                                                 gagagaTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCg                                                                      >  2:870786/1‑150 (MQ=255)
                                                                                     taTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAAt                                                                  <  2:1829662/149‑1 (MQ=255)
                                                                                         cGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAg                                                             <  2:568164/150‑1 (MQ=255)
                                                                                              gACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGtgtg                                                         <  2:1990713/150‑1 (MQ=255)
                                                                                              gACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGtgtg                                                         <  2:1723025/150‑1 (MQ=255)
                                                                                               aCATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGa                                                        <  2:68922/150‑1 (MQ=255)
                                                                                               aCATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGa                                                        <  2:417551/150‑1 (MQ=255)
                                                                                                cATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGat                                                       >  1:909765/1‑150 (MQ=255)
                                                                                                  tAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGatat                                                     >  1:1923191/1‑150 (MQ=255)
                                                                                                  tAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATAtt                                                    >  2:1927086/1‑150 (MQ=255)
                                                                                                          tCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTg                                             >  1:1012560/1‑150 (MQ=255)
                                                                                                           ccAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGt                                            >  1:1683442/1‑150 (MQ=255)
                                                                                                            cAATTTCCGTAAGGC‑‑AAAAAAAAAGAAATGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTc                                           <  2:1752195/150‑1 (MQ=255)
                                                                                                               tttCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCAAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCaaa                                        <  1:869625/150‑1 (MQ=255)
                                                                                                                 tCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAaata                                      >  1:330811/1‑150 (MQ=255)
                                                                                                                  ccGTAAGGCAAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAaata                                      >  2:112609/1‑150 (MQ=255)
                                                                                                                  ccGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAaataa                                     <  1:1128767/150‑1 (MQ=255)
                                                                                                                  ccGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCCTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAaataa                                     >  2:1673359/1‑150 (MQ=255)
                                                                                                                     tAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCAGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAAttt                                  <  1:125466/150‑1 (MQ=255)
                                                                                                                      aaGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGGGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAAtttt                                 <  2:1028538/150‑1 (MQ=255)
                                                                                                                       aGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTa                                <  1:1368400/150‑1 (MQ=255)
                                                                                                                        ggC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTaaa                               >  2:218476/1‑150 (MQ=255)
                                                                                                                         gC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTaaa                              <  1:1198410/150‑1 (MQ=255)
                                                                                                                         gC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTaaa                              >  2:535750/1‑150 (MQ=255)
                                                                                                                          c‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTaaaa                             >  2:1820812/1‑150 (MQ=255)
                                                                                                                              aaaaaaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAg                           >  2:1964132/1‑150 (MQ=255)
                                                                                                                               aaaaaaaGAAAGGTATCTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAAACACACTGTGATATTAGTATAGACAAATAATTTTAAAAAGc                          >  1:158772/1‑150 (MQ=255)
                                                                                                                               aaaaaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTAGCTTTCTCATTATCTCTGGTAACTGGAGATAATTTATAGCATAAGGGCATAAGAGCTAATGAGATCGTAATCACAGTGTGATATTCGTATTCTCAAATAATTTTAAAAAGc                          >  1:819863/1‑150 (MQ=255)
                                                                                                                                 aaaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAg                        <  2:398460/150‑1 (MQ=255)
                                                                                                                                   aaagaaagGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGaa                      >  1:1489402/1‑150 (MQ=255)
                                                                                                                                        aagGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTg                 >  2:488144/1‑150 (MQ=255)
                                                                                                                                         agGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGa                >  2:1663228/1‑150 (MQ=255)
                                                                                                                                          ggTAACTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAt               >  1:1719929/1‑150 (MQ=255)
                                                                                                                                          ggTAACTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAt               >  1:228434/1‑150 (MQ=255)
                                                                                                                                          ggTAACTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAt               >  2:1478946/1‑150 (MQ=255)
                                                                                                                                          ggTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAt               >  1:1516815/1‑150 (MQ=255)
                                                                                                                                           gTAACTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAtt              >  1:15317/1‑150 (MQ=255)
                                                                                                                                           gTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTAATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATAGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAtt              <  1:1244912/150‑1 (MQ=255)
                                                                                                                                            tAACTTCCCCCGCGAGAGTGAAGTTAACTTTCTCATTATCTTTGGTAACTTGATATAATTTTTAGAATAAGGTCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTa             <  1:1240864/150‑1 (MQ=255)
                                                                                                                                            tAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGTTAACTGGAGATAAGTTTTAGCAGAAGGGCATAACTGATAATGTGATCGTAATCACAGTGTGAGATTCGTATTGTCAAATAATGTTAAAAAGCAGAAGCGTGATTa             >  1:1974122/1‑150 (MQ=255)
                                                                                                                                              aCTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTAct           >  2:1343508/1‑150 (MQ=255)
                                                                                                                                                ttCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTActct         <  2:606439/150‑1 (MQ=255)
                                                                                                                                                 tCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACtctc        <  1:1820812/150‑1 (MQ=255)
                                                                                                                                                  cccccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTca       >  1:1943988/1‑150 (MQ=255)
                                                                                                                                                   ccccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTcac      <  1:535255/150‑1 (MQ=255)
                                                                                                                                                   ccccGCGAGAGGGAAGTGACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTAAAATAATTTTAAAAAGAAGAAGCGTGATTACTCTcac      >  1:1406931/1‑150 (MQ=255)
                                                                                                                                                       gcgAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATTTGATCGTAATCAAAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAATCGTGATTACTCTCAAacaa  >  1:317823/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
TGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGC‑‑AAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACACAA  >  NC_000913/2900175‑2900473

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