Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,760,118 G→A 45.3% W378* (TGG→TGA yfjI → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,760,1180GA45.3% 15.5 / 100.5 82W378* (TGG→TGAyfjICP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (26/19);  new base A (21/16);  total (47/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-01

GAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAGATTACT  >  NC_000913/2760019‑2760215
                                                                                                   |                                                                                                 
gaACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTgg                                                                                                   <  2:1260287/100‑1 (MQ=255)
gaACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTgg                                                                                                   <  1:457593/100‑1 (MQ=255)
gaACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTgg                                                                                                   <  2:463945/100‑1 (MQ=255)
  aCATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGta                                                                                                 >  2:978167/1‑100 (MQ=255)
  aCATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGta                                                                                                 <  1:2297141/100‑1 (MQ=255)
  aCATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGta                                                                                                 <  2:421101/100‑1 (MQ=255)
  aCATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGAta                                                                                                 <  2:1960082/100‑1 (MQ=255)
   cATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGatat                                                                                                >  2:1518647/1‑100 (MQ=255)
   cATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGatat                                                                                                >  2:2415000/1‑100 (MQ=255)
   cATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGatat                                                                                                >  2:810097/1‑100 (MQ=255)
    aTGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGtata                                                                                               <  1:57684/100‑1 (MQ=255)
      ggCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATAtt                                                                                             >  1:2477841/1‑100 (MQ=255)
       gCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTg                                                                                            <  1:126140/100‑1 (MQ=255)
       gCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTg                                                                                            <  1:877998/100‑1 (MQ=255)
        cGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTga                                                                                           <  2:688282/100‑1 (MQ=255)
          tggCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCTCAGTAAAAGCAGCAGTGGATATAGTGGCTTGATATAGTgaag                                                                                         <  1:386110/99‑1 (MQ=21)
          aGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTgaag                                                                                         <  1:18746/100‑1 (MQ=255)
             cTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAAt                                                                                      <  2:974081/100‑1 (MQ=255)
             cTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGGGGCTTGATATATTGAAGAAt                                                                                      <  2:2437440/100‑1 (MQ=255)
              ttGCTGCCTTATTTCTTTATTTTATCGGTGATGGAGGCGATATACCTGCTACCGCAGTAAAAGCAGCAGCGGAGTTAGTGGCTTGGTATTTTGAAGAAt                                                                                      <  1:153083/99‑1 (MQ=21)
               tgctgcGTTACTTCATTACTTCAGCGGTGATGGAGGTGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAAt                                                                                      <  1:530708/98‑1 (MQ=255)
               tgctgcCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATAc                                                                                    >  2:324008/1‑100 (MQ=255)
                 ctgcCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACAt                                                                                  >  1:1235030/1‑100 (MQ=255)
                   gcCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACATcc                                                                                >  1:908324/1‑100 (MQ=255)
                     cTTACTTCATTACTTCAGCGGTGAGGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGTAGTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGc                                                                              >  1:2145764/1‑100 (MQ=255)
                      ttacttcattacttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCt                                                                             <  1:2396629/100‑1 (MQ=255)
                       tacttcattacttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCtt                                                                            <  2:304241/100‑1 (MQ=255)
                         cttcattacttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGt                                                                          >  1:295511/1‑100 (MQ=255)
                         cttcattacttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGt                                                                          >  2:2212268/1‑100 (MQ=255)
                           tcattacttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGTtc                                                                        >  2:1327525/1‑100 (MQ=255)
                             attacttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGTtctc                                                                      <  1:2540429/100‑1 (MQ=255)
                              ttacttcgGCGTTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGGTTGATATATTGAAGAATACATCCGCTTGTTctct                                                                     >  1:1357653/1‑100 (MQ=38)
                              ttacttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGTTctct                                                                     >  2:372143/1‑100 (MQ=255)
                                cctgcaGCGGAGATTGAGGCGATATATCTGTTACCGCAGTATACGCAGCCGTGCAGATAGTGGTGTGGTCTATGGAAGTATGCATCCGCTTGTTgtct                                                                     >  1:2525538/5‑98 (MQ=16)
                                acttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTaa                                                                   >  2:359641/1‑100 (MQ=255)
                                 cttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGTTCTCTaaa                                                                  >  2:2302735/1‑100 (MQ=255)
                                  ttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTaaaa                                                                 >  2:1928886/1‑100 (MQ=255)
                                  ttcaGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGTTCTCTaaaa                                                                 >  2:837139/1‑100 (MQ=255)
                                    caGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTaaaaaa                                                               >  2:1721009/1‑100 (MQ=255)
                                     aGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGTTCTCTAAAAaag                                                              <  1:2533735/100‑1 (MQ=255)
                                       cGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAaagaa                                                            >  2:839093/1‑100 (MQ=255)
                                        ggtggtGGAGGCGATATATGTGTTACCGCATTAAAAGCAGGAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAaagaag                                                           <  1:291698/100‑1 (MQ=255)
                                          tGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAagag                                                         <  2:2581374/100‑1 (MQ=255)
                                             tGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGttt                                                      >  1:1208864/1‑100 (MQ=255)
                                              ggAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGtttt                                                     >  2:709900/1‑100 (MQ=255)
                                                aGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCt                                                   >  2:1252482/1‑100 (MQ=255)
                                                 ggCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGCAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCtt                                                  <  1:1870147/100‑1 (MQ=255)
                                                 ggCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCtt                                                  >  1:821858/1‑100 (MQ=255)
                                                  gCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCttt                                                 <  2:2099757/100‑1 (MQ=255)
                                                   cGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTa                                                <  1:1982043/100‑1 (MQ=255)
                                                     atatatCTTTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGa                                              <  2:2387295/100‑1 (MQ=255)
                                                       atatCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATg                                            >  1:326978/1‑100 (MQ=255)
                                                        tatCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGt                                           >  2:2492563/1‑100 (MQ=255)
                                                         atCTGTCACCGCAGTAAAAGCAGCAGTTGAGATAGTGGATTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGtt                                          <  1:611093/100‑1 (MQ=255)
                                                          tCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGttt                                         >  1:871860/1‑100 (MQ=255)
                                                            tGTTACCGCAGTAGAAGAAGCAGGGGAAAAAGTGGCTTGGAATATTGAAGAATACATGGGTTGGTAAACTAAAAAAGAAAAGGTATCATTAAATggt                                          >  1:1268730/1‑95 (MQ=255)
                                                            tGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCa                                       <  1:2018852/100‑1 (MQ=255)
                                                            tGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCa                                       >  1:2396352/1‑100 (MQ=255)
                                                             gTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAg                                      >  1:1682248/1‑100 (MQ=255)
                                                                aCCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAg                                   >  1:2005865/1‑100 (MQ=255)
                                                                  cGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCa                                 >  1:1940175/1‑100 (MQ=255)
                                                                   gCAGTAAAAGCAGCATTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAg                                <  2:2548488/100‑1 (MQ=255)
                                                                   gCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAg                                >  2:1959526/1‑100 (MQ=255)
                                                                       tAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGa                            >  2:467705/1‑100 (MQ=255)
                                                                          aaGCAGCAGTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCt                         >  1:1801389/1‑100 (MQ=255)
                                                                           agcagcAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCtt                        >  1:1736425/1‑100 (MQ=255)
                                                                            gcagcaGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCtt                        <  1:497584/99‑1 (MQ=255)
                                                                              agcagTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTAt                     <  2:1701137/100‑1 (MQ=255)
                                                                              agcagTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTAt                     <  2:1557354/100‑1 (MQ=255)
                                                                                 agTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTAttgt                  >  1:943785/1‑100 (MQ=255)
                                                                                  gTGGAGATAGTGGCTTGGTATATTGAAGTATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTAttg                   <  1:1452903/98‑1 (MQ=255)
                                                                                  gTGGAGATAGTGGCTTGATATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTAttgtt                 <  1:2212253/100‑1 (MQ=255)
                                                                                   tGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTAttgttg                <  1:2326797/100‑1 (MQ=255)
                                                                                    ggAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTgg               <  1:246453/100‑1 (MQ=255)
                                                                                       gaTAGTGGCTTGATATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATa            >  1:2518645/1‑100 (MQ=255)
                                                                                        aTAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAGGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATaa           >  1:279812/1‑100 (MQ=255)
                                                                                         tAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATaaa          >  1:1840120/1‑100 (MQ=255)
                                                                                         tAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATaaa          >  2:950991/1‑100 (MQ=255)
                                                                                         tAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATaaa          >  1:1020059/1‑100 (MQ=255)
                                                                                         tAGTGGCTTGATATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATaaa          >  1:2310716/1‑100 (MQ=255)
                                                                                          aGTGGCTTGATATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTACTTTAGATGTTTCAGAAGCAGATGAGCTTTAATGTTGGATaaa          >  1:296478/1‑99 (MQ=255)
                                                                                           gTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAg        >  2:2344742/1‑100 (MQ=255)
                                                                                           gTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAg        <  2:8802/100‑1 (MQ=255)
                                                                                           gTGGCTTGATATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAg        >  1:107216/1‑100 (MQ=255)
                                                                                              gCTTGATATATTGAAGAATACATCCGCTTGTTGTCTAAATAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAGAtt     <  1:1819146/100‑1 (MQ=255)
                                                                                               cTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAGATTa    >  2:2471511/1‑100 (MQ=255)
                                                                                               cTTGATATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAGATTa    >  1:1907002/1‑100 (MQ=255)
                                                                                                 tGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAGATTACt  >  2:6679/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
GAACATGGCGAGGCTTGCTGCCTTACTTCATTACTTCAGCGGTGATGGAGGCGATATATCTGTTACCGCAGTAAAAGCAGCAGTGGAGATAGTGGCTTGGTATATTGAAGAATACATCCGCTTGTTCTCTAAAAAAGAAGAGTTTTCTTTAGATGTTTCAGAAGCAGATGAGCTTTATTGTTGGATAAAAGATTACT  >  NC_000913/2760019‑2760215

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