Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,145,923 T→C intergenic (‑66/‑94) yceF ← / → yceD hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,145,9230TC99.0% 618.8 / ‑5.1 201intergenic (‑66/‑94)yceF/yceDhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (98/100);  total (100/101)
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.71e-01

Only 100 of 202 total aligned reads displayed.

GATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCG  >  NC_000913/1145833‑1146012
                                                                                          |                                                                                         
gATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCaa                                                                                         >  2:166523/1‑93 (MQ=255)
 atatAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAc                                                                                        >  1:4171140/1‑93 (MQ=255)
  tataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCATCACTGAACGAAAATAAAGCAGCAAcc                                                                                       <  1:3853265/93‑1 (MQ=255)
  tataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAcc                                                                                       >  2:7566630/1‑93 (MQ=255)
  tataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAcc                                                                                       >  2:3834885/1‑93 (MQ=255)
   ataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAccc                                                                                      <  2:3791386/93‑1 (MQ=255)
      gCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCa                                                                                   <  1:6652048/93‑1 (MQ=255)
      gCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCa                                                                                   >  2:3894337/1‑93 (MQ=255)
      gCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCa                                                                                   >  1:318491/1‑93 (MQ=255)
        gTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGAAGTAACCCGAAAt                                                                                 >  2:5593073/1‑93 (MQ=255)
         tATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtg                                                                                <  2:3213832/93‑1 (MQ=255)
         tATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtg                                                                                >  2:477471/1‑93 (MQ=255)
          aTTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgt                                                                               <  2:581138/93‑1 (MQ=255)
          aTTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgt                                                                               <  2:5797307/93‑1 (MQ=255)
           ttGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgtg                                                                              >  2:7898050/1‑93 (MQ=255)
           ttGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgtg                                                                              >  2:6540924/1‑93 (MQ=255)
            tgatgaTAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgtgt                                                                             <  2:539042/93‑1 (MQ=255)
            tgatgaTAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgtgt                                                                             <  2:6236314/93‑1 (MQ=255)
             gatgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgtgtg                                                                            <  2:4045223/93‑1 (MQ=255)
              atgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGc                                                                           <  2:495841/93‑1 (MQ=255)
              atgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGc                                                                           <  1:679410/93‑1 (MQ=255)
               tgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCg                                                                          >  2:1917388/1‑93 (MQ=255)
                gatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGAAACCCGCAATGTGTGCGa                                                                         >  1:945111/1‑93 (MQ=255)
                  tAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAat                                                                       >  1:2613701/1‑93 (MQ=255)
                      cATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAattatt                                                                   <  2:5765733/93‑1 (MQ=255)
                        ttttAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTgg                                                                 >  2:5087960/1‑93 (MQ=255)
                        ttttAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTgg                                                                 >  2:3133569/1‑93 (MQ=255)
                        ttttAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTgg                                                                 >  1:5919726/1‑93 (MQ=255)
                         tttAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGc                                                                >  2:5003429/1‑93 (MQ=255)
                            aaaCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCaaa                                                             >  1:1155380/1‑93 (MQ=255)
                             aaCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCaaaa                                                            >  2:4111917/1‑93 (MQ=255)
                             aaCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCaaaa                                                            <  1:7803471/93‑1 (MQ=255)
                             aaCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCaaaa                                                            <  1:2050587/93‑1 (MQ=255)
                             aaCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCaaaa                                                            <  1:1744263/93‑1 (MQ=255)
                              aCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAg                                                           <  1:5778101/93‑1 (MQ=255)
                              aCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAg                                                           <  2:556827/93‑1 (MQ=255)
                               cTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTAAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAgg                                                          >  1:1701461/1‑93 (MQ=255)
                                tATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGc                                                         <  1:4748841/93‑1 (MQ=255)
                                tATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGc                                                         <  1:1370839/93‑1 (MQ=255)
                                  tGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCaa                                                       <  1:7244224/93‑1 (MQ=255)
                                  tGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCaa                                                       <  2:2523211/93‑1 (MQ=255)
                                    cgcTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAAcc                                                     >  1:5567061/1‑93 (MQ=255)
                                     gcTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACca                                                    <  1:6012267/93‑1 (MQ=255)
                                     gcTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACca                                                    >  2:2468512/1‑93 (MQ=255)
                                     gcTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACca                                                    >  2:1489180/1‑93 (MQ=255)
                                       ttCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACcaca                                                  >  1:516382/1‑93 (MQ=255)
                                       ttCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACcaca                                                  <  1:5463964/93‑1 (MQ=255)
                                       ttCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACcaca                                                  <  1:5545761/93‑1 (MQ=255)
                                        tCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAg                                                 <  1:6553191/93‑1 (MQ=255)
                                         cGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAgg                                                <  2:4611421/93‑1 (MQ=255)
                                          gTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGc                                               >  2:6638042/1‑93 (MQ=255)
                                           ttttGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCt                                              <  1:6364408/93‑1 (MQ=255)
                                               gCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCt                                          >  1:4000251/1‑93 (MQ=255)
                                                 aGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATCATTGGCAAAAGGCAACCACAGGCTGCCttt                                        >  2:2917476/1‑93 (MQ=255)
                                                  ggTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCtttt                                       <  1:3303461/93‑1 (MQ=255)
                                                   gTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCttttt                                      >  1:5930763/1‑93 (MQ=255)
                                                   gTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCttttt                                      <  1:5479484/93‑1 (MQ=255)
                                                   gTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCttttt                                      >  2:7769758/1‑93 (MQ=255)
                                                    ttGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCtttttt                                     >  2:2170817/1‑92 (MQ=255)
                                                    ttGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTc                                     >  2:1290843/1‑93 (MQ=255)
                                                         gtttgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGa                                >  2:4068043/1‑93 (MQ=255)
                                                            tgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGActc                             >  1:2241452/1‑93 (MQ=255)
                                                            tgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGActc                             <  1:4589915/93‑1 (MQ=255)
                                                             gttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGActct                            >  1:5765020/1‑93 (MQ=255)
                                                              ttATCAGCACTGAACGAAAATAAAGCAGCCACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCAGCCTTTTTCTTTGACTCTa                           <  2:7240844/93‑1 (MQ=255)
                                                              ttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTa                           <  2:609377/93‑1 (MQ=255)
                                                              ttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTa                           >  2:2129397/1‑93 (MQ=255)
                                                              ttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTa                           <  1:5438174/93‑1 (MQ=255)
                                                              ttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTa                           >  1:1930421/1‑93 (MQ=255)
                                                               tATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTAt                          >  2:392594/1‑93 (MQ=255)
                                                                aTCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATg                         <  1:6885038/93‑1 (MQ=255)
                                                                  cAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGAc                       <  2:7461469/93‑1 (MQ=255)
                                                                   aGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACg                      >  2:4251426/1‑93 (MQ=255)
                                                                    gCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGt                     >  2:6780021/1‑93 (MQ=255)
                                                                     cACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGtt                    <  1:2500893/93‑1 (MQ=255)
                                                                      aCTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTa                   >  1:7385415/1‑93 (MQ=255)
                                                                       cTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTAc                  >  1:562486/1‑93 (MQ=255)
                                                                       cTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTAc                  <  1:5204785/93‑1 (MQ=255)
                                                                        tGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACa                 >  2:3383612/1‑93 (MQ=255)
                                                                        tGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACa                 >  2:4527817/1‑93 (MQ=255)
                                                                         gAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACaa                >  2:2026082/1‑93 (MQ=255)
                                                                             gAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGtt            >  1:3004655/1‑93 (MQ=255)
                                                                             gAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGtt            >  2:5217153/1‑93 (MQ=255)
                                                                             gAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGtt            >  2:6138842/1‑93 (MQ=255)
                                                                              aaaaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTa           <  1:5270857/93‑1 (MQ=255)
                                                                              aaaaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTa           <  1:7051209/93‑1 (MQ=255)
                                                                               aaaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTaa          <  2:1823767/93‑1 (MQ=255)
                                                                               aaaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTaa          >  2:302929/1‑93 (MQ=255)
                                                                                 aTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAata        <  1:7566631/93‑1 (MQ=255)
                                                                                 aTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAata        <  1:3783719/93‑1 (MQ=255)
                                                                                 aTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAata        <  1:5175970/93‑1 (MQ=255)
                                                                                 aTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAata        <  1:6820403/93‑1 (MQ=255)
                                                                                 aTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCGGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAata        >  2:3614898/1‑93 (MQ=255)
                                                                                  tAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAatat       <  2:4344266/93‑1 (MQ=255)
                                                                                  tAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAatat       <  2:2728499/93‑1 (MQ=255)
                                                                                   aaaGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATg      <  1:2798571/93‑1 (MQ=255)
                                                                                    aaGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgc     <  2:3650080/93‑1 (MQ=255)
                                                                                      gcagcaACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgcgt   <  1:5773738/93‑2 (MQ=255)
                                                                                       cAGCAACCCGCTATGTGTGCGAATTATTGGCAAAAGGCCACCACCTGCTGCCTTTTTCTTTGACTCTATGAAGTTACCAAGTTAATATgcgtg  >  1:1138444/1‑91 (MQ=255)
                                                                                       cAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgcgtg  >  1:5043138/1‑91 (MQ=255)
                                                                                          |                                                                                         
GATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCG  >  NC_000913/1145833‑1146012

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