Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,592,201 T→A 7.5% V182D (GTC→GAC)  tsr → methyl‑accepting chemotaxis protein I, serine sensor receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,592,2010TA 7.5% 307.1 / 11.9 121V182D (GTC→GAC) tsrmethyl‑accepting chemotaxis protein I, serine sensor receptor
Reads supporting (aligned to +/- strand):  ref base T (60/50);  new base A (6/3);  total (66/55)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

Only 100 of 121 total aligned reads displayed.

TGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCAT  >  NC_000913/4592058‑4592349
                                                                                                                                               |                                                                                                                                                    
tGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGa                                                                                                                                                <  2:606003/150‑1 (MQ=255)
 gcgcTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGAt                                                                                                                                               <  1:824009/150‑1 (MQ=255)
    cTGGCGGAGCTGATCCAACTTTTTGGTGTAGGCAAAATCAACTAGTTCTTTGATCAGCGGACCCAGGTATATGAGGACGGTTTCGTGAAGCAGTATGTGGCTTACATGGAGGAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAAc                                                                                                                                            <  2:318686/150‑1 (MQ=255)
      ggCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGTAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACaa                                                                                                                                          >  2:1172668/1‑150 (MQ=255)
        cGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGCCGGTTTCGAGCAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATg                                                                                                                                        <  2:112910/150‑1 (MQ=255)
        cGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATg                                                                                                                                        <  2:981340/150‑1 (MQ=255)
        cGGAGCTGATCCAACTGGTAGGGGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATg                                                                                                                                        <  2:668005/150‑1 (MQ=255)
         ggAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCTGACGGCGATATCatgga                                                                                                                                       <  2:262122/150‑5 (MQ=255)
          gAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGcc                                                                                                                                      <  1:1089251/150‑1 (MQ=255)
          gAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGcc                                                                                                                                      <  1:658117/150‑1 (MQ=255)
             cTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGcctcc                                                                                                                                   <  1:1425237/150‑1 (MQ=255)
             cTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGcctcc                                                                                                                                   >  2:1421427/1‑150 (MQ=255)
                 tCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACa                                                                                                                               >  1:1336497/1‑150 (MQ=255)
                   cAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCTGACGGCGATATCatgg                                                                                                                                        <  1:1230876/139‑4 (MQ=255)
                   cAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGc                                                                                                                             >  2:394152/1‑150 (MQ=255)
                     aCTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCa                                                                                                                           >  2:952056/1‑150 (MQ=255)
                       tGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAgg                                                                                                                         >  1:1581788/1‑150 (MQ=255)
                          tAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGa                                                                                                                      >  1:1411510/1‑150 (MQ=255)
                            ggTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGAtg                                                                                                                    >  2:574424/1‑150 (MQ=255)
                              tGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGAtgtg                                                                                                                  >  2:1613861/1‑150 (MQ=255)
                               gCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTgg                                                                                                                 <  2:1458917/150‑1 (MQ=255)
                                   gCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTc                                                                                                             <  1:1156844/150‑1 (MQ=255)
                                    cAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCt                                                                                                            <  2:1018650/150‑1 (MQ=255)
                                    cAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCt                                                                                                            >  1:971586/1‑150 (MQ=255)
                                       aaTCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCtggt                                                                                                         >  1:1354538/1‑150 (MQ=255)
                                        aTCAACGAGTTCTTTGATCTGCCGACCCACGGATATCAGGACGGTGTCGAGAAGCAGTAAGTGGCATACATGGAGCAAACCGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCtggtg                                                                                                        <  2:299733/150‑1 (MQ=255)
                                        aTCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCtggtg                                                                                                        <  1:1113791/150‑1 (MQ=255)
                                          cAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTggg                                                                                                      <  1:59005/150‑1 (MQ=255)
                                           aaCGCGTGCTTTTATCAGCCGCCCCAGGGTTAAAAGGACGGTTTCGAGAAGCAGTATTTGGCTTACATGGAGCAACAGGAGCGGCTCCATGATATCGCCTTCAGCGATAACAATGCCTCCTCCAGCGAGGCGATGTGGATTCTGGTGGGc                                                                                                     <  2:1158763/150‑1 (MQ=255)
                                             cGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGt                                                                                                   >  1:1553168/1‑150 (MQ=255)
                                              gAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGtg                                                                                                  >  1:777849/1‑150 (MQ=255)
                                               aGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGtga                                                                                                 >  1:1579706/1‑150 (MQ=255)
                                                gTTCTTTGATCTGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACAGGGAGCAAAACGATCGGCGCCATGATATCGCCGGCAGCGATAACAATGGCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGtgat                                                                                                <  1:1147699/150‑1 (MQ=255)
                                                 ttCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGtgatg                                                                                               <  2:565058/150‑1 (MQ=255)
                                                  tCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGtgatga                                                                                              >  1:550835/1‑150 (MQ=255)
                                                    tttGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGAtc                                                                                            >  1:252632/1‑150 (MQ=255)
                                                      tGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGAtcgt                                                                                          >  1:518560/1‑150 (MQ=255)
                                                       gATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGAtcgtc                                                                                         >  2:137766/1‑150 (MQ=255)
                                                       gATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGAtcgtc                                                                                         >  2:1220721/1‑150 (MQ=255)
                                                        aTCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGAtcgtcg                                                                                        <  1:137766/150‑1 (MQ=255)
                                                          cAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGAccgccagt                                                                                      >  2:498477/1‑145 (MQ=255)
                                                             ccGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTg                                                                                   >  1:87143/1‑150 (MQ=255)
                                                                aCCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCg                                                                                >  2:553736/1‑150 (MQ=255)
                                                                 cccAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCgg                                                                               >  1:463537/1‑150 (MQ=255)
                                                                 cccAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCgg                                                                               >  1:923500/1‑150 (MQ=255)
                                                                  ccAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGt                                                                              <  2:1258270/150‑1 (MQ=255)
                                                                    aGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCa                                                                            <  1:1448809/150‑1 (MQ=255)
                                                                    aGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCa                                                                            <  1:1456444/150‑1 (MQ=255)
                                                                       gATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGAACGTCGTACTGGCGGACATCt                                                                         >  1:344991/1‑150 (MQ=255)
                                                                          atCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCg                                                                      <  2:971586/150‑1 (MQ=255)
                                                                            cAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGcc                                                                    >  2:960958/1‑150 (MQ=255)
                                                                                   gTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGt                                                             >  1:1049103/1‑150 (MQ=255)
                                                                                    tttCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGtt                                                            >  1:1699961/1‑150 (MQ=255)
                                                                                     ttCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTc                                                           >  1:914247/1‑150 (MQ=255)
                                                                                     ttCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTc                                                           <  2:518560/150‑1 (MQ=255)
                                                                                       cGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGTTGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCgg                                                         >  1:42325/1‑150 (MQ=255)
                                                                                       cGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCgg                                                         <  2:252632/150‑1 (MQ=255)
                                                                                        gagaAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGt                                                        <  2:1336497/150‑1 (MQ=255)
                                                                                         agaAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTa                                                       >  2:1233793/1‑150 (MQ=255)
                                                                                           aaGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGTCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTAtt                                                     <  1:20691/150‑1 (MQ=255)
                                                                                            aGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTa                                                    >  2:1272042/1‑150 (MQ=255)
                                                                                             gCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTaa                                                   <  2:1049103/150‑1 (MQ=255)
                                                                                             gCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCACCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTaa                                                   >  2:702145/1‑150 (MQ=255)
                                                                                              cAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTaaa                                                  >  2:572337/1‑150 (MQ=255)
                                                                                              cAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTAAGCGATAACAATTCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTaaa                                                  >  2:158427/1‑150 (MQ=255)
                                                                                               aGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAg                                                 <  2:1581788/150‑1 (MQ=255)
                                                                                                gTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGc                                                <  1:1613861/150‑1 (MQ=255)
                                                                                                   tgtgGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACAGGCGGTCATCTTCGCCGTCTGGGTCGGTATTTAAGCCTc                                             >  1:768799/1‑150 (MQ=255)
                                                                                                     tgGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGc                                           >  2:609026/1‑150 (MQ=255)
                                                                                                     tgGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGc                                           <  1:574424/150‑1 (MQ=255)
                                                                                                      ggCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCt                                          >  2:278366/1‑150 (MQ=255)
                                                                                                      ggCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCt                                          <  2:1366727/150‑1 (MQ=255)
                                                                                                        cTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTgg                                        >  2:435929/1‑150 (MQ=255)
                                                                                                            cATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAgc                                    >  2:1126786/1‑150 (MQ=255)
                                                                                                              tGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTgg                                        >  1:1377156/1‑144 (MQ=255)
                                                                                                              tGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTgg                                        <  2:1377156/144‑1 (MQ=255)
                                                                                                              tGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAgcgc                                  >  1:878490/1‑150 (MQ=255)
                                                                                                              tGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAgcgc                                  >  2:483648/1‑150 (MQ=255)
                                                                                                                   cAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTTGCTGGTAGCGCCAATg                             >  1:746644/1‑150 (MQ=255)
                                                                                                                   cAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATg                             >  2:1452663/1‑150 (MQ=255)
                                                                                                                    aaaaCGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGa                            >  1:247883/1‑150 (MQ=255)
                                                                                                                    aaaaCGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGa                            >  2:417578/1‑150 (MQ=255)
                                                                                                                        cGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCg                        <  1:1220721/150‑1 (MQ=255)
                                                                                                                        cGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCg                        >  2:1035067/1‑150 (MQ=255)
                                                                                                                         gATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGc                       <  1:779380/150‑1 (MQ=255)
                                                                                                                           tCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCt                     <  1:1035067/150‑1 (MQ=255)
                                                                                                                                tCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTg                >  1:534429/1‑150 (MQ=255)
                                                                                                                                tCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTACCGCCAATGAATCGCCTGATTg                >  2:622863/1‑150 (MQ=255)
                                                                                                                                gcattgttatCGCTGACGGCGATATCAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTg                >  1:1570149/8‑150 (MQ=255)
                                                                                                                                gcattgttatCGCTGACGGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTg                >  2:238742/8‑150 (MQ=255)
                                                                                                                                gcattgttatCGCTGACGGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTg                >  2:579241/8‑150 (MQ=255)
                                                                                                                                  caggATATCGCCGGCAGCGATAACTATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGGGAGGATCGTCGTACTGGCGGTCATCTTCGCCGGCTGGTGCGGGATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGAc              <  1:224763/147‑1 (MQ=255)
                                                                                                                                   ttgttatCGCTGACGGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACa             >  2:1014976/5‑150 (MQ=255)
                                                                                                                                   aTGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGAcc             >  1:1591609/1‑149 (MQ=255)
                                                                                                                                     gATATCGCCGTCAGCGATAACAATGCCTCCTCCAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGCTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGc           <  2:465929/150‑1 (MQ=255)
                                                                                                                                      ttatCGCTGACGGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCa          >  1:100928/2‑150 (MQ=255)
                                                                                                                                       tatCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCAt         <  1:285899/150‑1 (MQ=255)
                                                                                                                                         tCGCCGTTAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTc       <  1:1096170/150‑1 (MQ=255)
                                                                                                                                            ccGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTc    <  2:1671177/150‑1 (MQ=255)
                                                                                                                                             cGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCa   <  2:1477739/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCAT  >  NC_000913/4592058‑4592349

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: