Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,672,460 G→A H83H (CAC→CAT yphB ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,672,4600GA98.8% 877.6 / ‑2.0 248H83H (CAC→CATyphBhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (140/105);  total (141/107)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

Only 100 of 262 total aligned reads displayed.

TGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGC  >  NC_000913/2672371‑2672549
                                                                                         |                                                                                         
tGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCatgt                                                                                        <  2:8515461/93‑1 (MQ=255)
   tAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGa                                                                                     >  1:5435610/1‑93 (MQ=255)
     gACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATa                                                                                   >  1:4935598/1‑93 (MQ=255)
     gACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATa                                                                                   >  1:1020849/1‑93 (MQ=255)
      acacCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAg                                                                                  <  2:1457324/93‑1 (MQ=255)
      acacCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAg                                                                                  <  1:1753606/93‑1 (MQ=255)
       cacCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGt                                                                                 <  1:2540252/93‑1 (MQ=255)
             tgggATGTTCATACCCCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCAt                                                                           <  1:7455786/90‑1 (MQ=255)
             tGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCAt                                                                           >  1:3602514/1‑93 (MQ=255)
              gCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGTCCAGCCAGCCATCGCCATGTAGATAGTGCGCATc                                                                          >  2:3458004/1‑93 (MQ=255)
               cGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATcc                                                                         >  2:3945941/1‑93 (MQ=255)
                gATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATccc                                                                        >  1:7852857/1‑93 (MQ=255)
                 aTGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCa                                                                       >  1:963960/1‑93 (MQ=255)
                 aTGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCa                                                                       >  2:7700619/1‑93 (MQ=255)
                     tCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCa                                                                   >  2:3530387/1‑93 (MQ=255)
                       aTACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAAc                                                                 >  2:5541612/1‑93 (MQ=255)
                        tACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACa                                                                >  2:2747754/1‑93 (MQ=255)
                         acacCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGACAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACAt                                                               >  1:6655752/1‑93 (MQ=255)
                         acacCAGACACAAACTATCATCGCTATGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACAt                                                               <  1:6514029/93‑1 (MQ=255)
                          ccccAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACAtt                                                              <  1:1683410/91‑1 (MQ=255)
                           acCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTc                                                             >  1:8418943/1‑93 (MQ=255)
                           acCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTc                                                             <  2:6904363/93‑1 (MQ=255)
                            ccAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCg                                                            <  1:456776/93‑1 (MQ=255)
                             cAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGCTAGTGCGCATCCCACTCAACATTCgg                                                           >  1:5301848/1‑93 (MQ=255)
                              aGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGt                                                          >  1:1051653/1‑93 (MQ=255)
                              aGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGt                                                          <  2:8697462/93‑1 (MQ=255)
                               gACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGtt                                                         <  1:72469/93‑1 (MQ=255)
                                  acaAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCa                                                      <  2:626051/93‑1 (MQ=255)
                                  acaAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCa                                                      <  1:7776990/93‑1 (MQ=255)
                                   caAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAg                                                     <  2:1437422/93‑1 (MQ=255)
                                   caAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAg                                                     <  1:7142652/93‑1 (MQ=255)
                                   caAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAg                                                     >  2:4752625/1‑93 (MQ=255)
                                    aaaCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGt                                                    <  2:3244135/93‑1 (MQ=255)
                                    aaaCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGt                                                    >  1:1396281/1‑93 (MQ=255)
                                     aaCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGtt                                                   <  1:5979277/93‑1 (MQ=255)
                                       cTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGa                                                 >  2:7898206/1‑93 (MQ=255)
                                          tcatcGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATAc                                              <  1:8545603/93‑1 (MQ=255)
                                          tcatcGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATAc                                              >  1:8047183/1‑93 (MQ=255)
                                          tcatcGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATAc                                              >  2:2646423/1‑93 (MQ=255)
                                           catcGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACt                                             >  1:7972170/1‑93 (MQ=255)
                                           catcGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACt                                             <  2:1175512/93‑1 (MQ=255)
                                            atcGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCAATCAACATTCGGTTGCAGTTGATACTc                                            >  2:3877455/1‑93 (MQ=255)
                                             tcGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCa                                           >  1:1585898/1‑93 (MQ=255)
                                             tcGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCa                                           >  2:8103125/1‑93 (MQ=255)
                                              cGCTGTGCGAGACACATTGCCTCTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCAc                                          <  2:5615462/93‑1 (MQ=255)
                                              cGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCAc                                          <  2:5329822/93‑1 (MQ=255)
                                              cGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCAc                                          <  1:388913/93‑1 (MQ=255)
                                               gCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACg                                         >  1:3203021/1‑93 (MQ=255)
                                               gCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACg                                         <  2:6414143/93‑1 (MQ=255)
                                               gCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACg                                         <  2:1363475/93‑1 (MQ=255)
                                                cTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGc                                        <  2:8164889/93‑1 (MQ=255)
                                                  gtgCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGccc                                      >  1:5897576/1‑93 (MQ=255)
                                                   tgCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGcccc                                     <  1:4255767/93‑1 (MQ=255)
                                                    gCGAGCCCCATTGCCACTCGCCCAGCGAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCt                                    <  1:3681890/93‑1 (MQ=255)
                                                    gCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCt                                    <  1:5170086/93‑1 (MQ=255)
                                                       agaCACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGcc                                 >  2:6351506/1‑93 (MQ=255)
                                                       agaCACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGcc                                 >  2:5052750/1‑93 (MQ=255)
                                                       agaCACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGcc                                 >  1:4868129/1‑93 (MQ=255)
                                                        gaCACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCa                                >  2:4344999/1‑93 (MQ=255)
                                                        gaCACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCa                                >  1:7539527/1‑93 (MQ=255)
                                                        gaCACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCa                                >  2:3629631/1‑93 (MQ=255)
                                                        gaCACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCa                                <  1:5746316/93‑1 (MQ=255)
                                                          cacaTTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGa                              <  2:719654/93‑1 (MQ=255)
                                                           acaTTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGAc                             >  1:5494217/1‑93 (MQ=255)
                                                             aTTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACaa                           <  1:3477336/93‑1 (MQ=255)
                                                              ttGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTAGGTTGCAGTTGATACTCACGCCCCTGCCAGACaaa                          >  1:7943322/1‑93 (MQ=255)
                                                                gCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAAc                        <  1:6953133/93‑1 (MQ=255)
                                                                gCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAAc                        >  2:2246241/1‑93 (MQ=255)
                                                                gCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAAc                        >  2:1036879/1‑93 (MQ=255)
                                                                gCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAAc                        >  1:3745421/1‑93 (MQ=255)
                                                                  cACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGa                      <  2:4248013/93‑1 (MQ=255)
                                                                   ccTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGAt                     <  2:3383228/92‑1 (MQ=255)
                                                                    cTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGAtt                    >  2:7734511/1‑93 (MQ=255)
                                                                     tCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTg                   >  1:823428/1‑93 (MQ=255)
                                                                      cGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACAATCGGTTGCGGTTGATACTCACGCCCCTGCCAGAAAAAACGATTGc                  >  2:936400/1‑93 (MQ=255)
                                                                        cgcagccagcCATCGCCATGTAGATAGTGCGCAGCCCACTCAACATGCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCg                <  1:5005086/91‑1 (MQ=255)
                                                                        cccAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCg                >  1:1142787/1‑93 (MQ=255)
                                                                         ccagccagCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGc               >  2:7579449/1‑93 (MQ=255)
                                                                         ccagccagCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGc               <  2:718415/93‑1 (MQ=255)
                                                                          cagccagcCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCt              >  1:4025912/1‑93 (MQ=255)
                                                                          cagccagcCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCt              >  1:872857/1‑93 (MQ=255)
                                                                          cagccagcCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCt              >  1:8351150/1‑93 (MQ=255)
                                                                           agccagccATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTc             <  2:1935196/93‑1 (MQ=255)
                                                                           agccagccATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTc             >  1:6817241/1‑93 (MQ=255)
                                                                            gccagccaTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCa            >  2:3016143/1‑93 (MQ=255)
                                                                              cagccaTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATCCTCACGCCCCTGCCAGACAAAACGATTGCCGCTCAcc          <  2:6796285/93‑1 (MQ=255)
                                                                              cagccaTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCAcc          <  2:4528864/93‑1 (MQ=255)
                                                                               agccaTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCAccc         <  2:7434973/93‑1 (MQ=255)
                                                                               agccaTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCAccc         >  1:1217917/1‑93 (MQ=255)
                                                                                gccaTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCg        >  2:3070852/1‑93 (MQ=255)
                                                                                 ccaTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCgg       <  1:2356607/93‑1 (MQ=255)
                                                                                  caTCGCTATGTAGTTAGTGCTCATCCCACTGACAATGCGGTTGCAGTGGATACACTCGCCCCTGCCAGACAAAACGATTGCCGCTCGCCCGGt      <  2:509985/93‑1 (MQ=255)
                                                                                  caTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGt      <  1:7870529/93‑1 (MQ=255)
                                                                                  caTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCCGt      >  1:4762483/1‑93 (MQ=255)
                                                                                   aTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGtt     <  1:3337164/93‑1 (MQ=255)
                                                                                    tCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGttt    >  1:2917243/1‑93 (MQ=255)
                                                                                     cGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTg   >  1:7298765/1‑93 (MQ=255)
                                                                                     cGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTg   <  2:2529567/93‑1 (MQ=255)
                                                                                      gCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGc  >  1:6229363/1‑93 (MQ=255)
                                                                                      gCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGc  >  1:1927659/1‑93 (MQ=255)
                                                                                         |                                                                                         
TGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGC  >  NC_000913/2672371‑2672549

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