Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,477,281 Δ1 bp intergenic (‑391/‑29) b4256 ← / → yjgN hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,477,2790T.99.4% 754.2 / ‑2.9 160intergenic (‑389/‑31)b4256/yjgNhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (98/61);  total (99/61)
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.99e-01

Only 100 of 168 total aligned reads displayed.

AGTTTTTGCCAGCCTATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATAGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTT  >  NC_000913/4477187‑4477368
                                                                                            |                                                                                         
aGTTTTTGCCAGCCTATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAAt                                                                                           >  2:382346/1‑93 (MQ=255)
aGTTTTTGCCAGCCTATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAAt                                                                                           <  1:3785487/93‑1 (MQ=255)
  tttttGCCAGCCTATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAAttc                                                                                         <  2:3794397/93‑2 (MQ=255)
     ttGCCAGCCTATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAAttcagg                                                                                      >  2:4205468/1‑89 (MQ=255)
     ttGCCAGCCTATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAAttcagg                                                                                      <  1:4841701/93‑5 (MQ=255)
      tGCCAGCCTATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCTATGGATTTTTTAATTCAggg                                                                                    >  1:6062069/1‑93 (MQ=255)
      tGCCAGCCTATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAggg                                                                                    <  2:6630304/93‑1 (MQ=255)
      tGCCAGCCTATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAggg                                                                                    <  1:6542955/93‑1 (MQ=255)
       gCCAGCCTATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGa                                                                                   <  1:5859826/93‑1 (MQ=255)
         cAGCCTATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAAt                                                                                 >  2:4561053/1‑93 (MQ=255)
         cAGCCTATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAAt                                                                                 >  2:1013429/1‑93 (MQ=255)
          aGCCTATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAAtt                                                                                <  1:2087591/93‑1 (MQ=255)
          aGCCTATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAAtt                                                                                >  2:5363127/1‑93 (MQ=255)
           gCCTATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAAttt                                                                               <  2:3986190/93‑1 (MQ=255)
             cTATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAAttttt                                                                             >  2:3688228/1‑93 (MQ=255)
              tATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTa                                                                            >  1:1467606/1‑93 (MQ=255)
              tATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTa                                                                            >  1:1824459/1‑93 (MQ=255)
              tATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTa                                                                            >  1:6248426/1‑93 (MQ=255)
              tATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTa                                                                            <  2:4484515/93‑1 (MQ=255)
              tATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTa                                                                            >  1:3323872/1‑93 (MQ=255)
               aTTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTAt                                                                           >  1:5414711/1‑93 (MQ=255)
               aTTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTAt                                                                           <  1:3630249/93‑1 (MQ=255)
                ttACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATg                                                                          >  2:6452453/1‑93 (MQ=255)
                ttACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATg                                                                          >  2:3015749/1‑93 (MQ=255)
                ttACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATg                                                                          >  2:4205348/1‑93 (MQ=255)
                ttACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATg                                                                          >  2:5581432/1‑93 (MQ=255)
                   cTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGCATTTTTATTGCt                                                                       >  1:5735050/1‑93 (MQ=255)
                   cTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCt                                                                       >  1:3969592/1‑93 (MQ=255)
                    tGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTc                                                                      >  1:6328655/1‑93 (MQ=255)
                       tAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAg                                                                   >  2:4347370/1‑93 (MQ=255)
                       tAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATAGATTTTTTAATTCAGGGAATTTTTATGGCTCAAg                                                                   <  2:6328655/93‑1 (MQ=255)
                         agagTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGtt                                                                 <  1:4334637/93‑1 (MQ=255)
                         agagTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGtt                                                                 >  1:4420501/1‑93 (MQ=255)
                         agagTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGtt                                                                 >  2:2254575/1‑93 (MQ=255)
                         agagTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGtt                                                                 <  2:190809/93‑1 (MQ=255)
                          gagTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGtta                                                                <  2:1829370/93‑1 (MQ=255)
                           agTATTGTAACTTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGttat                                                               >  2:1589826/1‑93 (MQ=255)
                           agTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGttat                                                               <  2:3499059/93‑1 (MQ=255)
                            gTATTGTACCGTGGCAGAAAATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGttatt                                                              >  1:6441140/1‑93 (MQ=255)
                               ttGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAAt                                                           >  2:5438301/1‑93 (MQ=255)
                                tGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATg                                                          <  2:6077310/93‑1 (MQ=255)
                                 gTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGa                                                         <  2:5089006/93‑1 (MQ=255)
                                  tAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGaa                                                        >  2:3405135/1‑93 (MQ=255)
                                   aaCGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGaaa                                                       <  1:5583785/93‑1 (MQ=255)
                                   aaCGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGaaa                                                       <  1:3432273/93‑1 (MQ=255)
                                    aCGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAAt                                                      <  2:3968407/93‑1 (MQ=255)
                                    aCGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAAt                                                      <  1:1781496/93‑1 (MQ=255)
                                    aCGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAAt                                                      <  1:1031929/93‑1 (MQ=255)
                                     cGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATg                                                     <  1:2252182/93‑1 (MQ=255)
                                      gTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATgg                                                    >  1:6366295/1‑93 (MQ=255)
                                      gTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATgg                                                    <  1:2217295/93‑1 (MQ=255)
                                      atGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATgg                                                    <  1:5235966/92‑1 (MQ=255)
                                        ggCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAt                                                  <  1:2479792/93‑1 (MQ=255)
                                        ggCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAt                                                  <  1:2980539/93‑1 (MQ=255)
                                            gAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAtgttc                                              >  1:5985372/1‑89 (MQ=255)
                                            gAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAtgttc                                              >  2:2063249/1‑89 (MQ=255)
                                            gAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAtgttc                                              >  2:2825772/1‑89 (MQ=255)
                                            gAACATTTTACTGTCTGTATATATAGGTATATCAATGGGTTTTTTTATTTAGGGGATTTTTATTGCTCCAGGTATTAATTAAATTGGTGGTTc                                             >  2:4944013/1‑93 (MQ=255)
                                             aaCATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTcc                                            >  2:1125964/1‑93 (MQ=255)
                                              aCATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCg                                           <  1:5715606/93‑1 (MQ=255)
                                                 ttttACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTcc                                        >  2:3640745/1‑93 (MQ=255)
                                                 ttttACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTcc                                        >  1:2583956/1‑93 (MQ=255)
                                                   ttACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCa                                      >  2:3765943/1‑93 (MQ=255)
                                                     tctgtctgtATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCAtt                                    <  2:4673232/92‑1 (MQ=255)
                                                     aCTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCAtt                                    >  2:3044050/1‑93 (MQ=255)
                                                     aCTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCAtt                                    <  2:6080665/93‑1 (MQ=255)
                                                     aCTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCAtt                                    <  1:2688880/93‑1 (MQ=255)
                                                      ctgtctgtATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTc                                   <  2:5894396/93‑1 (MQ=255)
                                                         tctgtATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCgtt                                <  2:2017514/93‑1 (MQ=255)
                                                            gtATATATAAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCgtttgt                             <  1:303503/93‑1 (MQ=255)
                                                               tatataAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGtttt                          >  2:5193504/1‑93 (MQ=255)
                                                                 tataAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTc                        >  2:999243/1‑93 (MQ=255)
                                                                 tataAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTc                        >  1:3238471/1‑93 (MQ=255)
                                                                 tataAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTc                        >  2:2171081/1‑93 (MQ=255)
                                                                  ataAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCa                       >  1:1677752/1‑93 (MQ=255)
                                                                   taAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCAt                      <  1:942515/93‑1 (MQ=255)
                                                                   taAGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCAt                      >  1:4432203/1‑93 (MQ=255)
                                                                     aGTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATgg                    <  2:5944165/93‑1 (MQ=255)
                                                                      gTTATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGt                   >  2:346957/1‑93 (MQ=255)
                                                                       ttATATCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTa                  <  1:24688/93‑1 (MQ=255)
                                                                         atatCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACa                <  2:2134673/93‑1 (MQ=255)
                                                                         atatCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACa                >  1:74592/1‑93 (MQ=255)
                                                                          tatCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAg               >  2:3493358/1‑93 (MQ=255)
                                                                          tatCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAg               <  1:3323019/93‑1 (MQ=255)
                                                                          tatCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAg               >  2:4994054/1‑93 (MQ=255)
                                                                          tatCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAg               <  2:3320473/93‑1 (MQ=255)
                                                                            tCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGt             >  1:5998634/1‑93 (MQ=255)
                                                                            tCAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGt             <  2:6450411/93‑1 (MQ=255)
                                                                             cAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTg            >  2:3190409/1‑93 (MQ=255)
                                                                             cAATGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTg            >  1:2672867/1‑93 (MQ=255)
                                                                              aaTGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTga           >  2:5340101/1‑93 (MQ=255)
                                                                              aaTGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTga           >  1:3960336/1‑93 (MQ=255)
                                                                              aaTGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTga           >  2:3516671/1‑93 (MQ=255)
                                                                              aaTGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTga           >  1:1135269/1‑93 (MQ=255)
                                                                               aTGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTgag          >  1:3777355/1‑93 (MQ=255)
                                                                               aTGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTgag          >  1:3080413/1‑93 (MQ=255)
                                                                                tGGATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTgaga         >  2:4968/1‑93 (MQ=255)
                                                                                 ggATTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTgagag        >  1:868530/1‑93 (MQ=255)
                                                                                   aTTTTTTAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGat      <  1:3856788/93‑1 (MQ=255)
                                                                                      ttttAATTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATAtt   >  2:5117597/1‑93 (MQ=255)
                                                                                            |                                                                                         
AGTTTTTGCCAGCCTATTACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATAGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTT  >  NC_000913/4477187‑4477368

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