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.6% 815.6 / ‑3.9 252intergenic (‑66/‑94)yceF/yceDhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (117/134);  minor base . (1/0);  total (118/134)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 257 total aligned reads displayed.

GATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCG  >  NC_000913/1145833‑1146012
                                                                                          |                                                                                         
 atatAGCTTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAc                                                                                        <  1:3792390/93‑1 (MQ=255)
  tataGCGTATTGATGATAGGCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAcc                                                                                       >  1:6475819/1‑93 (MQ=255)
  tataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAcc                                                                                       >  2:2492050/1‑93 (MQ=255)
  tataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAcc                                                                                       <  2:1767908/93‑1 (MQ=255)
  tataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAcc                                                                                       >  2:5794842/1‑93 (MQ=255)
  tataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAcc                                                                                       >  1:3633967/1‑93 (MQ=255)
   ataGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAAccc                                                                                      >  1:2721806/1‑93 (MQ=255)
     aGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGACAATAAAGCAGCAACCCGc                                                                                    >  1:6600428/1‑93 (MQ=255)
     aGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGc                                                                                    >  1:1512407/1‑93 (MQ=255)
     aGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGc                                                                                    <  2:1882293/93‑1 (MQ=255)
     aGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGc                                                                                    >  1:3923085/1‑93 (MQ=255)
      gCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCa                                                                                   <  2:62313/93‑1 (MQ=255)
       cGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCATGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCaa                                                                                  <  2:2992390/93‑1 (MQ=255)
       cGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCaa                                                                                  <  1:5017191/93‑1 (MQ=255)
        gTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGGAACCCGCAAt                                                                                 >  1:209231/1‑93 (MQ=255)
        gTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAt                                                                                 <  1:3464695/93‑1 (MQ=255)
          aTTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgt                                                                               <  1:4227766/93‑1 (MQ=255)
            tgatgaTAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgtgt                                                                             <  1:2932417/93‑1 (MQ=255)
            tgatgaTAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgtgt                                                                             <  2:23931/93‑1 (MQ=255)
            tgatgaTAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgtgt                                                                             <  2:3991381/93‑1 (MQ=255)
             gatgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgtgtg                                                                            <  2:3756365/93‑1 (MQ=255)
             gatgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAAtgtgtg                                                                            <  2:2729247/93‑1 (MQ=255)
              atgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGc                                                                           >  2:1573562/1‑93 (MQ=255)
               tgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCg                                                                          >  2:2185936/1‑93 (MQ=255)
               tgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCg                                                                          >  2:2807901/1‑93 (MQ=255)
               tgatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCg                                                                          <  1:91514/93‑1 (MQ=255)
                gatAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGa                                                                         >  1:1077012/1‑93 (MQ=255)
                  tAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAat                                                                       >  1:5376726/1‑93 (MQ=255)
                  tAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAat                                                                       <  1:813423/93‑1 (MQ=255)
                   aGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAatt                                                                      <  2:1239487/93‑1 (MQ=255)
                   aGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAatt                                                                      >  1:3581964/1‑93 (MQ=255)
                        ttttAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTgg                                                                 >  1:4500509/1‑93 (MQ=255)
                           tAAACTATGCGCTTCGTTTTGTAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCaa                                                              >  2:5891925/1‑93 (MQ=255)
                            aaaCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCaaa                                                             <  1:377095/93‑1 (MQ=255)
                             aaCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCaaaa                                                            <  1:5975993/93‑1 (MQ=255)
                              aCTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAATCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAg                                                           <  2:4754265/93‑1 (MQ=255)
                               cTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAgg                                                          >  2:411673/1‑93 (MQ=255)
                                tATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGc                                                         >  1:2781567/1‑93 (MQ=255)
                                tATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGc                                                         <  2:4250348/93‑1 (MQ=255)
                                 aTGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCa                                                        >  1:5334567/1‑93 (MQ=255)
                                  tGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCaa                                                       <  2:3527305/93‑1 (MQ=255)
                                    cgcTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAAcc                                                     <  1:4722335/93‑1 (MQ=255)
                                    cgcTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAAcc                                                     >  1:2440433/1‑93 (MQ=255)
                                     gcTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACca                                                    >  2:3221403/1‑93 (MQ=255)
                                     gcTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACca                                                    >  2:6430316/1‑93 (MQ=255)
                                     gcTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACca                                                    >  1:5771419/1‑93 (MQ=255)
                                     gcTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACca                                                    >  2:4259692/1‑93 (MQ=255)
                                      cTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACcac                                                   >  1:47700/1‑93 (MQ=255)
                                        tCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAg                                                 <  1:585595/93‑1 (MQ=255)
                                         cGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAGTGTGTGCGAATTATTGGCAAAAGGCAACCACAgg                                                >  1:6565296/1‑93 (MQ=255)
                                          gTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGc                                               <  2:2910602/93‑1 (MQ=255)
                                          gTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGc                                               <  1:3440897/93‑1 (MQ=255)
                                                  ggTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGACTTATTGGCAAAAGGCAACCACAGGCTGCCtttt                                       >  2:5480413/1‑93 (MQ=255)
                                                  ggTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCtttt                                       <  2:1054268/93‑1 (MQ=255)
                                                    ttGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTc                                     <  1:1041200/93‑1 (MQ=255)
                                                      gATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCtt                                   <  1:996444/93‑1 (MQ=255)
                                                      gATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCtt                                   <  2:1876627/93‑1 (MQ=255)
                                                       aTGTTTGTTATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCttt                                  <  2:6723006/93‑1 (MQ=255)
                                                        tgtttgttATCAGCACTGAACGAAAATAAAGCATCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTg                                 <  1:5605954/93‑1 (MQ=255)
                                                        tgtttgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTg                                 >  1:761023/1‑93 (MQ=255)
                                                        tgtttgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTg                                 >  1:2228342/1‑93 (MQ=255)
                                                        tgtttgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTg                                 >  2:5323302/1‑93 (MQ=255)
                                                         gtttgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGa                                <  1:1081754/93‑1 (MQ=255)
                                                         gtttgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGa                                <  2:3430958/93‑1 (MQ=255)
                                                          tttgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGAc                               <  2:4682948/93‑1 (MQ=255)
                                                           ttgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGAct                              <  2:4855905/93‑1 (MQ=255)
                                                            tgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGActc                             <  1:3015205/93‑1 (MQ=255)
                                                            tgttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGActc                             <  1:1001224/93‑1 (MQ=255)
                                                             gttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTTGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGActct                            >  1:1224359/1‑93 (MQ=255)
                                                             gttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGActct                            <  1:1272234/93‑1 (MQ=255)
                                                             gttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGActct                            >  2:6287939/1‑93 (MQ=255)
                                                              ttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTa                           >  2:1176493/1‑93 (MQ=255)
                                                              ttATCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTa                           >  1:1727395/1‑93 (MQ=255)
                                                                aTCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATg                         <  1:2035931/93‑1 (MQ=255)
                                                                aTCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATg                         <  1:4993685/93‑1 (MQ=255)
                                                                 tCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGa                        >  2:651556/1‑93 (MQ=255)
                                                                 tCAGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGa                        <  2:4788054/93‑1 (MQ=255)
                                                                   aGCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACg                      <  2:5617589/93‑1 (MQ=255)
                                                                    gCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGt                     >  2:705510/1‑93 (MQ=255)
                                                                    gCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGt                     <  2:5083611/93‑1 (MQ=255)
                                                                    gCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGt                     <  2:5124132/93‑1 (MQ=255)
                                                                    gCACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCATTTTTCTTTGACTCTATGACGt                     >  2:724894/1‑93 (MQ=255)
                                                                     cACTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGtt                    <  1:1865870/93‑1 (MQ=255)
                                                                      aCTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTa                   <  1:385894/93‑1 (MQ=255)
                                                                       cTGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTAc                  <  2:6058613/93‑1 (MQ=255)
                                                                        tGAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACa                 >  1:6385669/1‑93 (MQ=255)
                                                                         gAACGAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACaa                >  1:4570959/1‑93 (MQ=255)
                                                                             gAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGtt            <  1:5484788/93‑1 (MQ=255)
                                                                             gAAAATAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGtt            <  1:1365062/93‑1 (MQ=255)
                                                                               aaaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTaa          <  2:5671722/93‑1 (MQ=255)
                                                                                aaTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAat         <  1:3492171/93‑1 (MQ=255)
                                                                                 aTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAata        >  1:11613/1‑93 (MQ=255)
                                                                                 aTAAAGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAata        <  1:5210039/93‑1 (MQ=255)
                                                                                   aaaGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATg      <  1:3023367/93‑1 (MQ=255)
                                                                                   aaaGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATg      >  2:3606911/1‑93 (MQ=255)
                                                                                    aaGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgc     >  2:374533/1‑93 (MQ=255)
                                                                                    aaGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgc     >  2:5941369/1‑93 (MQ=255)
                                                                                    aaGCAGCAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgc     <  2:1814946/93‑1 (MQ=255)
                                                                                      gcagcaACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgcgc   >  2:5013475/1‑93 (MQ=255)
                                                                                      gcagcaACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATgcgc   <  2:3491883/93‑1 (MQ=255)
                                                                                          |                                                                                         
GATATAGCGTATTGATGATAGCCATTTTAAACTATGCGCTTCGTTTTGCAGGTTGATGTTTGTTATCAGCACTGAACGAAAATAAAGCAGTAACCCGCAATGTGTGCGAATTATTGGCAAAAGGCAACCACAGGCTGCCTTTTTCTTTGACTCTATGACGTTACAAAGTTAATATGCGCG  >  NC_000913/1145833‑1146012

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