Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 573,627 T→C intergenic (+65/‑77) ybcQ → / → b4509 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913573,6270TC99.6% 862.0 / ‑4.8 264intergenic (+65/‑77)ybcQ/b4509hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (154/109);  total (155/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

Only 100 of 266 total aligned reads displayed.

ATATCGTTGTTAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCA  >  NC_000913/573537‑573717
                                                                                          |                                                                                          
atatCGTTGTTAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCaa                                                                                          >  1:3413817/1‑93 (MQ=255)
 tatCGTTGTTAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAat                                                                                         >  2:218094/1‑93 (MQ=255)
  atCGTTGTTAATAACTCTAATTAATAGGCTATTTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAata                                                                                        <  2:2549987/93‑1 (MQ=255)
    cGTTGTTAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAatata                                                                                      <  2:2705577/93‑1 (MQ=255)
    cGTTGTTAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAatata                                                                                      <  2:3773094/93‑1 (MQ=255)
      ttgttAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAAtatata                                                                                    >  2:6654077/1‑93 (MQ=255)
      ttgttAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAAtatata                                                                                    >  1:3320546/1‑93 (MQ=255)
      ttgttAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAAtatata                                                                                    >  2:32436/1‑93 (MQ=255)
      ttgttAATAACTCGCATTAATATGCCAATTGTTTACTAAAAATTATTAAAAAAGGCGCGTTGAGACGCCCCCAAGTATAAACGGCACtatata                                                                                    >  1:1719294/1‑93 (MQ=255)
       tgttAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATaa                                                                                   <  2:6389402/93‑1 (MQ=255)
          taataaCTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAg                                                                                >  1:1975294/1‑93 (MQ=255)
           aataaCTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGa                                                                               >  2:5468508/1‑93 (MQ=255)
             taaCTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAg                                                                             <  2:1518528/93‑1 (MQ=255)
             taaCTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAg                                                                             <  1:4965612/93‑1 (MQ=255)
              aaCTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAgg                                                                            >  2:3357974/1‑93 (MQ=255)
              aaCTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAgg                                                                            <  1:983262/93‑1 (MQ=255)
               aCTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGt                                                                           >  1:4097455/1‑93 (MQ=255)
                ctctAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGtt                                                                          <  2:838054/93‑1 (MQ=255)
                ctctAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGtt                                                                          >  1:3610999/1‑93 (MQ=255)
                ctctAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGtt                                                                          >  1:212068/1‑93 (MQ=255)
                 tctAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGttt                                                                         >  2:4909281/1‑93 (MQ=255)
                 tctAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGttt                                                                         >  1:1656384/1‑93 (MQ=255)
                   taattaatATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTtat                                                                       >  2:4024133/1‑93 (MQ=255)
                    aattaatATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTtata                                                                      >  1:5540844/1‑93 (MQ=255)
                    aattaatATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTtata                                                                      <  1:2802386/93‑1 (MQ=255)
                     attaatATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTtatat                                                                     >  1:966864/1‑93 (MQ=255)
                     attaatATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTtatat                                                                     <  2:308029/93‑1 (MQ=255)
                     agtaatATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTtatat                                                                     >  1:6618589/3‑93 (MQ=255)
                          tatGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTa                                                                >  2:5568931/1‑93 (MQ=255)
                          tatGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTa                                                                >  2:4370557/1‑93 (MQ=255)
                          tatGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTa                                                                >  1:425800/1‑93 (MQ=255)
                            tGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGa                                                              >  2:1629838/1‑93 (MQ=255)
                             gCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGaa                                                             <  2:722432/93‑1 (MQ=255)
                                 aTTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCaa                                                         >  1:5450041/1‑93 (MQ=255)
                                  ttGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAg                                                        <  2:143209/93‑1 (MQ=255)
                                  ttGTTTACAAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAg                                                        <  1:5813563/93‑1 (MQ=255)
                                   tGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAgg                                                       <  2:6116748/93‑1 (MQ=255)
                                    gTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGt                                                      <  1:5295310/93‑1 (MQ=255)
                                    gTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGt                                                      <  2:6092007/93‑1 (MQ=255)
                                    gTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGt                                                      >  1:3121146/1‑93 (MQ=255)
                                     tttACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAATGTTTATATAGTTAGAAGCAAGGtt                                                     >  1:4373215/1‑93 (MQ=255)
                                     tttACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGtt                                                     <  1:4648572/93‑1 (MQ=255)
                                         cTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGc                                                 >  2:3962314/1‑93 (MQ=255)
                                         cTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGc                                                 >  2:761527/1‑93 (MQ=255)
                                          tAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCt                                                <  1:5509303/93‑1 (MQ=255)
                                          tAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCt                                                >  2:381488/1‑93 (MQ=255)
                                           aaaaaTTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCtt                                               >  2:5459644/1‑93 (MQ=255)
                                            aaaaTTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTc                                              >  1:243389/1‑93 (MQ=255)
                                            aaaaTTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTc                                              >  2:2936793/1‑93 (MQ=255)
                                             aaaTTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCt                                             <  2:6237691/93‑1 (MQ=255)
                                              aaTTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTa                                            >  1:6351617/1‑93 (MQ=255)
                                                ttattaAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTaaa                                          >  1:1019794/1‑93 (MQ=255)
                                                ttattaAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTaaa                                          >  1:6511760/1‑93 (MQ=255)
                                                  attaAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAgg                                        >  1:4525253/1‑93 (MQ=255)
                                                  attaAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAgg                                        >  2:2277495/1‑93 (MQ=255)
                                                  attaAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAgg                                        >  2:4402304/1‑93 (MQ=255)
                                                   ttaAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGa                                       >  1:2587743/1‑93 (MQ=255)
                                                     aaaaaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAg                                     <  1:6465361/93‑1 (MQ=255)
                                                       aaaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTg                                   >  1:5927407/1‑93 (MQ=255)
                                                        aaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTgg                                  >  2:5753841/1‑93 (MQ=255)
                                                         aTGGGGCGTTGAGATGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGc                                 >  1:6262714/1‑93 (MQ=255)
                                                         aTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGc                                 >  2:1485932/1‑93 (MQ=255)
                                                         aTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGc                                 >  2:3171919/1‑93 (MQ=255)
                                                         aTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGc                                 >  1:2180817/1‑93 (MQ=255)
                                                          tGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCt                                >  2:2246419/1‑93 (MQ=255)
                                                           ggggCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCATGGTTGTGCTTCTAAAGGAAGTGGCtt                               <  1:851152/93‑1 (MQ=255)
                                                           ggggCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCtt                               >  1:3715205/1‑93 (MQ=255)
                                                           ggggCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCtt                               >  2:3551240/1‑93 (MQ=255)
                                                            gggCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTg                              <  2:4527995/93‑1 (MQ=255)
                                                            gggCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTg                              <  2:1319194/93‑1 (MQ=255)
                                                            gggCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTg                              >  1:4770123/1‑93 (MQ=255)
                                                            gggCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGCTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTg                              >  2:5860292/1‑93 (MQ=255)
                                                                gTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAggg                          <  2:4881903/93‑1 (MQ=255)
                                                                gTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAggg                          >  2:4260762/1‑93 (MQ=255)
                                                                 ttGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGa                         <  2:1907160/93‑1 (MQ=255)
                                                                   gagaCGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGc                       >  2:2583917/1‑93 (MQ=255)
                                                                     gaCGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCa                     >  1:6176090/1‑93 (MQ=255)
                                                                      aCGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCAc                    <  1:5684753/93‑1 (MQ=255)
                                                                        gCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACtt                  >  2:1259880/1‑93 (MQ=255)
                                                                         cccccAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTta                 >  1:2233812/1‑93 (MQ=255)
                                                                          ccccAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTtat                >  2:5989791/1‑93 (MQ=255)
                                                                           cccAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTtata               >  1:3398382/1‑93 (MQ=255)
                                                                           cccAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTtata               <  2:6440802/93‑1 (MQ=255)
                                                                           cccAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTtata               <  2:5141525/93‑1 (MQ=255)
                                                                            ccAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTatat              >  2:5054284/1‑93 (MQ=255)
                                                                             cAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATg             >  2:1837739/1‑93 (MQ=255)
                                                                              aaaaaTAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGt            <  2:344253/93‑1 (MQ=255)
                                                                               aaaaTAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGtt           >  2:2085889/1‑93 (MQ=255)
                                                                               aaaaTAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGtt           <  2:2754751/93‑1 (MQ=255)
                                                                               aaaaTAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGtt           >  1:12081/1‑93 (MQ=255)
                                                                                aaaTAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTg          >  2:467620/1‑93 (MQ=255)
                                                                                 aaTAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTgg         <  1:5183186/93‑1 (MQ=255)
                                                                                 aaTAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTgg         <  2:5340774/93‑1 (MQ=255)
                                                                                  aTAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTggg        <  1:4347837/93‑1 (MQ=255)
                                                                                  aTAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTggg        <  1:4145132/93‑1 (MQ=255)
                                                                                    aaaGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTggggg      >  1:3467561/1‑92 (MQ=255)
                                                                                    aaaGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGa      >  1:3068698/1‑93 (MQ=255)
                                                                                    aaaGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGa      <  2:3676791/93‑1 (MQ=255)
                                                                                     aaGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAg     <  2:248202/93‑1 (MQ=255)
                                                                                       gggCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGGGGc   >  1:3876419/1‑93 (MQ=255)
                                                                                          |                                                                                          
ATATCGTTGTTAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCA  >  NC_000913/573537‑573717

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