Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,866,718 T→C A74A (GCA→GCG)  yidR ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,866,7180TC83.5% 743.1 / 139.3 339A74A (GCA→GCG) yidRhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/22);  new base C (155/128);  total (189/150)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 349 total aligned reads displayed.

AAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGTGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATGGTCTCGCCGGTAAACGACGCGCCAGAAG  >  NC_000913/3866626‑3866809
                                                                                            |                                                                                           
 aGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGtgtg                                                                                            <  2:2305331/93‑1 (MQ=255)
  gTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTgcgc                                                                                           <  1:6290898/93‑1 (MQ=255)
   tCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGcgcg                                                                                          >  2:4196348/1‑93 (MQ=255)
     gTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGTGCGcc                                                                                        >  1:6674503/1‑93 (MQ=255)
      tAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGccc                                                                                       >  2:5782478/1‑93 (MQ=255)
      tAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGccc                                                                                       <  1:7512474/93‑1 (MQ=255)
       aaTGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCt                                                                                      >  1:988525/1‑93 (MQ=255)
        aTGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCtg                                                                                     >  1:5254603/1‑93 (MQ=255)
        aTGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCtg                                                                                     <  2:6763954/93‑1 (MQ=255)
        aTGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCtg                                                                                     >  1:144062/1‑93 (MQ=255)
         tGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCtgt                                                                                    <  2:5345716/93‑1 (MQ=255)
          gCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCtgtg                                                                                   <  2:6531184/93‑1 (MQ=255)
           ccATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGTGCGCCCTGTGa                                                                                  >  1:1225901/1‑93 (MQ=255)
           ccATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGa                                                                                  <  2:3198243/93‑1 (MQ=255)
            cATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGAc                                                                                 >  1:6088578/1‑93 (MQ=255)
            cATGTTGCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGAc                                                                                 <  1:2647158/93‑1 (MQ=255)
              tGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGAcgc                                                                               >  2:3322239/1‑93 (MQ=255)
                tttCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGTGCGCCCTGTGAcgcgc                                                                             <  2:646745/93‑1 (MQ=255)
                tttCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGAcgcgc                                                                             >  2:2221642/1‑93 (MQ=255)
                    atcaGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCgata                                                                         >  1:2962486/1‑93 (MQ=255)
                    atcaGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCgata                                                                         >  1:3694846/1‑93 (MQ=255)
                    atcaGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCgata                                                                         >  2:3661409/1‑93 (MQ=255)
                     tcaGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCgatag                                                                        >  2:5077658/1‑93 (MQ=255)
                      caGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCgataga                                                                       <  2:858019/93‑1 (MQ=255)
                       aGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCgatagat                                                                      >  2:2067670/1‑93 (MQ=255)
                       aGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCgatagat                                                                      <  2:4656846/93‑1 (MQ=255)
                       aGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCgatagat                                                                      >  1:1600052/1‑93 (MQ=255)
                        ggATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCgatagata                                                                     <  2:3511013/93‑1 (MQ=255)
                        ggATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCgacagata                                                                     >  1:767116/1‑93 (MQ=255)
                         gATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATaa                                                                    <  2:3868841/93‑1 (MQ=255)
                         gATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATaa                                                                    >  2:2157143/1‑93 (MQ=255)
                           ttttCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAacc                                                                  <  1:2832463/93‑1 (MQ=255)
                            tttCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGTGCGCCCTGTGACGCGCGATAGATAacct                                                                 >  2:5203177/1‑93 (MQ=255)
                            tttCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAacct                                                                 <  1:6747955/93‑1 (MQ=255)
                               ccGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAacctcga                                                              <  1:2171053/93‑1 (MQ=255)
                               ccGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAacctcga                                                              >  2:2278334/1‑93 (MQ=255)
                                cGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGTGCGCCCTGTGACGCGCGATAGATAacctcgac                                                             >  1:7401205/1‑93 (MQ=255)
                                     cGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGTGCGCCCAGTGACGCGCGATAGATAACCTCGACCTCGc                                                        >  2:7434888/1‑93 (MQ=255)
                                      gTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGTGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGcc                                                       <  2:7305947/93‑1 (MQ=255)
                                      gTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGcc                                                       >  2:185038/1‑93 (MQ=255)
                                       tGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGTGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGcnn                                                      >  2:2829441/1‑91 (MQ=255)
                                        ggATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGTGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCgg                                                     >  2:3079685/1‑93 (MQ=255)
                                        ggATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCgg                                                     >  1:5765868/1‑93 (MQ=255)
                                         gATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGTGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGt                                                    <  1:4087083/93‑1 (MQ=255)
                                         gATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGt                                                    >  2:3265468/1‑93 (MQ=255)
                                          aTGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTa                                                   <  1:1944292/93‑1 (MQ=255)
                                          aTGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTa                                                   >  1:6627167/1‑93 (MQ=255)
                                           tGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGTGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTAt                                                  >  1:1807995/1‑93 (MQ=255)
                                           tGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTAt                                                  >  2:5997356/1‑93 (MQ=255)
                                           tGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTAt                                                  >  2:5913755/1‑93 (MQ=255)
                                            gAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATg                                                 >  2:4905862/1‑93 (MQ=255)
                                              aaaCATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGa                                               <  2:6069578/93‑1 (MQ=255)
                                               aaCATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGat                                              <  2:2190543/93‑1 (MQ=255)
                                               aaCATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGat                                              >  1:5040173/1‑93 (MQ=255)
                                                 cATATTTCTCTGACTTCGGATGAACGGTCGCCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGatat                                            <  2:2599849/93‑1 (MQ=255)
                                                 cATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGatat                                            >  1:3205209/1‑93 (MQ=255)
                                                 cATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGatat                                            >  2:5650662/1‑93 (MQ=255)
                                                 cATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGatat                                            >  2:1494119/1‑93 (MQ=255)
                                                 cATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGatat                                            >  2:796201/1‑93 (MQ=255)
                                                  atatTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATAtt                                           <  2:6694126/93‑1 (MQ=255)
                                                  atatTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATAtt                                           >  1:4024397/1‑93 (MQ=255)
                                                    atTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGTGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTca                                         >  1:2523450/1‑93 (MQ=255)
                                                    atTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTca                                         >  2:7275750/1‑93 (MQ=255)
                                                     tttCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTcac                                        >  2:3844812/1‑93 (MQ=255)
                                                      ttCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTcaca                                       <  1:2301352/93‑1 (MQ=255)
                                                       tctcTGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCacac                                      <  1:101959/93‑1 (MQ=255)
                                                        nnctGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACg                                     <  1:1352594/91‑1 (MQ=255)
                                                        ctctGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACg                                     >  2:7006236/1‑93 (MQ=255)
                                                         tctGACTTCGGATGAACGGTCACCACGCCGACGTGTGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGc                                    >  1:5395628/1‑93 (MQ=255)
                                                         tctGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGc                                    >  1:865878/1‑93 (MQ=255)
                                                         tctGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGc                                    <  1:4838284/93‑1 (MQ=255)
                                                         tctGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGc                                    >  2:4271502/1‑93 (MQ=255)
                                                         tctGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGc                                    >  2:1299450/1‑93 (MQ=255)
                                                          ctGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCt                                   <  2:249759/93‑1 (MQ=255)
                                                           tGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTc                                  >  1:892282/1‑93 (MQ=255)
                                                           tGACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTc                                  >  1:6910735/1‑93 (MQ=255)
                                                            gACTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCa                                 >  1:563197/1‑93 (MQ=255)
                                                              cTTCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAAt                               <  2:1728205/93‑1 (MQ=255)
                                                               ttCGGATGAACGGTCACCACGCCGACGTGTGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATg                              >  2:2697265/1‑93 (MQ=255)
                                                               ttCGGATGAACGGTCACCACGCCGACGTGTGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATg                              >  1:327281/1‑93 (MQ=255)
                                                               ttCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATg                              >  2:1643934/1‑93 (MQ=255)
                                                               ttCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATg                              >  1:584740/1‑93 (MQ=255)
                                                               ttCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATg                              >  1:6567664/1‑93 (MQ=255)
                                                               ttCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATg                              <  2:6296355/93‑1 (MQ=255)
                                                               ttCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATg                              >  1:813905/1‑93 (MQ=255)
                                                                tCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATgg                             <  2:4800554/93‑1 (MQ=255)
                                                                tCGGATGAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATgg                             <  1:6810638/93‑1 (MQ=255)
                                                                      gAACGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATGGTCTCGc                       <  1:4715722/93‑1 (MQ=255)
                                                                         cGGTCACCACGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATGGTCTCGCCgg                    <  1:4993747/93‑1 (MQ=255)
                                                                            tCACCACGCCGACGTGTGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATGGTCTCGCCGGTaa                 <  2:1212873/93‑1 (MQ=255)
                                                                             caccacGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATGGTCTCGCCGGTaaa                <  1:645591/93‑1 (MQ=255)
                                                                             caccacGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATGGTCTCGCCGGTaaa                <  1:4355181/93‑1 (MQ=255)
                                                                              accacGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATGGTCTCGCCGGTAAac               >  2:5566108/1‑93 (MQ=255)
                                                                              accacGCCGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATGGTCTCGCCGGTAAac               >  1:3624351/1‑93 (MQ=255)
                                                                                 acGCCGACGTGTGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATGGTCTCGCCGGTAAacgac            >  1:6044088/1‑93 (MQ=255)
                                                                                    ccGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATGGTCTCGCCGGTAAACGAcgcg         >  1:4212026/1‑93 (MQ=255)
                                                                                     cGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATGGTCTCGCCGGTAAACGACgcgc        <  1:5691126/93‑1 (MQ=255)
                                                                                     cGACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATGGTCTCGCCGGTAAACGACgcgc        >  1:2380804/1‑93 (MQ=255)
                                                                                      gACGTGCGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATGGTCTCGCCGGTAAACGACGCGcc       <  1:3790340/93‑1 (MQ=255)
                                                                                          tgtgCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATGGTCTCGCCGGTAAACGACGCGCCAGaa   <  2:4046983/93‑1 (MQ=255)
                                                                                            |                                                                                           
AAGTCGTAATGCCATGTTTCATCAGGATTTTCCGGGCCGTGGATGAAAACATATTTCTCTGACTTCGGATGAACGGTCACCACGCCGACGTGTGCGCCCTGTGACGCGCGATAGATAACCTCGACCTCGCCGGTATGGATATTCACACGCTCAATGGTCTCGCCGGTAAACGACGCGCCAGAAG  >  NC_000913/3866626‑3866809

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