Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,862,688 A→G 5.0% intergenic (‑259/‑83) msrB ← / → gapA methionine sulfoxide reductase B/glyceraldehyde‑3‑phosphate dehydrogenase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,862,6880AG 5.0% 313.5 / 5.2 99intergenic (‑259/‑83)msrB/gapAmethionine sulfoxide reductase B/glyceraldehyde‑3‑phosphate dehydrogenase A
Reads supporting (aligned to +/- strand):  ref base A (49/45);  new base G (2/3);  total (51/48)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

Only 100 of 105 total aligned reads displayed.

CCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGTTTTGGCCGTATCGGTCGCATTGTTTTCCGTGCT  >  NC_000913/1862541‑1862830
                                                                                                                                                   |                                                                                                                                              
ccTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGa                                                                                                                                                >  1:206354/1‑148 (MQ=255)
 cTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGac                                                                                                                                               <  6:100028/148‑1 (MQ=255)
  tAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGaca                                                                                                                                              >  6:281168/1‑148 (MQ=255)
  tAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGaca                                                                                                                                              >  7:140604/1‑148 (MQ=255)
            tgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGAt                                                                                                                                          >  8:316571/1‑142 (MQ=255)
            tgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGAt                                                                                                                                          <  7:316571/142‑1 (MQ=255)
             gaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTTATTTGTCGCAATGATTGACACGATTCCGCt                                                                                                                                    <  6:206261/147‑1 (MQ=255)
              aaaaTCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCtt                                                                                                                                   <  2:206354/147‑1 (MQ=255)
                 aTCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                                                                                               >  8:353703/1‑148 (MQ=255)
                         tttATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTa                                                                                                                        <  5:281168/147‑1 (MQ=255)
                          ttATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGGCACGATTCCGCTTGACGCTGCGTaa                                                                                                                       <  1:407299/147‑1 (MQ=255)
                             tCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTa                                                                                                                        >  5:390422/1‑143 (MQ=255)
                             tCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTa                                                                                                                        <  6:390422/143‑1 (MQ=255)
                               gTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGtttt                                                                                                                 <  5:78885/148‑1 (MQ=255)
                                      cccTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGAAGCTGCGTAAGGTTTTTGTAAtt                                                                                                          >  4:383655/1‑148 (MQ=255)
                                        cTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAAtttt                                                                                                        <  3:383655/148‑1 (MQ=255)
                                          ttAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTa                                                                                                       >  6:192476/1‑147 (MQ=255)
                                          ttAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTa                                                                                                       >  1:203318/1‑147 (MQ=255)
                                           tAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTa                                                                                                       >  7:247161/1‑146 (MQ=255)
                                           tAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTAc                                                                                                      <  8:140604/147‑1 (MQ=255)
                                            aaaaTTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAg                                                                                                    >  3:202404/1‑148 (MQ=255)
                                              aaTTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTa                                                                                                       >  1:7922/1‑143 (MQ=255)
                                                ttCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCaa                                                                                                <  2:203318/148‑1 (MQ=255)
                                                   ggggCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCt                                                                                             <  5:410940/148‑1 (MQ=255)
                                                     ggCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCttt                                                                                           <  2:221352/148‑1 (MQ=255)
                                                     ggCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCttt                                                                                           <  5:192476/148‑1 (MQ=255)
                                                      gcgcCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCtttt                                                                                          <  8:199326/148‑1 (MQ=255)
                                                        gcCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTAt                                                                                        <  5:85690/148‑1 (MQ=255)
                                                               ccATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTaa                                                                                 <  7:307962/148‑1 (MQ=255)
                                                                 aTGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAAc                                                                                >  2:1072/1‑147 (MQ=255)
                                                                   gtgGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCtt                                                                                            >  1:12551/1‑133 (MQ=255)
                                                                   gtgGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCtt                                                                                            <  2:12551/133‑1 (MQ=255)
                                                                      gTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATa                                                                           <  7:230430/147‑1 (MQ=255)
                                                                       tctcAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGc                                                                         >  2:190997/1‑148 (MQ=255)
                                                                         tcAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtg                                                                       >  6:171966/1‑148 (MQ=255)
                                                                          cAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGCCGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtgg                                                                      >  5:304796/1‑148 (MQ=255)
                                                                          cAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtgg                                                                      >  5:346993/1‑148 (MQ=255)
                                                                          cAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtgg                                                                      >  3:109598/1‑148 (MQ=255)
                                                                           aaGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtggt                                                                     <  1:1072/148‑1 (MQ=255)
                                                                           aaGCACAAAGGAAGCGTGAGGCGAGTCAGTCGCGTAATGCGTAGGCAAAGGATTGATTTGTCGCAAAGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtggt                                                                     <  4:219556/148‑1 (MQ=255)
                                                                            aGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtggtg                                                                    <  1:190997/148‑1 (MQ=255)
                                                                              cccAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGTATTGATTTGTAGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGa                                                                  <  5:344939/148‑1 (MQ=255)
                                                                              cccAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtggtgg                                                                   >  6:344939/1‑147 (MQ=255)
                                                                              cccAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtggtgg                                                                   >  3:367223/1‑147 (MQ=255)
                                                                               ccAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGaa                                                                 <  7:438665/148‑1 (MQ=255)
                                                                                  aaGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAata                                                               >  4:136024/1‑147 (MQ=255)
                                                                                  aaGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAata                                                               >  5:40912/1‑147 (MQ=255)
                                                                                   aGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTACGCTTGCCGCTGCGTAAGGTTTTTGTAATTTTACAGGCCACCTTTTATTCACTAACAAATAGCTGGTGGAatata                                                             >  2:142491/1‑148 (MQ=255)
                                                                                     gAAGATTTAGGCGAGGCAGTCGCGTAATGCTTACGCACAGGATTGATTTGTCGCAATGATAGACACGATTCCGCTTGACGCTGCGTAATTTTTTTGTAATTTTACATGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATg                                                           <  4:398320/148‑1 (MQ=255)
                                                                                     gAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATg                                                           <  4:202404/148‑1 (MQ=255)
                                                                                            gAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCTATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATc                                                     >  2:274824/1‑147 (MQ=255)
                                                                                              ggCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCAGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCaaa                                                  >  5:23058/1‑148 (MQ=255)
                                                                                              gccGAGTCAGGCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGAGGCTGCGTTAGGTTGTTGTAATTTTACAGGCACCCTTTTATTCACTATCAAATCGCGGGTGGAATATATGACTATCa                                                    >  6:139322/3‑146 (MQ=255)
                                                                                                 gAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGAGTATCAAAGTa                                               <  7:53841/148‑1 (MQ=255)
                                                                                                  agtcagtcGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAg                                              <  3:6102/148‑1 (MQ=255)
                                                                                                   gtcagtcGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAgg                                             <  5:171966/148‑1 (MQ=255)
                                                                                                    tcagtcGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTaa                                                                                 <  6:428383/111‑1 (MQ=255)
                                                                                                    tcagtcGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTaa                                                                                 >  5:428383/1‑111 (MQ=255)
                                                                                                    tcagtcGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGt                                            <  4:109598/148‑1 (MQ=255)
                                                                                                       gtcGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTAt                                          >  5:250590/1‑147 (MQ=255)
                                                                                                        tcGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTa                                                                                  <  3:46751/106‑1 (MQ=255)
                                                                                                        tcGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTa                                                                                  >  4:46751/1‑106 (MQ=255)
                                                                                                         cgcgTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACt                                                        <  5:376029/131‑1 (MQ=255)
                                                                                                         cgcgTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACt                                                        >  6:376029/1‑131 (MQ=255)
                                                                                                         cgcgTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCaa                                       >  1:231973/1‑148 (MQ=255)
                                                                                                         cgcgTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCaa                                       >  2:28171/1‑148 (MQ=255)
                                                                                                           cgTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACg                                     >  6:102788/1‑148 (MQ=255)
                                                                                                             tAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGt                                   >  1:284556/1‑148 (MQ=255)
                                                                                                              aaTGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTc                                                                                      <  5:255588/96‑1 (MQ=255)
                                                                                                              aaTGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTc                                                                                      >  6:255588/1‑96 (MQ=255)
                                                                                                              aaTGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACaa                                                                              >  6:63774/1‑104 (MQ=255)
                                                                                                              aaTGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACaa                                                                              <  5:63774/104‑1 (MQ=255)
                                                                                                               aTGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGttt                                 <  6:40912/148‑1 (MQ=255)
                                                                                                                 gCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTa                                                                                                             >  4:164638/1‑70 (MQ=255)
                                                                                                                 gCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTa                                                                                                             <  3:164638/70‑1 (MQ=255)
                                                                                                                 gCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTa                                                                                                       <  3:68443/76‑1 (MQ=255)
                                                                                                                 gCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTa                                                                                                       >  4:68443/1‑76 (MQ=255)
                                                                                                                 gCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCttt                                                                                           >  1:236990/1‑88 (MQ=255)
                                                                                                                 gCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCttt                                                                                           <  2:236990/88‑1 (MQ=255)
                                                                                                                    tAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGtt                                  >  3:235098/1‑142 (MQ=255)
                                                                                                                        cacaGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGTTTTGGCCGTAt                        >  6:372391/1‑148 (MQ=255)
                                                                                                                        cacaGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGTTTTGGCCGTAt                        >  4:239005/1‑148 (MQ=255)
                                                                                                                           aGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAatat                                                              >  6:203137/1‑107 (MQ=255)
                                                                                                                           aGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAatat                                                              <  5:203137/107‑1 (MQ=255)
                                                                                                                           aGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGTTTTGGCCGTATCgg                     >  3:325528/1‑148 (MQ=255)
                                                                                                                            ggATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGTTTTGGCCGTATCGGt                    <  3:136024/148‑1 (MQ=255)
                                                                                                                              attgattTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGTTTTGGCCGTATCGGTCg                  >  5:415248/1‑148 (MQ=255)
                                                                                                                              attgattTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGTTTTGGCCGTATCGGTCg                  >  5:211396/1‑148 (MQ=255)
                                                                                                                              attgattTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGTTTTGGCCGTATCGGTCg                  >  7:79878/1‑148 (MQ=255)
                                                                                                                                  attTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGTTTTGGCCGTATCGGTCGCAtt              >  8:340862/1‑148 (MQ=255)
                                                                                                                                  attTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGTTTTGGCCGTATCGGTCGCAt               <  3:239005/147‑1 (MQ=255)
                                                                                                                                        cGCAATGATTGGCACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGTTTTGGCCGTATCGGTCGCATTGTTTTc        >  3:244864/1‑148 (MQ=255)
                                                                                                                                        cGCAATGATTGGCACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGTTTTGGCCGTATCGGTCGCATTGTTTTc        >  3:250349/1‑148 (MQ=255)
                                                                                                                                          cAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGtt                                  >  4:323615/1‑120 (MQ=255)
                                                                                                                                          cAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGtt                                  <  3:323615/120‑1 (MQ=255)
                                                                                                                                           aaTGATTGGCACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGTTTTGGCCGTATCGGTCGCATTGTTTTCCGt     <  4:244864/148‑1 (MQ=255)
                                                                                                                                           aaTGATTGGCACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGTTTTGGCCGTATCGGTCGCATTGTTTTCCGt     <  4:250349/148‑1 (MQ=255)
                                                                                                                                              gATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGTTTTGGCCGTAGTGGTCGCATTGTTTTCCGtgct  >  5:223350/1‑148 (MQ=255)
                                                                                                                                                   acacGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATg                                                           >  7:23360/1‑86 (MQ=255)
                                                                                                                                                   acacGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATg                                                           <  8:23360/86‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGAATATATGACTATCAAAGTAGGTATCAACGGTTTTGGCCGTATCGGTCGCATTGTTTTCCGTGCT  >  NC_000913/1862541‑1862830

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