Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,616,015 G→T V85F (GTT→TTT)  yneK → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,616,0150GT90.7% 758.3 / 52.5 269V85F (GTT→TTT) yneKuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (12/13);  new base T (133/110);  total (146/123)
Fisher's exact test for biased strand distribution p-value = 5.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 276 total aligned reads displayed.

TCAACTCAGCCACAAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGAT  >  NC_000913/1615917‑1616110
                                                                                                  |                                                                                               
   aCTCAGCCACAAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTc                                                                                              <  1:1270616/99‑1 (MQ=255)
   aCTCAGCCACAAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTc                                                                                              >  2:3440373/1‑99 (MQ=255)
    cTCAGCCACAAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCa                                                                                             <  1:130527/99‑1 (MQ=255)
    cTCAGCCACAAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCa                                                                                             <  2:1382672/99‑1 (MQ=255)
       aGCCACAAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGtt                                                                                          >  1:512492/1‑99 (MQ=255)
         ccACAAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATTCAATGAAAAGTTTCAGTTTa                                                                                        <  2:373813/99‑1 (MQ=255)
         caacaAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTa                                                                                        <  2:1823454/97‑1 (MQ=255)
          cacaAAAAAGATATTTAATAACCTTGACAACTTATCATAGATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTaa                                                                                       <  1:999745/99‑1 (MQ=255)
           ataaataaGATGTCTAATCACCTGGAGAGATTAGCATATATTGATAGTACTTCTTTTCGGGGTAACTGCACGCTATACAATCAAAAGTTTCATTTTaaa                                                                                      <  2:4268691/97‑1 (MQ=255)
           acaAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTaaa                                                                                      >  2:436513/1‑99 (MQ=255)
            caAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAAc                                                                                     >  1:618331/1‑99 (MQ=255)
             aaaaaaGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACt                                                                                    >  1:3293167/1‑99 (MQ=255)
             aaaaaaGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACt                                                                                    >  1:4684785/1‑99 (MQ=255)
                aaaGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCAc                                                                                 <  2:3824557/99‑1 (MQ=255)
                aaaGATATTTAATAACCTTGACAACTTATAATATATTCTGAGGAGTTCTTTTCGGGATAACGACACGCTATACAAGGAAAAGTTTCAGTCTAAACTCAc                                                                                 >  1:4587759/1‑99 (MQ=255)
                 aaGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACt                                                                                >  2:2664685/1‑99 (MQ=255)
                 aaGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGAGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACt                                                                                >  1:2523441/1‑99 (MQ=255)
                  aGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACtt                                                                               <  2:2146419/99‑1 (MQ=255)
                     tatTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATc                                                                            <  1:1887488/99‑1 (MQ=255)
                      atTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATcc                                                                           >  1:4984975/1‑99 (MQ=255)
                      atTTAATAACCTTGACAACTTAGCATATATTCAGATGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCATTTATcc                                                                           <  2:2451751/99‑1 (MQ=255)
                       tttAATACCCTTGACAACTTATCATATATTCATAGGACTTCTTTTCGGGGTAACGATACGCTATACGATGAAAAGTTTTAGTTTAAACTCACTTATccc                                                                          <  2:3258634/99‑1 (MQ=255)
                        ttAATAACCTTGACAACTTATCGTATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCg                                                                         >  2:3128185/1‑99 (MQ=255)
                          aataaCCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGcc                                                                       >  2:1524351/1‑99 (MQ=255)
                          aataaCCTTGACAACTTATAATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGcc                                                                       >  2:2509979/1‑99 (MQ=255)
                             aaCCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAgg                                                                    >  1:2728499/1‑99 (MQ=255)
                              aCCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGa                                                                   <  1:2059043/99‑1 (MQ=255)
                              aCCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGa                                                                   <  2:674418/99‑1 (MQ=255)
                               ccTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGCTTAAACTCACTTATCCCGCCAGGaa                                                                  >  2:4880444/1‑99 (MQ=255)
                                cTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAAt                                                                 >  2:1187086/1‑99 (MQ=255)
                                cTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAAt                                                                 >  1:3135516/1‑99 (MQ=255)
                                  tGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGTATgg                                                               <  2:902649/99‑1 (MQ=255)
                                  tGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATgg                                                               <  1:1005433/99‑1 (MQ=255)
                                    aCAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTTACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGGCAGGAATGGaa                                                             <  2:4981912/99‑1 (MQ=255)
                                    aCAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGaa                                                             <  1:4135683/99‑1 (MQ=255)
                                     cAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAg                                                            >  1:1408780/1‑99 (MQ=255)
                                      aaCTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGa                                                           <  2:1895806/99‑1 (MQ=255)
                                      aaCTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGa                                                           >  1:1103571/1‑99 (MQ=255)
                                       aCTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGac                                                          <  2:4606748/99‑1 (MQ=255)
                                       aCTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGac                                                          <  1:2640880/99‑1 (MQ=255)
                                        cTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGaca                                                         <  1:870098/99‑1 (MQ=255)
                                         ttATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGacac                                                        <  1:3926461/99‑1 (MQ=255)
                                         ttATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGacac                                                        >  1:1952528/1‑99 (MQ=255)
                                          tATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGAcaca                                                       <  1:4546710/99‑1 (MQ=255)
                                           aTCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACaa                                                      >  2:605408/1‑99 (MQ=255)
                                            tCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAGGACACaaa                                                     >  2:3995450/1‑99 (MQ=255)
                                              atatatTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGa                                                   >  2:2325386/1‑99 (MQ=255)
                                               tatatTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGaa                                                  >  2:1537971/1‑99 (MQ=255)
                                               tatatTCAGAGGACTTCTTTTAGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGaa                                                  >  2:3671879/1‑99 (MQ=255)
                                               tatatTCAGAGCACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGaa                                                  <  2:3895727/99‑1 (MQ=255)
                                                atatTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGaaa                                                 >  1:1962576/1‑99 (MQ=255)
                                                 tatTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGaaaa                                                <  1:1034647/99‑1 (MQ=255)
                                                   ttCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAata                                              >  2:5111783/1‑99 (MQ=255)
                                                   ttCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAata                                              >  1:4095818/1‑99 (MQ=255)
                                                   ttCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAata                                              <  2:1259919/99‑1 (MQ=255)
                                                    tCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAatat                                             >  1:4739407/1‑99 (MQ=255)
                                                     cAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATAtt                                            <  1:74092/99‑1 (MQ=255)
                                                      agagGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCACGCCAGGAATGGAAGACACAAAGAAAATATTg                                           >  2:280116/1‑99 (MQ=255)
                                                        agGACTCCTTTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCTCCAGGAATGGAAGACACAAAGAAAATATTGaa                                         >  2:5180457/1‑99 (MQ=255)
                                                           aCTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAAttt                                      <  1:3007764/99‑1 (MQ=255)
                                                             ttcttTTCGGGGTAACGTCATGCTATACAATGATAAGTTTAAGTTTAAACTCACTTAGCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCa                                    <  1:1143186/99‑1 (MQ=255)
                                                             ttcttTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCa                                    >  2:3150073/1‑99 (MQ=255)
                                                             ttcttTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCa                                    >  1:1563886/1‑99 (MQ=255)
                                                             ttcttTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCa                                    >  2:2109554/1‑99 (MQ=255)
                                                             ttcttTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCa                                    >  1:2957124/1‑99 (MQ=255)
                                                             ttcttTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCa                                    >  1:320656/1‑99 (MQ=255)
                                                              tcttTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAg                                   <  2:2063987/99‑1 (MQ=255)
                                                               cttTTCGGTGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAgg                                  <  1:3391308/99‑1 (MQ=255)
                                                               cttTTCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAgg                                  <  1:1761706/99‑1 (MQ=255)
                                                                 tttCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGgta                                <  1:4246071/99‑1 (MQ=255)
                                                                   tCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGgtagt                              >  1:3615112/1‑99 (MQ=255)
                                                                   tCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGgtagt                              >  2:2973384/1‑99 (MQ=255)
                                                                   tCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGgtagt                              <  2:2833539/99‑1 (MQ=255)
                                                                   tCGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGgtagt                              <  1:2815839/99‑1 (MQ=255)
                                                                    cGGGGTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGgtagta                             >  1:2317744/1‑99 (MQ=255)
                                                                     ggggTAACGACACGCTATGCAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGGATTTCAGGTAGtaa                            >  2:4715981/1‑99 (MQ=255)
                                                                      gggTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGtaat                           >  2:795293/1‑99 (MQ=255)
                                                                      gggTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGtaat                           >  2:1636033/1‑99 (MQ=255)
                                                                        gTAACGACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGtaataa                         >  2:5084664/1‑99 (MQ=255)
                                                                        gTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGtaataa                         <  2:4810848/99‑1 (MQ=255)
                                                                        gTAACGACACGCTATACAATAAAAAAATTAAGTTTAAACTCACTTATCACGACAATAATAAAAAAAACAAAGAAAATATTAAATTTCAAAtaataataa                         >  2:4327057/1‑99 (MQ=255)
                                                                         tAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATaaa                        <  2:723565/99‑1 (MQ=255)
                                                                             gACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTa                    >  1:4198384/1‑99 (MQ=255)
                                                                             gACACGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGAAATGGAAGACACAAAGAAAATATTGAATTTAAGGTAGTAATAAATTTa                    >  1:5278499/1‑99 (MQ=255)
                                                                              acacGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTaa                   >  1:2704967/1‑99 (MQ=255)
                                                                               cacGCTATACAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAg                  <  1:67919/99‑1 (MQ=255)
                                                                               cacGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAg                  >  2:4028595/1‑99 (MQ=255)
                                                                                    tataCAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCta             >  2:1097103/1‑99 (MQ=255)
                                                                                      taCAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCtatt           <  2:223738/99‑1 (MQ=255)
                                                                                      taCAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCtatt           >  1:951366/1‑99 (MQ=255)
                                                                                       aCAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCtattt          <  1:4818912/99‑1 (MQ=255)
                                                                                       aCAATGAAAAGTTTCAGTTTAAACTAACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCtattt          >  2:865577/1‑99 (MQ=255)
                                                                                        cAATGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCtagtta         <  2:4392694/99‑1 (MQ=255)
                                                                                         aaTGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCtatttat        >  1:2125605/1‑99 (MQ=255)
                                                                                         aaTGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCtatttat        >  2:2399986/1‑99 (MQ=255)
                                                                                         aaTGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTAAATTTCAGGTAGTAATAAATTTAAGCCCtatttat        >  1:3787433/1‑99 (MQ=255)
                                                                                           tGAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTatttattt      <  2:451151/99‑1 (MQ=255)
                                                                                            gAAAAGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTAtttattta     >  2:4277006/1‑99 (MQ=255)
                                                                                              aaaGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTCGa   >  2:3076671/1‑99 (MQ=255)
                                                                                              aaaGTTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGa   >  2:1452744/1‑99 (MQ=255)
                                                                                                  |                                                                                               
TCAACTCAGCCACAAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGAT  >  NC_000913/1615917‑1616110

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