Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,291,109 A→T 5.2% K766* (AAA→TAA) ‡ yraJ → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,291,1090AT 5.2% 895.5 / 27.9 251K766* (AAA→TAA) ‡yraJputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (114/123);  new base T (5/8);  total (120/131)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.42e-01

Only 100 of 255 total aligned reads displayed.

TCGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAAT  >  NC_000913/3290962‑3291256
                                                                                                                                                   |                                                                                                                                                   
tcGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTaa                                                                                                                                                    >  1:348539/1‑149 (MQ=255)
 cGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTaaa                                                                                                                                                   <  1:455093/149‑1 (MQ=255)
   gTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAAcc                                                                                                                                                 <  2:1332758/149‑1 (MQ=255)
      gATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGtt                                                                                                                                              <  2:854048/149‑1 (MQ=255)
        tactacGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGc                                                                                                                                            <  1:1286555/149‑1 (MQ=255)
        tactacGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGc                                                                                                                                            <  1:428516/149‑1 (MQ=255)
          ctacGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTTTACCGTTGcga                                                                                                                                          <  3:732047/149‑3 (MQ=255)
             cGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCttt                                                                                                                                       <  1:25426/149‑1 (MQ=255)
              gCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTg                                                                                                                                      <  3:1003225/149‑1 (MQ=255)
                      ccGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCtt                                                                                                                              <  3:1068715/149‑1 (MQ=255)
                        gATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGc                                                                                                                            <  2:167523/149‑1 (MQ=255)
                        gATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGc                                                                                                                            >  2:733401/1‑149 (MQ=255)
                        gATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGc                                                                                                                            <  1:59763/149‑1 (MQ=255)
                              gtagATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCaa                                                                                                                      >  4:1128733/4‑149 (MQ=255)
                              gTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCaa                                                                                                                      >  2:230877/1‑149 (MQ=255)
                              gTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCaa                                                                                                                      <  3:738860/149‑1 (MQ=255)
                                 gATCTTGAACAAACAACCCAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACga                                                                                                                   <  2:101522/149‑1 (MQ=255)
                                 gATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACga                                                                                                                   <  1:1692099/149‑1 (MQ=255)
                                 gATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACga                                                                                                                   >  1:67752/1‑149 (MQ=255)
                                   tCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACgatg                                                                                                                 >  4:76718/1‑149 (MQ=255)
                                       gAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAAc                                                                                                                                                  >  2:1271770/1‑112 (MQ=255)
                                        aacaaacaACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACg                                                                                                            >  3:1460998/1‑149 (MQ=255)
                                         acaaacaaCACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACgg                                                                                                           >  3:811083/1‑149 (MQ=255)
                                           aaacaaCACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACgggg                                                                                                         >  4:370518/1‑149 (MQ=255)
                                            aacaacACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTCCCAGCAACGATGATACGGGGc                                                                                                        >  4:738803/1‑149 (MQ=255)
                                             acaacaCAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACgggg                                                                                                         >  4:656849/1‑147 (MQ=255)
                                             acaacaCAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGca                                                                                                       <  4:1043088/149‑1 (MQ=255)
                                               aacaCAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGcaac                                                                                                     >  4:1294195/1‑149 (MQ=255)
                                                acacAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGcaaca                                                                                                    >  3:157254/1‑149 (MQ=255)
                                                  acaGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCa                                                                                                                          <  2:99131/125‑1 (MQ=255)
                                                  acaGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCa                                                                                                                          >  1:99131/1‑125 (MQ=255)
                                                  acaGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACgg                                                                                                           <  1:1529470/140‑1 (MQ=255)
                                                  acaGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACgg                                                                                                           >  2:1529470/1‑140 (MQ=255)
                                                  acaGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACaaa                                                                                                  >  1:1721994/1‑149 (MQ=255)
                                                  acaGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACaaa                                                                                                  <  2:185497/149‑1 (MQ=255)
                                                  acaGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACaaa                                                                                                  <  2:70125/149‑1 (MQ=255)
                                                  acaGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACaaa                                                                                                  >  1:1183575/1‑149 (MQ=255)
                                                   caGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAg                                                                                                 <  3:984530/149‑1 (MQ=255)
                                                     gTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTa                                                                                               >  3:880507/1‑149 (MQ=255)
                                                      tttGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTAt                                                                                              >  3:227534/1‑149 (MQ=255)
                                                       ttGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGctct                                                                                                                               >  3:665223/1‑115 (MQ=255)
                                                         gtggtgCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTAtcgt                                                                                           <  1:101233/149‑1 (MQ=255)
                                                          tggtgCCCAACAGAGGTGCAATGGTAGCGGCGC‑‑TTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGa                                                                                        <  2:1576651/149‑1 (MQ=255)
                                                           ggtgCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTg                                                                                                                                             <  2:65648/97‑1 (MQ=255)
                                                               cccAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGa                                                                                     >  2:216718/1‑149 (MQ=255)
                                                                ccAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATATGGGGCAACGAAGTATCGTCGATGAg                                                                                    >  3:1070628/1‑149 (MQ=255)
                                                                ccAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAg                                                                                    >  3:1612844/1‑149 (MQ=255)
                                                                 cAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGa                                                                                        >  1:636153/1‑144 (MQ=255)
                                                                     agagGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTTTACCGTTGc                                                                                                                                            <  3:1682945/88‑1 (MQ=255)
                                                                          tGCAATGGTAGCGGCGCGTTTCAACGCCCATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATAc                                                                          >  2:1119068/1‑149 (MQ=255)
                                                                            cAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACta                                                                        >  1:1436324/1‑149 (MQ=255)
                                                                               tGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTatat                                                                     >  1:1146917/1‑149 (MQ=255)
                                                                                 gTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATAtct                                                                   >  2:484666/1‑149 (MQ=255)
                                                                                       gcgcGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGa                                                                                        >  2:255097/1‑122 (MQ=255)
                                                                                       gcgcGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATTAGGGCGGCATACTATATCTCTCTgg                                                             >  4:1101890/1‑149 (MQ=255)
                                                                                         gcgTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGa                                                           <  3:1294195/149‑1 (MQ=255)
                                                                                         gcgTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGa                                                           >  1:36688/1‑149 (MQ=255)
                                                                                           gTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACaaa                                                                                                  >  2:937292/1‑108 (MQ=255)
                                                                                              tCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTGGATGAGGGCGGCATACTATATCTCTCTGGGATATCg                                                      <  4:917860/149‑1 (MQ=255)
                                                                                              tCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCg                                                      >  2:909246/1‑149 (MQ=255)
                                                                                                aaCGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAg                                                    <  4:1070268/149‑1 (MQ=255)
                                                                                                   gCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTaa                                                 <  3:617744/149‑1 (MQ=255)
                                                                                                   gCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTaa                                                 >  4:1679145/1‑149 (MQ=255)
                                                                                                   gCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTaa                                                 >  2:391297/1‑149 (MQ=255)
                                                                                                     cAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAAt                                               >  1:867024/1‑149 (MQ=255)
                                                                                                       atatCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACgat                                                                                                                  <  3:1118085/80‑1 (MQ=255)
                                                                                                           cGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACaaa                                         <  1:1249756/149‑1 (MQ=255)
                                                                                                            ggTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAg                                        >  3:285743/1‑149 (MQ=255)
                                                                                                             gTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGc                                       >  4:812849/1‑149 (MQ=255)
                                                                                                              ttATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCt                                      <  3:350672/149‑1 (MQ=255)
                                                                                                                 tCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAg                                                    >  4:261201/1‑132 (MQ=255)
                                                                                                                   gcgTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACt                                 >  4:1491904/1‑149 (MQ=255)
                                                                                                                     gTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACggg                                                                                                          <  2:678435/74‑1 (MQ=255)
                                                                                                                        cTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACg                            >  3:1274863/1‑149 (MQ=255)
                                                                                                                         ttGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGc                           >  1:1094613/1‑149 (MQ=255)
                                                                                                                         ttGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGc                           >  3:259696/1‑149 (MQ=255)
                                                                                                                         ttGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGc                           >  3:404040/1‑149 (MQ=255)
                                                                                                                         ttGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGc                           >  1:471984/1‑149 (MQ=255)
                                                                                                                          tGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCt                          <  1:909246/149‑1 (MQ=255)
                                                                                                                            ttACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTgg                        >  2:63322/1‑149 (MQ=255)
                                                                                                                             tACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTggg                       <  4:500377/149‑1 (MQ=255)
                                                                                                                              aCAATCAGCGATCGCAACGGTAAACCGTTGCCCTTCGGCGCTCTTGCCAGCAACgac                                                                                                                  >  3:1491192/1‑56 (MQ=255)
                                                                                                                                aGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGaa                    >  2:382004/1‑149 (MQ=255)
                                                                                                                                  tCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGa                                   >  1:1499773/1‑132 (MQ=255)
                                                                                                                                  tCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGa                                   <  2:1499773/132‑1 (MQ=255)
                                                                                                                                  tCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAAt                  >  2:349058/1‑149 (MQ=255)
                                                                                                                                  tCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAAt                  >  1:281830/1‑149 (MQ=255)
                                                                                                                                  tCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAAt                  >  1:803917/1‑149 (MQ=255)
                                                                                                                                    aGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCa                <  4:157254/149‑1 (MQ=255)
                                                                                                                                     gCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAg               >  4:98270/1‑149 (MQ=255)
                                                                                                                                     caaagggCAACGGTTTACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTggg                       <  1:78758/135‑1 (MQ=255)
                                                                                                                                     caaagggCAACGGTTTACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTggg                       >  2:78758/7‑141 (MQ=255)
                                                                                                                                       gATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGc             >  3:1563762/1‑149 (MQ=255)
                                                                                                                                           gCAACGGTTTACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCa                <  1:754982/142‑1 (MQ=255)
                                                                                                                                           gCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGAt         >  4:1454947/1‑149 (MQ=255)
                                                                                                                                               cGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCATGCAGATcaac     >  2:302290/1‑149 (MQ=255)
                                                                                                                                                 gTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATcaacaa   >  4:1498546/1‑149 (MQ=255)
                                                                                                                                                 gTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATcaacaa   <  1:928511/149‑1 (MQ=255)
                                                                                                                                                  tAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTgg                        <  4:1556627/127‑1 (MQ=255)
                                                                                                                                                  tAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAAt  <  1:325053/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAAT  >  NC_000913/3290962‑3291256

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