Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,291,110 A→T 7.6% K766I (AAA→ATA) ‡ yraJ → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,291,1100AT 7.6% 226.6 / 7.1 88K766I (AAA→ATA) ‡yraJputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (32/47);  new base T (4/3);  total (37/51)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CTCGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGT  >  NC_000913/3290961‑3291258
                                                                                                                                                     |                                                                                                                                                    
ctcGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTaa                                                                                                                                                      >  2:1663726/1‑150 (MQ=255)
        aTACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGc                                                                                                                                              <  1:875111/150‑1 (MQ=255)
         tactacGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGcc                                                                                                                                             >  1:1459945/1‑150 (MQ=255)
           ctacGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTTTACCGTTGcgat                                                                                                                                           <  1:728314/150‑4 (MQ=255)
           ctacGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCt                                                                                                                                           >  2:285613/1‑150 (MQ=255)
            aacGCAGCTGCCCGATAACGGCGATCTTGAACAACCAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGGTACAGTCAGCGATCGCAACGGTTTACCGTTGcgatc                                                                                                                                          <  1:875526/149‑5 (MQ=255)
             acGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCttt                                                                                                                                         <  2:583482/150‑1 (MQ=255)
              cGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTTTACCGTTGcgatcgc                                                                                                                                        <  1:506607/150‑7 (MQ=255)
                  gCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCGCAACAGAGGTGCAATGGTTGCGGCGCGGTTCACCGCCAATATCGGTTATCGCGTAGTTGTTACAGTCAGCGATCGCAACGGTAAGCCGTTGCCCTTTGgcgc                                                                                                                                    >  2:518348/1‑150 (MQ=255)
                  gCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGgcgc                                                                                                                                    >  2:593830/1‑150 (MQ=255)
                  gCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGgcgc                                                                                                                                    <  2:1274919/150‑1 (MQ=255)
                    tGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGCTCGCAACGGTAAACCGTTGCCCTTTGCCGctc                                                                                                                                  >  2:789442/1‑150 (MQ=255)
                     gCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCACCAGAGGTGCAATGGTAGCGGCGCGTTTCAACTCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGCTCGCAACGGTAAACCGTTGCGCTTTGGCGctct                                                                                                                                 >  2:275450/1‑150 (MQ=255)
                     gCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCCGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGctct                                                                                                                                 >  2:1620329/1‑150 (MQ=255)
                      cccGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCtt                                                                                                                                >  2:1346149/1‑150 (MQ=255)
                       ccGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTg                                                                                                                               >  2:1708983/1‑150 (MQ=255)
                        cGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGc                                                                                                                              <  1:238709/150‑1 (MQ=255)
                         gATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTTTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGcc                                                                                                                             >  1:298109/1‑150 (MQ=255)
                          aTAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCTCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCCCAACGGTAAACCCTTGCCCTTTGGCGCTCTTGCCa                                                                                                                            >  1:334987/1‑150 (MQ=255)
                           tcaCGTCGATCTTGACCAAACTACACAGTTTGTGGTGCCCATCAGAGGTGCAATGGTAGCGGATCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAg                                                                                                                           <  1:669035/148‑1 (MQ=255)
                           tAAAGTCTATCTTGAACAAAAAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTATCGGCGCGTTTCAACGCCAATATCTGTTATCGCGTACTTGTTACCGGCAGCGATCGCAACGGTAAACCGTTGCGCTTTGGCGCTCTTGCCAg                                                                                                                           <  2:1209803/150‑1 (MQ=255)
                              cGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCGAACAGAGGTGCAATGGTAGCGGCGCGGTTCGACGCCAATAGCGGTTATCGCGTCCTTGTTAGAGTCAGCGATCGCAACGGTACACCGTTGCCCTTTGGCGCTCTTGGCAGccc                                                                                                                        >  2:1113305/1‑148 (MQ=255)
                              cGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCaa                                                                                                                        >  1:1314783/1‑150 (MQ=255)
                                tCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACg                                                                                                                      >  1:58532/1‑150 (MQ=255)
                                  gATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACgat                                                                                                                    <  2:1179421/150‑1 (MQ=255)
                                   aTCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACgatg                                                                                                                   <  2:414683/150‑1 (MQ=255)
                                   aTCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACgatg                                                                                                                   <  1:1318694/150‑1 (MQ=255)
                                        gAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTAACGCGTACTTGTTACAGACAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACg                                                                                                              >  1:1390108/1‑150 (MQ=255)
                                        gAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCACGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACg                                                                                                              >  2:1642886/1‑150 (MQ=255)
                                          acaaacaaCACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGATTCAACGCCAAGATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAAACGTTGCCCTTTGGCGCTCTTGCCATCAACGATGATACggg                                                                                                            >  2:912545/1‑150 (MQ=255)
                                          acaaacaaCACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGAGTCACCGCCAATATCGGTTATCGCGTACTTGTTACAGTCCGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACggg                                                                                                            >  1:978827/1‑150 (MQ=255)
                                           caaacaaCACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACgggg                                                                                                           <  1:105935/150‑1 (MQ=255)
                                             aacaacACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGca                                                                                                         <  2:157648/150‑1 (MQ=255)
                                             aacaacACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGCTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGca                                                                                                         <  2:1668701/150‑1 (MQ=255)
                                               caacaCAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGTTGATACGGGGgaac                                                                                                       >  2:126087/1‑150 (MQ=255)
                                                 acacAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACATTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGcaacaa                                                                                                     >  2:827045/1‑150 (MQ=255)
                                                  cacaGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACaaa                                                                                                    <  1:5545/150‑1 (MQ=255)
                                                   acaGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAg                                                                                                   <  2:1190850/150‑1 (MQ=255)
                                                    caGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGt                                                                                                  <  1:1346149/150‑1 (MQ=255)
                                                     aGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTa                                                                                                 <  1:925715/150‑1 (MQ=255)
                                                      gTTTGTGGTGCCCAACAGAGGAGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTAt                                                                                                <  1:282288/150‑1 (MQ=255)
                                                        ttGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCCATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTAtcg                                                                                              <  1:1088711/150‑1 (MQ=255)
                                                         tgtgGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGACCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTAtcgt                                                                                             <  2:669161/150‑1 (MQ=255)
                                                         tgtgGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTGACAGTCAGCGCTCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGCTGATACGGGGCAACAAAGTAtcgg                                                                                             >  1:1309968/1‑149 (MQ=255)
                                                          gtggtgCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTAtcgtc                                                                                            <  1:714843/150‑1 (MQ=255)
                                                          gtggtgCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTAtcgtc                                                                                            <  1:567977/150‑1 (MQ=255)
                                                                cccAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAg                                                                                      >  2:919480/1‑150 (MQ=255)
                                                                cccAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAg                                                                                      >  1:548024/1‑150 (MQ=255)
                                                                       gagGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCa                                                                               <  1:1663726/150‑1 (MQ=255)
                                                                       gagGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCa                                                                               <  1:159786/150‑1 (MQ=255)
                                                                         ggTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATa                                                                             >  2:1038961/1‑150 (MQ=255)
                                                                           tGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACt                                                                           <  2:45627/150‑1 (MQ=255)
                                                                              aaTGGTAGCGTCGCGTTTCCACGCCAATATCTGTTCTCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCACCGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACtata                                                                        <  1:310839/150‑1 (MQ=255)
                                                                                    aGCGGCGCGTTTCAACGCCAAAATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATAtctctc                                                                  <  1:1708983/150‑1 (MQ=255)
                                                                                     gcggcgCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATctctct                                                                 >  1:620981/1‑150 (MQ=255)
                                                                                         cgcgTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGggcggc                                                                                <  1:382630/131‑1 (MQ=255)
                                                                                         cgcgTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGggcggc                                                                                >  2:382630/1‑131 (MQ=255)
                                                                                            gTTTCACCGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGatat                                                          >  2:1647284/1‑150 (MQ=255)
                                                                                              ttCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCg                                                        <  2:1681929/150‑1 (MQ=255)
                                                                                              ttCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGAGGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCg                                                        <  2:497774/150‑1 (MQ=255)
                                                                                                cAACGCCAATATCGGTTATCGCGAACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTCGCGCTCTTGCCAGCAACGATCATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCAGACTATATCTCTCTGGGATATCGAg                                                      >  2:1664918/1‑150 (MQ=255)
                                                                                                    gcgattatGGGTTATCGCGTACTTGTTACAGTCAGCGTTCGCAACGGTAAACCGTTGCCCTTTGGCGATCTTGCCAGCAACGATGATACGGGGCAACAAATTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTaaa                                                  <  2:1196295/145‑1 (MQ=255)
                                                                                                     ccAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAAt                                                 <  2:1117003/150‑1 (MQ=255)
                                                                                                      cAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGAAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATc                                                <  2:794339/150‑1 (MQ=255)
                                                                                                      cAATATCGGTTATCGCGTACTTGGAACAGTCAGCGATCGCACCGGTAAACCGCGGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAATGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATc                                                <  2:165026/150‑1 (MQ=255)
                                                                                                        atatCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATcac                                              <  2:1459945/150‑1 (MQ=255)
                                                                                                           tCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACaaa                                           <  2:636846/150‑1 (MQ=255)
                                                                                                           tCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACaaa                                           <  2:488575/150‑1 (MQ=255)
                                                                                                            cGGCTATCGCGTACATGTTACACTCAGCGCTCGCAACGGTAACCCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGATGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAg                                          <  1:1113305/150‑1 (MQ=255)
                                                                                                              gTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGACGGCGGCATACTATATCTCTGAGGGATAACGAGTAAATCACAAAGCt                                        >  2:1590371/1‑150 (MQ=255)
                                                                                                                         cTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGc                             >  1:1599864/1‑150 (MQ=255)
                                                                                                                             ttACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTggg                         <  1:518348/150‑1 (MQ=255)
                                                                                                                             ttACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTggg                         <  1:593830/150‑1 (MQ=255)
                                                                                                                              tACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTgggg                        <  2:1081204/150‑1 (MQ=255)
                                                                                                                              aaCAGTCAGCGCTCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTgggg                        <  1:1715029/149‑1 (MQ=255)
                                                                                                                                cAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGaa                      >  1:1420668/1‑150 (MQ=255)
                                                                                                                                    cAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGCGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCa                  >  2:336960/1‑150 (MQ=255)
                                                                                                                                       cgctCGCAACTGTAAACCGTAGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGc               <  1:275450/147‑1 (MQ=255)
                                                                                                                                       cGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGc               <  2:999732/150‑1 (MQ=255)
                                                                                                                                       cGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGc               >  1:853576/1‑150 (MQ=255)
                                                                                                                                       cGATCGCAACGGTAAACCGTTGCCCTATGGCGCTCTTGCCAGCAACGATGATACGGGGCAACATAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGc               <  2:194925/150‑1 (MQ=255)
                                                                                                                                        aagggCAACGGTTTACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCa              >  1:674855/5‑150 (MQ=255)
                                                                                                                                          gggCAACGGTTTACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCata            >  1:692306/3‑148 (MQ=255)
                                                                                                                                          gggCAACGGTTTACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGa            >  1:962167/3‑150 (MQ=255)
                                                                                                                                           ggCAACGGTTTACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGAt           >  2:1059809/2‑150 (MQ=255)
                                                                                                                                            gCAACGGTAAACCGGTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATc          <  1:1627055/150‑1 (MQ=255)
                                                                                                                                             cAACGGTACACCGTTGCCCTTTGGCGCTCTTGCCATCAACGATGATAGGGGGCAACAATGTTTCGTCTATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCGGGACTGTACGCGGGGGAAATCAGGCAGATca         <  2:903411/150‑1 (MQ=255)
                                                                                                                                             cAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATca         <  1:1642886/150‑1 (MQ=255)
                                                                                                                                                ccgTAAACCGGTGCCCCTTGGCGCTCTAGCCAGCAACGATGCTCCGGGGCAACAAAGTATCGTAGATGATGTCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATcaaca      <  1:41143/148‑1 (MQ=255)
                                                                                                                                                    aaaCCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTTTCGTCGATGAGGGCGGCATACTATATCTATCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGt  >  2:89700/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTCGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAATGT  >  NC_000913/3290961‑3291258

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