Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,620,798 Δ1 bp intergenic (+255/+186) ydeF → / ← ydeH putative transport protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,620,7950C.99.4% 796.7 / ‑3.5 176intergenic (+252/+189)ydeF/ydeHputative transport protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (86/89);  total (87/89)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

Only 100 of 199 total aligned reads displayed.

CTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAT  >  NC_000913/1620703‑1620883
                                                                                            |                                                                                        
cTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGc                                                                                          <  1:1969604/93‑1 (MQ=255)
 tGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGcc                                                                                         >  2:2758020/1‑93 (MQ=255)
  gggAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGccc                                                                                        >  1:1584810/1‑93 (MQ=255)
   ggAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGccct                                                                                       >  2:1021094/1‑92 (MQ=255)
    gAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGccctt                                                                                      >  1:3009976/1‑91 (MQ=255)
     aataatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGcccttg                                                                                     >  1:3417891/1‑90 (MQ=255)
     aataatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGcccttg                                                                                     >  2:3320323/1‑90 (MQ=255)
       taatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGtt                                                                                  >  1:4360846/1‑93 (MQ=255)
       taatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGtt                                                                                  <  1:692203/93‑1 (MQ=255)
       taatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGtt                                                                                  <  2:3797658/93‑1 (MQ=255)
         attttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTAc                                                                                <  1:2959687/93‑1 (MQ=255)
          ttttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACa                                                                               >  1:1534718/1‑93 (MQ=255)
            ttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGt                                                                             >  1:4515066/1‑93 (MQ=255)
              attttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTtg                                                                           >  2:2170438/1‑93 (MQ=255)
               ttttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTtgt                                                                          >  2:2630540/1‑93 (MQ=255)
               ttttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTtgt                                                                          <  2:4261519/93‑1 (MQ=255)
               ttttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTtgt                                                                          >  1:2642088/1‑93 (MQ=255)
                 ttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTgtgt                                                                        >  1:2775820/1‑93 (MQ=255)
                 ttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTgtgt                                                                        <  1:2484138/93‑1 (MQ=255)
                  tctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTa                                                                       <  2:183047/93‑1 (MQ=255)
                  tctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTa                                                                       <  2:4244157/93‑1 (MQ=255)
                   ctggttttctggTCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTAc                                                                      <  1:2311713/93‑1 (MQ=255)
                   ctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTAc                                                                      <  1:2597705/93‑1 (MQ=255)
                    tggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACa                                                                     >  1:3409652/1‑93 (MQ=255)
                     ggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACaa                                                                    >  2:1619392/1‑93 (MQ=255)
                     ggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACaa                                                                    >  1:3539756/1‑93 (MQ=255)
                            tggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGcc                                                              >  2:4103865/1‑93 (MQ=255)
                            tggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCt                                                             >  1:1058637/1‑93 (MQ=255)
                            tggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCt                                                             >  2:4455939/1‑93 (MQ=255)
                            tggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCt                                                             <  1:2351641/93‑1 (MQ=255)
                            tggCCGCGTATTTCAGCCACAAAGGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCt                                                             <  2:3412894/93‑1 (MQ=255)
                             ggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTg                                                            >  1:2091806/1‑93 (MQ=255)
                                cgcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGAtt                                                         >  1:3551955/1‑93 (MQ=255)
                                cgcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGAtt                                                         >  1:1020922/1‑93 (MQ=255)
                                 gcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGAttt                                                        >  1:426382/1‑93 (MQ=255)
                                 gcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGAttt                                                        >  2:3204190/1‑93 (MQ=255)
                                  cgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGAtttt                                                       <  1:2145349/93‑1 (MQ=255)
                                  cgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGAtttt                                                       >  1:1313606/1‑93 (MQ=255)
                                    tATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTa                                                     <  1:24121/93‑1 (MQ=255)
                                    tATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTa                                                     <  2:3682961/93‑1 (MQ=255)
                                    tATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTa                                                     >  2:3962103/1‑93 (MQ=255)
                                    tATTTCAGCCACAAATGGGATAACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTa                                                     <  2:3479139/93‑1 (MQ=255)
                                     aTTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTAt                                                    <  1:1672758/93‑1 (MQ=255)
                                     aTTTAAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTAt                                                    <  1:2449281/93‑1 (MQ=255)
                                      tttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATg                                                   >  2:899139/1‑93 (MQ=255)
                                      tttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATg                                                   >  1:1851913/1‑93 (MQ=255)
                                      tttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATg                                                   <  1:1960260/93‑1 (MQ=255)
                                      tttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATg                                                   >  2:4669241/1‑93 (MQ=255)
                                      tttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATg                                                   <  2:869919/93‑1 (MQ=255)
                                       ttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGa                                                  >  1:1816604/1‑93 (MQ=255)
                                        tCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGAc                                                 <  1:843009/93‑1 (MQ=255)
                                        tCAGCCACAAATGAGATTACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTTATTTTTATGAc                                                 <  1:4239928/93‑1 (MQ=255)
                                           gCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGc                                              <  1:4122237/93‑1 (MQ=255)
                                             cacaAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGa                                            >  2:4160314/1‑93 (MQ=255)
                                             cacaAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGa                                            <  1:2693153/93‑1 (MQ=255)
                                                aaaTGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaa                                         >  1:1947422/1‑93 (MQ=255)
                                                  aTGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaaaa                                       >  1:4013597/1‑93 (MQ=255)
                                                  aTGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaaaa                                       >  1:3245035/1‑93 (MQ=255)
                                                   tGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaaaaa                                      <  2:2317716/93‑1 (MQ=255)
                                                    gggATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaaaaaa                                     <  2:4299423/93‑1 (MQ=255)
                                                     ggATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAAc                                    <  1:3189682/93‑1 (MQ=255)
                                                      gATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAAcc                                   <  2:1082861/93‑1 (MQ=255)
                                                        tGACTAATGAACGGAGATAATCCCTCCCCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGc                                 <  1:3506538/93‑1 (MQ=255)
                                                        tGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATGTGTATGACGGCGAAAAAAAAACGc                                 >  1:4254950/1‑93 (MQ=255)
                                                         gACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGcc                                >  1:3589140/1‑93 (MQ=255)
                                                         gACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGcc                                >  2:2464115/1‑93 (MQ=255)
                                                          aCTAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCa                               <  1:3306861/93‑1 (MQ=255)
                                                            tAATGAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGt                             <  2:1226116/93‑1 (MQ=255)
                                                             aaTGAGCGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTa                            <  1:1952771/93‑1 (MQ=255)
                                                               agAACGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTaaa                          <  2:1306869/92‑1 (MQ=255)
                                                                 aaCGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAAcc                        <  2:1656530/93‑1 (MQ=255)
                                                                 aaCGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAAcc                        >  2:2300177/1‑93 (MQ=255)
                                                                 aaCGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAAcc                        <  2:1516182/93‑1 (MQ=255)
                                                                  aCGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcg                       <  2:3666289/93‑1 (MQ=255)
                                                                  aCGGAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACAGCGAAAAAAAACCGCCAGTAAACcg                       >  2:783861/1‑93 (MQ=255)
                                                                    ggAGATAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggc                     <  1:665322/93‑1 (MQ=255)
                                                                    agagaTAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggc                     <  1:4374269/92‑1 (MQ=255)
                                                                      agaTAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggcgg                   <  2:710772/93‑1 (MQ=255)
                                                                       gaTAATCCCTCCCCTACCCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGt                  <  2:1752970/93‑1 (MQ=255)
                                                                       gaTAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggcggg                  >  1:1192213/1‑92 (MQ=255)
                                                                       gaTAATCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGt                  <  2:631344/93‑1 (MQ=255)
                                                                          aaTCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGaa               <  1:3842258/93‑1 (MQ=255)
                                                                          aaTCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGaa               <  1:2870859/93‑1 (MQ=255)
                                                                          aaTCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGaa               >  1:486291/1‑93 (MQ=255)
                                                                          aaTCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGaa               <  1:1325326/93‑1 (MQ=255)
                                                                          aaTCCCTCACCTAACCGGCCCTTGTTACAGTTGTGGACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCCGTAAACCGGCGGCGaa               >  1:4522217/1‑93 (MQ=255)
                                                                            tCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATg             >  2:4555220/1‑93 (MQ=255)
                                                                            tCCCTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATg             >  1:4694617/1‑93 (MQ=255)
                                                                            ccccTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATg             <  2:3786582/92‑1 (MQ=255)
                                                                              ccTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGCAAAAAAACCGCCAGTAAACCGGCGGTGAATGCt           >  2:2788762/1‑93 (MQ=255)
                                                                              ccTCACCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCt           <  2:3099834/93‑1 (MQ=255)
                                                                                  gccTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCa       <  1:3695447/92‑1 (MQ=255)
                                                                                  aCCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCa       >  2:4638051/1‑93 (MQ=255)
                                                                                  aCCTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCa       >  2:2986765/1‑93 (MQ=255)
                                                                                   ccTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCAt      >  1:2842308/1‑93 (MQ=255)
                                                                                    cTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg     <  2:2842045/93‑1 (MQ=255)
                                                                                    cTAACCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg     >  1:3161539/1‑93 (MQ=255)
                                                                                      aaCCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGa   >  1:1988276/1‑93 (MQ=255)
                                                                                      aaCCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGa   >  1:3502631/1‑93 (MQ=255)
                                                                                       aCCGGCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt  >  2:1002576/1‑93 (MQ=255)
                                                                                            |                                                                                        
CTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAT  >  NC_000913/1620703‑1620883

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