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,6270TC97.9% 447.6 / ‑1.0 144intergenic (+65/‑77)ybcQ/b4509hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (96/45);  total (97/47)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

Only 100 of 145 total aligned reads displayed.

TATCGTTGTTAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGC  >  NC_000913/573538‑573716
                                                                                         |                                                                                         
tatCGTTGTTAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAat                                                                                        >  1:2109373/1‑93 (MQ=255)
tatCGTTGTTAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAat                                                                                        <  1:1896953/93‑1 (MQ=255)
 atCGTTGTTAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCAAAAATAAAGGGCAatat                                                                                      >  2:1534167/1‑93 (MQ=255)
 atCGTTGTTAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATAAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAata                                                                                       <  1:3358472/93‑1 (MQ=255)
   cGTTGTTAACAACTCTAATTAATATGACAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAatata                                                                                     <  1:1098719/93‑1 (MQ=255)
    gttgttAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAatatat                                                                                    >  1:3601408/1‑93 (MQ=255)
    gttgttAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAatatat                                                                                    >  1:846699/1‑93 (MQ=255)
    gttgttAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAatatat                                                                                    <  2:442247/93‑1 (MQ=255)
    gttgttAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAatatat                                                                                    <  2:3238486/93‑1 (MQ=255)
     ttgttAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAAtatata                                                                                   <  2:2369476/93‑1 (MQ=255)
       gttAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAAc                                                                                 >  2:1677538/1‑93 (MQ=255)
       gttAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAAc                                                                                 <  2:237947/93‑1 (MQ=255)
       gttAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAAc                                                                                 >  2:1677539/1‑93 (MQ=255)
          aataaCTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGa                                                                              <  2:3905482/93‑1 (MQ=255)
            taaCTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAg                                                                            >  1:4676149/1‑93 (MQ=255)
              aCTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGt                                                                          >  2:1801568/1‑93 (MQ=255)
              aCTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGt                                                                          >  2:2847217/1‑93 (MQ=255)
              aCTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGt                                                                          >  2:4419497/1‑93 (MQ=255)
               ctctAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGtt                                                                         >  2:3648656/1‑93 (MQ=255)
                tctAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGttt                                                                        >  2:4393036/1‑93 (MQ=255)
                 ctAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTta                                                                       >  1:4583390/1‑93 (MQ=255)
                  taattaatATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTtat                                                                      <  2:171395/93‑1 (MQ=255)
                  taattaatATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTtat                                                                      <  2:171396/93‑1 (MQ=255)
                  taattaatATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTtat                                                                      <  1:2734057/93‑1 (MQ=255)
                  taattaatATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTtat                                                                      <  2:3131676/93‑1 (MQ=255)
                  taattaatATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTtat                                                                      >  2:4175161/1‑93 (MQ=255)
                  taattaatATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGAAGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTtat                                                                      >  2:3228602/1‑93 (MQ=255)
                   aattaatATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTtata                                                                     >  2:506714/1‑93 (MQ=255)
                    attaatATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTtatat                                                                    >  1:3054585/1‑93 (MQ=255)
                      taatATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAg                                                                  <  2:809074/93‑1 (MQ=255)
                       aatATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGt                                                                 >  1:4548774/1‑93 (MQ=255)
                        atatGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGtt                                                                >  1:674428/1‑93 (MQ=255)
                        atatGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGtt                                                                <  1:3113525/93‑1 (MQ=255)
                         tatGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTa                                                               >  2:1204911/1‑93 (MQ=255)
                          atGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAg                                                              <  2:1990726/93‑1 (MQ=255)
                          atGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAg                                                              >  1:1643176/1‑93 (MQ=255)
                          atGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAg                                                              >  2:1362806/1‑93 (MQ=255)
                           tGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCACAAATAAAGGCCAATGTATAACAGAAGGTTGATATAGTTAGa                                                             >  1:1520722/1‑93 (MQ=255)
                           tGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGa                                                             >  1:1730719/1‑93 (MQ=255)
                           tGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGa                                                             >  2:1213648/1‑93 (MQ=255)
                            gCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGaa                                                            <  2:3156801/93‑1 (MQ=255)
                              cAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGc                                                          >  1:996655/1‑93 (MQ=255)
                              cAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGc                                                          >  2:636000/1‑93 (MQ=255)
                                  tGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAgg                                                      >  1:2751980/1‑93 (MQ=255)
                                    tttACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGtt                                                    >  2:1387120/1‑93 (MQ=255)
                                     ttACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTtg                                                   >  1:1012429/1‑93 (MQ=255)
                                      tACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTtgt                                                  <  1:1839594/93‑1 (MQ=255)
                                       aCTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTtgtg                                                 >  1:1061820/1‑93 (MQ=255)
                                           aaaaTTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTc                                             >  1:1408947/1‑93 (MQ=255)
                                            aaaTTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCt                                            >  2:3972099/1‑93 (MQ=255)
                                                tattaAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAg                                        >  1:2773837/1‑93 (MQ=255)
                                                tattaAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAg                                        <  1:4250835/93‑1 (MQ=255)
                                                 attaAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAgg                                       >  2:4186131/1‑93 (MQ=255)
                                                 attaAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAgg                                       <  2:277087/93‑1 (MQ=255)
                                                   taAAAATGGGGCGTTGAGACGCCACCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGaa                                     <  1:2204399/93‑1 (MQ=255)
                                                     aaaaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGt                                   <  2:2152828/93‑1 (MQ=255)
                                                     aaaaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGt                                   >  1:2108595/1‑93 (MQ=255)
                                                     aaaaTGGGGCGTTGAGACGCCCCAAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCGAGGTTGTGCTTCTAAAGGAAGt                                   <  2:2706221/93‑1 (MQ=255)
                                                      aaaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTg                                  >  1:3177422/1‑93 (MQ=255)
                                                      aaaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTg                                  >  1:4141916/1‑93 (MQ=255)
                                                       aaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGCAGCAAGGTTGTGCTTCTAAAGGAAGTgg                                 >  2:278763/1‑93 (MQ=255)
                                                       aaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTgg                                 >  1:3698359/1‑93 (MQ=255)
                                                       aaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTgg                                 >  2:1883747/1‑93 (MQ=255)
                                                       aaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTgg                                 >  1:632382/1‑93 (MQ=255)
                                                       aaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTgc                                 >  2:4029075/1‑92 (MQ=255)
                                                        aTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGc                                >  2:2829792/1‑93 (MQ=255)
                                                         tGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGCTGTGCTTCTAAAGGAAGTGGCt                               >  2:383534/1‑93 (MQ=255)
                                                          ggggCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCtt                              >  2:1477083/1‑93 (MQ=255)
                                                          ggggCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCtt                              >  2:3067842/1‑93 (MQ=255)
                                                          ggggCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCtt                              >  1:2409308/1‑93 (MQ=255)
                                                            ggCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGa                            >  2:2753459/1‑93 (MQ=255)
                                                               gTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAggg                         <  1:4705005/93‑1 (MQ=255)
                                                               gTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAggg                         >  1:1829218/1‑93 (MQ=255)
                                                                ttGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGa                        <  2:3718897/93‑1 (MQ=255)
                                                                    gaCGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGccc                    >  2:3472514/1‑92 (MQ=255)
                                                                       gCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACtt                 >  2:4063611/1‑93 (MQ=255)
                                                                       gCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACtt                 >  1:2672200/1‑93 (MQ=255)
                                                                         ccccAAAAATAAAGGGCAATATATAACAGAAGGTTTGTATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTtat               <  1:255994/93‑1 (MQ=255)
                                                                         ccccAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTtat               >  2:517295/1‑93 (MQ=255)
                                                                         ccccAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTtat               >  1:4667861/1‑93 (MQ=255)
                                                                          cccAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTtata              >  2:2474223/1‑93 (MQ=255)
                                                                            cAAAAATAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATg            >  1:3142213/1‑93 (MQ=255)
                                                                               aaaTAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTg         >  1:2612839/1‑93 (MQ=255)
                                                                                aaTAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTgg        <  2:2170007/93‑1 (MQ=255)
                                                                                 aTAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTggg       <  1:1368334/93‑1 (MQ=255)
                                                                                 aTAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTggg       >  2:4303563/1‑93 (MQ=255)
                                                                                 aTAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTggg       <  2:1667574/93‑1 (MQ=255)
                                                                                 aTAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTggg       <  1:2859188/93‑1 (MQ=255)
                                                                                  tAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGTAGCCACTTATATGTTgggg      >  1:521655/1‑93 (MQ=255)
                                                                                  tAAAGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTgggg      >  2:2501130/1‑93 (MQ=255)
                                                                                   aaaGGGCAATATATAACAGAAGGTTTATATAGTTCGAAGCAAGGTTGTGCTTCTAACGGAAGTGGCTTGAGGGAGCCACTTATATGTTGAGGa     >  1:1925372/1‑93 (MQ=255)
                                                                                    aaGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGNCTTGNGGGAGCCACTTATATGTTGGGGAg    >  2:2684088/1‑93 (MQ=255)
                                                                                    aaGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAg    <  1:4750109/93‑1 (MQ=255)
                                                                                    aaGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAg    <  1:1229849/93‑1 (MQ=255)
                                                                                    aaGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAg    <  1:1229708/93‑1 (MQ=255)
                                                                                     aGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAgg   <  2:2974664/93‑1 (MQ=255)
                                                                                     aGGGCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGATGTGCTTCTAAAGGAAGTGGCTTGAGGGGGCCACTTATATGTTGgggggg   >  1:3759818/1‑93 (MQ=255)
                                                                                      gggCAATATATAACAGAAGGTTTTTATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGc  >  2:3511393/1‑93 (MQ=255)
                                                                                      gggCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGGGGc  >  1:391813/1‑93 (MQ=255)
                                                                                      gggCAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGc  >  2:3613210/1‑93 (MQ=255)
                                                                                         |                                                                                         
TATCGTTGTTAATAACTCTAATTAATATGCCAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGC  >  NC_000913/573538‑573716

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: