Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 889,923 A→G intergenic (‑104/‑166) ybjL ← / → ybjM putative transport protein YbjL/putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913889,9230AG98.6% 270.0 / ‑3.7 72intergenic (‑104/‑166)ybjL/ybjMputative transport protein YbjL/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (25/46);  minor base T (0/1);  total (25/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CCATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAACGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAA  >  NC_000913/889788‑890055
                                                                                                                                       |                                                                                                                                    
ccATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGAt                                                                                                                                  <  2:25526/140‑1 (MQ=255)
 cATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATa                                                                                                                                 <  2:900301/140‑1 (MQ=255)
  aTTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATag                                                                                                                                >  1:682615/1‑140 (MQ=255)
   tttAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagt                                                                                                                               >  2:381821/1‑140 (MQ=255)
    ttAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagta                                                                                                                              >  2:331578/1‑140 (MQ=255)
     tAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagtaa                                                                                                                             <  2:45744/140‑1 (MQ=255)
       aCAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagtaaag                                                                                                                           <  1:1399700/140‑1 (MQ=255)
           ttCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAAc                                                                                                                       >  2:1178729/1‑140 (MQ=255)
             cGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACaa                                                                                                                     <  1:1427251/140‑1 (MQ=255)
             cGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACaa                                                                                                                     <  1:656375/140‑1 (MQ=255)
             cGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACaa                                                                                                                     <  1:727064/140‑1 (MQ=255)
                 gACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAatata                                                                                                                 <  2:1316767/140‑1 (MQ=255)
                 gACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAatata                                                                                                                 <  1:593248/140‑1 (MQ=255)
                 gACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAatata                                                                                                                 <  2:682615/140‑1 (MQ=255)
                          atTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAg                                                                                                        <  1:276538/140‑1 (MQ=255)
                             cACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAgtgt                                                                                                     >  2:533337/1‑140 (MQ=255)
                              aCGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTa                                                                                                    >  2:483977/1‑140 (MQ=255)
                              aCGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTa                                                                                                    <  1:805310/140‑1 (MQ=255)
                              aCGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTa                                                                                                    <  1:77019/140‑1 (MQ=255)
                              aCGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTa                                                                                                    >  2:926653/1‑140 (MQ=255)
                                ggAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAAt                                                                                                  >  2:266045/1‑140 (MQ=255)
                                 gAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATc                                                                                                 <  1:205287/140‑1 (MQ=255)
                                 gAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATc                                                                                                 <  1:1261069/140‑1 (MQ=255)
                                  aGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATca                                                                                                <  2:289126/140‑1 (MQ=255)
                                    gCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATcaca                                                                                              <  2:138669/140‑1 (MQ=255)
                                     cTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACAt                                                                                             >  1:978518/1‑140 (MQ=255)
                                     cTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACAt                                                                                             <  2:685237/140‑1 (MQ=255)
                                      tAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACAtt                                                                                            >  1:1307137/1‑140 (MQ=255)
                                        aCTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAg                                                                                          >  1:806384/1‑140 (MQ=255)
                                        aCTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAg                                                                                          <  2:1121268/140‑1 (MQ=255)
                                            cttGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTAt                                                                                      >  1:785609/1‑140 (MQ=255)
                                              tGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCa                                                                                    >  1:230154/1‑140 (MQ=255)
                                               gTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCaa                                                                                   <  1:236371/140‑1 (MQ=255)
                                                 ttACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACg                                                                                 <  1:460516/140‑1 (MQ=255)
                                                 ttACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACg                                                                                 <  2:806384/140‑1 (MQ=255)
                                                  tACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACgg                                                                                >  1:988972/1‑140 (MQ=255)
                                                   aCTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGc                                                                               <  1:80339/140‑1 (MQ=255)
                                                    cTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCt                                                                              >  1:831833/1‑140 (MQ=255)
                                                      aGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCtat                                                                            <  2:310528/140‑1 (MQ=255)
                                                       gTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCtata                                                                           >  2:145856/1‑140 (MQ=255)
                                                         aaGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATg                                                                         <  1:65319/140‑1 (MQ=255)
                                                                tGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCg                                                                  <  1:806486/140‑1 (MQ=255)
                                                                tGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCg                                                                  <  1:483977/140‑1 (MQ=255)
                                                                 gaaagaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAATCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGt                                                                 <  1:803606/140‑1 (MQ=255)
                                                                 gaaagaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGt                                                                 <  1:926653/140‑1 (MQ=255)
                                                                    agaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTgg                                                              <  1:1139610/140‑1 (MQ=255)
                                                                     gaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGc                                                             <  1:961324/140‑1 (MQ=255)
                                                                     gaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGc                                                             <  1:886336/140‑1 (MQ=255)
                                                                      aaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCTCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGcc                                                            <  1:8939/140‑1 (MQ=255)
                                                                      aaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGcc                                                            <  1:777656/140‑1 (MQ=255)
                                                                      aaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGcc                                                            <  1:757928/140‑1 (MQ=255)
                                                                         tGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCtat                                                         >  1:686128/1‑140 (MQ=255)
                                                                          ggTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCtata                                                        <  2:683116/140‑1 (MQ=255)
                                                                          ggTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCtata                                                        <  2:831833/140‑1 (MQ=255)
                                                                          ggTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCtata                                                        <  2:978518/140‑1 (MQ=255)
                                                                          ggTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCtata                                                        <  1:1242430/140‑1 (MQ=255)
                                                                          ggTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCtata                                                        <  1:266045/140‑1 (MQ=255)
                                                                                 tACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCAt                                                 >  2:1038522/1‑140 (MQ=255)
                                                                                     aTAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGaa                                             >  1:950184/1‑140 (MQ=255)
                                                                                     aTAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGaa                                             >  1:905662/1‑140 (MQ=255)
                                                                                          gTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGa                                        <  1:381821/140‑1 (MQ=255)
                                                                                           ttttttACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGaa                                       <  2:929275/140‑1 (MQ=255)
                                                                                                   cAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCttt                               >  2:356320/1‑140 (MQ=255)
                                                                                                    aGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTa                              <  1:145856/140‑1 (MQ=255)
                                                                                                           aGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAGCGGTTATATGAATTGCGTTGGCATATATTAGCGTGGAATGCGAAGCGGCTTTATCTTACt                       <  1:357846/140‑1 (MQ=255)
                                                                                                           aGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACt                       <  1:249677/140‑1 (MQ=255)
                                                                                                             ggCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGa                     >  2:1055784/1‑140 (MQ=255)
                                                                                                                 cAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGc                 <  2:901811/140‑1 (MQ=255)
                                                                                                                    attcattcAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCcac              >  1:803380/1‑140 (MQ=255)
                                                                                                                     ttcattcaGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCcaca             <  1:1178729/140‑1 (MQ=255)
                                                                                                                                cACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGaaaaa  >  1:65753/1‑140 (MQ=255)
                                                                                                                                cACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGaaaaa  >  2:1185259/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                    
CCATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAACGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAA  >  NC_000913/889788‑890055

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