Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,422,664 A→G intergenic (‑295/‑32) yjfY ← / → rpsF hypothetical protein/30S ribosomal subunit protein S6

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,422,6640AG87.5% 756.9 / 90.3 303intergenic (‑295/‑32)yjfY/rpsFhypothetical protein/30S ribosomal subunit protein S6
Reads supporting (aligned to +/- strand):  ref base A (15/23);  new base G (117/148);  total (132/171)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

Only 100 of 319 total aligned reads displayed.

AGTTTAGCGAGGCATACAAGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACAGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCAGAGCGAACAGGTTCCGGGCA  >  NC_000913/4422572‑4422756
                                                                                            |                                                                                            
  tttAGCGAGGCATACAAGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACggg                                                                                            >  2:2903043/1‑93 (MQ=255)
  tttAGCGAGGCATACAAGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACggg                                                                                            <  1:370769/93‑1 (MQ=255)
   ttAGCGAGGCATACAAGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACggga                                                                                           <  1:6824271/93‑1 (MQ=255)
      gCGAGGCATACAAGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGc                                                                                        >  2:7023190/1‑93 (MQ=255)
      gCGAGGCATACAAGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGc                                                                                        <  2:475502/93‑1 (MQ=255)
       cGAGGCATACAAGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCt                                                                                       <  2:5849310/93‑1 (MQ=255)
       cGAGGCATACAAGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCt                                                                                       >  1:3190124/1‑93 (MQ=255)
        gAGGCATACAAGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTg                                                                                      >  2:117490/1‑93 (MQ=255)
           gCATACAAGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGaat                                                                                   >  1:6013940/1‑93 (MQ=255)
            cATACAAGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGaata                                                                                  >  2:2834239/1‑93 (MQ=255)
            cATACAAGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACAGGAGGCTGaata                                                                                  >  2:5741075/1‑93 (MQ=255)
              tACAAGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGaataat                                                                                >  1:1417017/1‑93 (MQ=255)
              tACAAGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGaataat                                                                                >  2:1203809/1‑93 (MQ=255)
               aCAAGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATc                                                                               <  1:2585969/93‑1 (MQ=255)
                 aaGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCg                                                                             <  2:3497017/93‑1 (MQ=255)
                  aGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGt                                                                            <  1:6997751/93‑1 (MQ=255)
                  aGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGt                                                                            <  1:932445/93‑1 (MQ=255)
                   gTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTa                                                                           >  2:3341125/1‑93 (MQ=255)
                   gTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTa                                                                           >  2:1621768/1‑93 (MQ=255)
                   gTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTa                                                                           <  1:1937057/93‑1 (MQ=255)
                   gTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTa                                                                           <  2:1340401/93‑1 (MQ=255)
                   gTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTa                                                                           <  1:4932348/93‑1 (MQ=255)
                   gTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTa                                                                           <  2:3084568/93‑1 (MQ=255)
                     aCTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAg                                                                         <  1:1576559/93‑1 (MQ=255)
                       tataACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGa                                                                       <  2:1577410/93‑1 (MQ=255)
                         taACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGc                                                                     >  2:5914031/1‑93 (MQ=255)
                         taACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACAGGAGGCTGAATAATCCGTAAGGAGc                                                                     <  2:6494265/93‑1 (MQ=255)
                           aCGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCaa                                                                   <  1:1241522/93‑1 (MQ=255)
                           aCGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCaa                                                                   >  1:4375864/1‑93 (MQ=255)
                           aCGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACAGGAGGCTGAATAATCCGTAAGGAGCaa                                                                   <  1:6117014/93‑1 (MQ=255)
                           aCGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACAGGAGGCTGAATAATCCGTAAGGAGCaa                                                                   <  2:2295162/93‑1 (MQ=255)
                              cgTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTc                                                                <  2:7406161/93‑1 (MQ=255)
                                tCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGa                                                              <  2:3630160/93‑1 (MQ=255)
                                tCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGa                                                              >  1:3299626/1‑93 (MQ=255)
                                 cATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGAt                                                             >  2:3589569/1‑93 (MQ=255)
                                  aTTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATg                                                            >  2:6368605/1‑93 (MQ=255)
                                  aTTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATg                                                            <  2:470582/93‑1 (MQ=255)
                                   tttttCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGc                                                           >  2:1053701/1‑93 (MQ=255)
                                   tttttCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGc                                                           >  2:2438076/1‑93 (MQ=255)
                                    ttttCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCg                                                          <  1:1325269/93‑1 (MQ=255)
                                       tCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCa                                                       >  2:1416995/1‑93 (MQ=255)
                                         aGCCGTTCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACAGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCAtt                                                     <  2:4635849/93‑1 (MQ=255)
                                          gCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTa                                                    >  1:7391421/1‑93 (MQ=255)
                                          gCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTa                                                    >  2:3392025/1‑93 (MQ=255)
                                          gCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTa                                                    >  2:683321/1‑93 (MQ=255)
                                            cGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACg                                                  >  2:6574866/1‑93 (MQ=255)
                                               ccTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGaaa                                               >  2:5596825/1‑93 (MQ=255)
                                               ccTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGaaa                                               <  2:1079168/93‑1 (MQ=255)
                                                cTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAAt                                              <  1:4788006/93‑1 (MQ=255)
                                                 tttAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATc                                             >  1:6185690/1‑93 (MQ=255)
                                                 tttAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATc                                             <  1:7029927/93‑1 (MQ=255)
                                                 tttAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATc                                             >  1:2507050/1‑93 (MQ=255)
                                                  ttAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCg                                            <  1:5215357/93‑1 (MQ=255)
                                                  ttAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACAGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCg                                            <  1:2888910/93‑1 (MQ=255)
                                                   tAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGt                                           <  1:5327544/93‑1 (MQ=255)
                                                    aaCACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGtt                                          <  2:4075252/93‑1 (MQ=255)
                                                    aaCACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGtt                                          <  2:7455781/93‑1 (MQ=255)
                                                    aaCACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGtt                                          <  1:6484798/93‑1 (MQ=255)
                                                    aaCACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGtt                                          <  1:4493513/93‑1 (MQ=255)
                                                    aaCACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACAGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGtt                                          >  1:6287882/1‑93 (MQ=255)
                                                     acacGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGttt                                         >  2:5509732/1‑93 (MQ=255)
                                                     acacGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGttt                                         <  2:5171898/93‑1 (MQ=255)
                                                      cacGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGtttt                                        >  1:3519341/1‑93 (MQ=255)
                                                       acGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGttttt                                       <  1:1751803/93‑1 (MQ=255)
                                                        cGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACAGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTa                                      <  1:2079183/93‑1 (MQ=255)
                                                          ttCCTTGCCTCCCCGGGATTTGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTAGATGCGTCATTACGAAATCGTTTTTATg                                    >  1:3954594/1‑93 (MQ=255)
                                                          ttCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATg                                    <  1:4235563/93‑1 (MQ=255)
                                                          ttCCTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATg                                    >  2:3842172/1‑93 (MQ=255)
                                                            ccTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGt                                  >  2:7345515/1‑93 (MQ=255)
                                                            ccTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGt                                  <  2:5642892/93‑1 (MQ=255)
                                                             cTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTc                                 <  1:1005942/93‑1 (MQ=255)
                                                             cTTGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTc                                 <  1:563751/93‑1 (MQ=255)
                                                              ttGCCTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTcc                                >  1:1672869/1‑93 (MQ=255)
                                                              ttGCCTCCCCGGGATTCGGCTGACCCAGACAGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTcc                                <  2:7173128/93‑1 (MQ=255)
                                                                 ccTCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATc                             <  1:5985979/93‑1 (MQ=255)
                                                                   tCCCCGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCt                           <  1:5603526/93‑1 (MQ=255)
                                                                     cccGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGa                         >  2:5280294/1‑93 (MQ=255)
                                                                     cccGGGATTCGGCTGACCCAGACAGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGa                         <  1:7088351/93‑1 (MQ=255)
                                                                      ccGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGAc                        >  2:6208411/1‑93 (MQ=255)
                                                                      ccGGGATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGAc                        >  1:2700622/1‑93 (MQ=255)
                                                                        gggATTCGGCTGACCCAGACAGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTATGAAATCGTTTTTATGGTCCATACTGACCa                      >  1:6532382/1‑93 (MQ=255)
                                                                        gggATTCGGCTGACCCAGACAGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCa                      >  2:4244538/1‑93 (MQ=255)
                                                                         ggATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCag                     >  2:6871885/1‑93 (MQ=255)
                                                                          gATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCaga                    >  2:5933404/1‑93 (MQ=255)
                                                                          gATTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCaga                    >  1:4667102/1‑93 (MQ=255)
                                                                           aTTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCagag                   <  2:1800163/93‑1 (MQ=255)
                                                                           aTTCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCagag                   <  2:6671989/93‑1 (MQ=255)
                                                                            ttCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCAGAGc                  >  1:5887433/1‑93 (MQ=255)
                                                                            ttCGGCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCATCATTACGAAATCGTTTTTATGGTCCATCCTGACCAGAGc                  <  2:4525566/93‑1 (MQ=255)
                                                                                gCTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCAGAGCGAAc              >  1:5111668/1‑93 (MQ=255)
                                                                                 cTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCAGAGCGAACa             <  1:4538030/93‑1 (MQ=255)
                                                                                 cTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCAGAGCGAACa             <  1:1451798/93‑1 (MQ=255)
                                                                                 cTGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCAGAGCGAACa             <  2:4035614/93‑1 (MQ=255)
                                                                                  tGACCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCAGAGCGAACAg            >  2:414502/1‑93 (MQ=255)
                                                                                    aCCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCAGAGCGAACAGGt          <  1:6335772/93‑1 (MQ=255)
                                                                                    aCCCAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCAGAGCGAACAGGt          <  1:2250837/93‑1 (MQ=255)
                                                                                    aCCCAGACAGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCAGAGCGAACAGGt          >  1:4957013/1‑93 (MQ=255)
                                                                                    aCCCAGACAGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCAGAGCGAACAGGt          >  1:1147755/1‑93 (MQ=255)
                                                                                     cccAGACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCAGAGCGAACAGGtt         <  2:2098874/93‑1 (MQ=255)
                                                                                         gACGGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCAGAGCGAACAGGTTCCgg     <  2:2763610/93‑1 (MQ=255)
                                                                                            |                                                                                            
AGTTTAGCGAGGCATACAAGTACTATAACGCGTCATTTTTCAGCCGACCTTTAACACGTTCCTTGCCTCCCCGGGATTCGGCTGACCCAGACAGGAGGCTGAATAATCCGTAAGGAGCAATTCGATGCGTCATTACGAAATCGTTTTTATGGTCCATCCTGACCAGAGCGAACAGGTTCCGGGCA  >  NC_000913/4422572‑4422756

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