Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,367,342 C→T 6.7% E191K (GAA→AAA) ‡ pspF ← psp operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,367,3420CT 6.7% 231.0 / 4.7 74E191K (GAA→AAA) ‡pspFpsp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (32/37);  new base T (3/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01

TGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACA  >  NC_000913/1367194‑1367490
                                                                                                                                                    |                                                                                                                                                          
tGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTc                                                                                                                                                            <  1:1568732/149‑1 (MQ=255)
 ggCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCt                                                                                                                                                           <  1:1516545/149‑1 (MQ=255)
  gCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTg                                                                                                                                                          >  1:604252/1‑149 (MQ=255)
  gCGATACACTGAACGTTCCACCACGTTTTTCAATTAAAGAAAATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGAGGGCAAAGTATTcgg                                                                                                                                                          >  1:811361/1‑147 (MQ=255)
   cGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGc                                                                                                                                                         >  1:1064317/1‑149 (MQ=255)
    gATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGcc                                                                                                                                                        >  1:503694/1‑149 (MQ=255)
     aTACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCTCTCCGTAAACCCCGGGCACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCa                                                                                                                                                       >  1:334259/1‑149 (MQ=255)
         acTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCaa                                                                                                                                                   >  1:457395/1‑149 (MQ=255)
         acTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCaa                                                                                                                                                   <  1:221793/149‑1 (MQ=255)
          cTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAAc                                                                                                                                                  >  1:12177/1‑149 (MQ=255)
            gAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAGAATACTTTGCCATccagat                                                                                                                                                <  1:1588318/149‑6 (MQ=255)
             aaCGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACAtt                                                                                                                                               <  1:140059/149‑1 (MQ=255)
             aaCGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACAtt                                                                                                                                               <  2:1447777/149‑1 (MQ=255)
             aaCGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACAtt                                                                                                                                               >  1:172405/1‑149 (MQ=255)
               cGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTAt                                                                                                                                             >  2:497377/1‑149 (MQ=255)
                gTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTAtt                                                                                                                                            >  1:893800/1‑148 (MQ=255)
                gTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATg                                                                                                                                            >  1:68430/1‑149 (MQ=255)
                gTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATg                                                                                                                                            >  1:935135/1‑149 (MQ=255)
                  tCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTc                                                                                                                                          <  2:1064317/149‑1 (MQ=255)
                      ccacGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCtt                                                                                                                                      >  2:843346/1‑149 (MQ=255)
                       cacGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCttt                                                                                                                                     >  1:1477356/1‑149 (MQ=255)
                              ttCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCCTCAACATTATTTCGATTtcgcgctc                                                                                                                              >  2:186931/1‑149 (MQ=255)
                                    tCACGAATATTTCCCGGCCTACGATAATTCAGCTATGTATCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAg                                                                                                                        <  1:866790/149‑1 (MQ=255)
                                    tCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAg                                                                                                                        >  2:910405/1‑149 (MQ=255)
                                        gAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCTCGCAgtggt                                                                                                                    >  1:705803/1‑149 (MQ=255)
                                        gAAAATTTCCCGGCCAACGATAATTCAGCAATGATTCTCTGGCGCGCTCAGTAAACACCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAgtggt                                                                                                                    >  1:201967/1‑149 (MQ=255)
                                         aaTATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAgtggtg                                                                                                                   <  2:604252/149‑1 (MQ=255)
                                             tttCCCGGCCAACGATATTTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAg                                                                                                               <  1:1577113/149‑1 (MQ=255)
                                             tttCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAg                                                                                                               >  1:915029/1‑149 (MQ=255)
                                                  cGGCCAACGATATTTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTa                                                                                                          <  2:412200/149‑1 (MQ=255)
                                                  cGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTa                                                                                                          <  2:935135/149‑1 (MQ=255)
                                                   ggCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCGGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTac                                                                                                         <  1:246935/149‑1 (MQ=255)
                                                   ggCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTac                                                                                                         <  2:43896/149‑1 (MQ=255)
                                                    gCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTaca                                                                                                        >  2:579385/1‑149 (MQ=255)
                                                      cAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGAGACATCTGGATGGGAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTagaac                                                                                                      <  1:165856/149‑1 (MQ=255)
                                                      cAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTacaac                                                                                                      <  1:718919/149‑1 (MQ=255)
                                                       aaCGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTacaaca                                                                                                     >  1:898052/1‑149 (MQ=255)
                                                        aCGATAATTCATCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACAt                                                                                                    <  2:941127/149‑1 (MQ=255)
                                                        aCGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACAt                                                                                                    <  2:941119/149‑1 (MQ=255)
                                                        aCGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACAt                                                                                                    <  2:946996/149‑1 (MQ=255)
                                                        aCGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGAAGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACAt                                                                                                    <  2:946999/149‑1 (MQ=255)
                                                             aaTTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCaaaa                                                                                               <  1:1190/149‑1 (MQ=255)
                                                              aTTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAg                                                                                              <  2:834763/149‑1 (MQ=255)
                                                               ttCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGc                                                                                             <  2:1486801/149‑1 (MQ=255)
                                                                     aaTGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑gtcg                                                                                 <  2:48236/149‑1 (MQ=255)
                                                                           tctcTGGCTCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATTGCTAAGTATTCTTCCATCAACATAATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAg                                                                           <  1:256388/149‑1 (MQ=255)
                                                                              ctGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAg                                                                        <  2:854509/149‑1 (MQ=255)
                                                                                ggCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGt                                                                      >  2:150055/1‑149 (MQ=255)
                                                                                         cGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGaa                                                             >  2:212007/1‑149 (MQ=255)
                                                                                           tAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGaaaa                                                           <  2:12177/149‑1 (MQ=255)
                                                                                             aaCCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGGGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGt                                                         <  1:768088/149‑1 (MQ=255)
                                                                                              aCCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTg                                                        <  2:898052/149‑1 (MQ=255)
                                                                                               ccccGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGGGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGc                                                       <  2:334259/149‑1 (MQ=255)
                                                                                                cccGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGcc                                                      <  1:459558/149‑1 (MQ=255)
                                                                                                             gCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCTCGAAAAGTGCCTTCATTGACCATc                                         <  2:915029/149‑1 (MQ=255)
                                                                                                               aGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGc                                       <  2:201967/149‑1 (MQ=255)
                                                                                                                    gATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGga                                  <  2:503694/149‑1 (MQ=255)
                                                                                                                    gATTTCCCGACACATCTGGAAGGCAAAGTATTCTGCCATCTACATTATGTCTCTTTCGCGCTCGCGCAGTGGTGGCATTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGga                                  <  1:1478789/149‑1 (MQ=255)
                                                                                                                       ttCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGGGGTGGCAGTTGTACAACATCAAACGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGAt                               <  1:443940/149‑1 (MQ=255)
                                                                                                                          ccGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCgg                            >  2:802826/1‑149 (MQ=255)
                                                                                                                           cGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGc                           >  1:1578612/1‑149 (MQ=255)
                                                                                                                             acacATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCAt                         >  2:1232895/1‑149 (MQ=255)
                                                                                                                                  tCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATtcgtc                    >  2:1558472/1‑149 (MQ=255)
                                                                                                                                    ttgATGGCAGAATACTTTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTcgc                  >  1:1080375/3‑149 (MQ=255)
                                                                                                                                    ttgATGGCAGAATACTTTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTcgc                  <  2:1080375/147‑1 (MQ=255)
                                                                                                                                    ttgATGGCAGAATACTTTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTcgc                  >  1:811319/3‑149 (MQ=255)
                                                                                                                                     ggATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTcgcg                 <  1:312189/149‑1 (MQ=255)
                                                                                                                                      gATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCgcgc                >  2:1018700/1‑149 (MQ=255)
                                                                                                                                      gATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTCTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCgcgc                >  1:356627/1‑149 (MQ=255)
                                                                                                                                       aTGGCAGAATACTTTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCa               >  2:1335861/1‑149 (MQ=255)
                                                                                                                                         ggCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGATAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATa             <  2:242446/149‑1 (MQ=255)
                                                                                                                                         ggCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATa             >  2:1366027/1‑149 (MQ=255)
                                                                                                                                           cAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATAcc           >  2:203492/1‑149 (MQ=255)
                                                                                                                                           cAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATAcc           >  2:1551379/1‑149 (MQ=255)
                                                                                                                                            aaaGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCCGACCa          >  1:1006837/1‑149 (MQ=255)
                                                                                                                                                tATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCg      <  1:843346/149‑1 (MQ=255)
                                                                                                                                                    cTGCCATCAACATTATGTCGCTTTCTCGCTCGCGCAGTGGGGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGcaca  <  2:68430/149‑1 (MQ=255)
                                                                                                                                                    cTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCACCACTGCGCGAGCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAAc        >  1:823139/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                          
TGGCGATACACTGAACGTTCCACCACGTTTTTCAATTCACGAATATTTCCCGGCCAACGATAATTCAGCAATGTTTCTCTGGCGCGCTCCGTAAACCCCGGGAACAGAGGCAGCTTGATTTCCCGACACATCTGGATGGCAAAGTATTCTGCCATCAACATTATGTCGCTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCA‑‑‑‑‑GTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACA  >  NC_000913/1367194‑1367490

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