Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,722,300 G→T 16.2% intergenic (+251/+28) insK → / ← glyS IS150 transposase B/glycine tRNA synthetase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,722,3000GT16.2% 611.9 / 96.7 232intergenic (+251/+28)insK/glySIS150 transposase B/glycine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (87/105);  new base T (20/17);  total (110/122)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-01

Only 100 of 237 total aligned reads displayed.

CCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTT  >  NC_000913/3722201‑3722399
                                                                                                   |                                                                                                   
ccTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAg                                                                                                     <  2:2993224/100‑1 (MQ=255)
  tAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGAt                                                                                                   >  1:3100360/1‑100 (MQ=255)
  tAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCATGCCAGAt                                                                                                   >  1:2944980/1‑100 (MQ=255)
    aaaaaaTCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATg                                                                                                  >  1:2144228/1‑99 (MQ=255)
      aaaaTCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCa                                                                                               <  1:1178698/100‑1 (MQ=255)
       aaaTCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAg                                                                                              <  2:2288855/100‑1 (MQ=255)
         aTCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAgn                                                                                             <  1:1929620/99‑2 (MQ=255)
         aTCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGc                                                                                            >  2:3989703/1‑100 (MQ=255)
         aTCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGc                                                                                            <  1:1817356/100‑1 (MQ=255)
         aTCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGc                                                                                            <  1:219242/100‑1 (MQ=255)
            tcATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCttt                                                                                         <  2:1937066/100‑3 (MQ=255)
             cATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCtttt                                                                                        <  2:932994/100‑1 (MQ=255)
             cATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAtt                                                                                        >  1:887823/1‑100 (MQ=255)
             cATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAtt                                                                                        <  2:2411259/100‑1 (MQ=255)
              aTAAAATTAATATATGGGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGTTTTATTATTTAGCGATAAAAAAAGCTTGCAATCTGGCAGGCTtttt                                                                                       <  1:3387112/100‑1 (MQ=255)
              aTAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCTtttt                                                                                       >  1:1526978/1‑100 (MQ=255)
              aTAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCTtttt                                                                                       <  2:4164482/100‑1 (MQ=255)
              aTAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAttt                                                                                       <  2:2063419/100‑1 (MQ=255)
              aTAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAttt                                                                                       >  1:1785934/1‑100 (MQ=255)
              aTAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAttt                                                                                       <  1:1599582/100‑1 (MQ=255)
              aTAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAt                                                                                         <  1:1849732/98‑1 (MQ=255)
              aTAAAAATAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAAATAGCAGGCTAttt                                                                                       >  1:296636/1‑100 (MQ=255)
               tAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTta                                                                                      >  2:3756406/1‑100 (MQ=255)
                aaaaTTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTtaa                                                                                     <  2:3390726/100‑1 (MQ=255)
                aaaaTTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTtaa                                                                                     <  2:1849057/100‑1 (MQ=255)
                aaaaTTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCAACATGACAGGCTATTtaa                                                                                     >  1:2566938/1‑100 (MQ=255)
                 aaaTTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTtaat                                                                                    >  2:1187472/1‑100 (MQ=255)
                  aaTTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCTttttttat                                                                                   <  2:2787915/100‑5 (MQ=38)
                  aaTTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTtaata                                                                                   <  2:1656211/100‑1 (MQ=255)
                    ttAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAAc                                                                                 <  1:2008070/100‑1 (MQ=255)
                     tAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCTttttttatcgc                                                                                <  1:2085937/100‑8 (MQ=37)
                     tAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCTttttttatcgc                                                                                <  2:1526983/100‑8 (MQ=37)
                        tatatGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCg                                                                             <  1:748134/100‑1 (MQ=255)
                          tatGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGtt                                                                           >  2:3613902/1‑100 (MQ=255)
                          tatGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGtt                                                                           <  1:2871360/100‑1 (MQ=255)
                           atGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGtta                                                                          >  2:2637036/1‑100 (MQ=255)
                              agaTAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGttatta                                                                       >  2:3364708/1‑100 (MQ=255)
                              agaTAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGttatta                                                                       <  1:2176671/100‑1 (MQ=255)
                              agaTAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGttatta                                                                       <  2:4088773/100‑1 (MQ=255)
                                 tAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGCCGttattaaat                                                                    <  2:3946447/100‑4 (MQ=255)
                                  aaTCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGCCGttattaaata                                                                   <  2:2508170/100‑5 (MQ=255)
                                    tCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCaa                                                                 >  2:2676811/1‑100 (MQ=255)
                                     cTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGCCGttattaaatagcc                                                                <  2:3129329/100‑8 (MQ=255)
                                       ttATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGCCGttattaaatagcctg                                                              <  2:215635/100‑10 (MQ=39)
                                       ttATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGCCGttattaaatagcc                                                                <  1:2662530/98‑8 (MQ=255)
                                        tATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGc                                                             >  2:2202410/1‑100 (MQ=255)
                                        tATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGCCGttattaaatagcctgc                                                             <  1:3901404/100‑11 (MQ=39)
                                            cagcagAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGaaa                                                         >  2:2313594/1‑100 (MQ=255)
                                            cagcagAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGaaa                                                         <  1:4215741/100‑1 (MQ=255)
                                            cagcagAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGaaa                                                         >  2:2959100/1‑100 (MQ=255)
                                                 gaagaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATcc                                                    <  2:1127432/100‑1 (MQ=255)
                                                 gaagaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATcc                                                    <  1:65078/100‑1 (MQ=255)
                                                   agaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGc                                                  >  1:3553842/1‑100 (MQ=255)
                                                    gaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCa                                                 >  2:1022025/1‑100 (MQ=255)
                                                      ttattaAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAAc                                               >  1:3678913/1‑100 (MQ=255)
                                                        attaAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAAcgc                                             <  2:3876360/100‑1 (MQ=255)
                                                           aaaGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCa                                          >  2:1677986/1‑100 (MQ=255)
                                                             aGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATGTCCGCAACGCGCAgg                                        >  2:4338280/1‑100 (MQ=255)
                                                             aGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGACTATTTCATAACGGCGTTAATATTGCAACAGCGAAATATCCACAACGCGCAgg                                        <  1:3073274/100‑1 (MQ=255)
                                                                ttGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAAc                                     <  2:672412/100‑1 (MQ=255)
                                                                 tGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACa                                    <  2:3255417/100‑1 (MQ=255)
                                                                  gCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAg                                   >  2:1402160/1‑100 (MQ=255)
                                                                  gCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAg                                   <  1:10077/100‑1 (MQ=255)
                                                                   cTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGt                                  >  2:2195304/1‑100 (MQ=255)
                                                                   cTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGt                                  >  1:4020763/1‑100 (MQ=255)
                                                                     gTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATACGCAACGCGCAGGAACAGtt                                 >  2:4386935/1‑99 (MQ=255)
                                                                      tattatTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTcg                               >  1:1394962/1‑100 (MQ=255)
                                                                         tattTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCgcgc                            <  2:592646/100‑1 (MQ=255)
                                                                         tattTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCgcgc                            >  2:2015726/1‑100 (MQ=255)
                                                                                  aTAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTcc                   <  2:3532616/100‑1 (MQ=255)
                                                                                  aTAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTcc                   >  1:44300/1‑100 (MQ=255)
                                                                                    aaaaaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAg                 <  2:1563822/100‑1 (MQ=255)
                                                                                    aaaaaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAg                 <  2:104518/100‑1 (MQ=255)
                                                                                     taaataGCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGc                >  1:4077106/2‑100 (MQ=255)
                                                                                     taaataGCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATACGCAACGCGCCGGAACAGTTCGCGCAGTTTCTCCAGc                >  1:777263/2‑100 (MQ=38)
                                                                                     aacaaaGCCTGCCAGATGGAAGGCTATTTAATAACAGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGATTCTCCa                  >  1:1461300/4‑98 (MQ=255)
                                                                                     aaaaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGc                >  1:3258510/1‑100 (MQ=255)
                                                                                      aaataGCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCa               <  2:609467/96‑1 (MQ=255)
                                                                                      aaataGCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCa               >  1:609465/1‑100 (MQ=255)
                                                                                      aaaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCa               >  2:80573/1‑100 (MQ=255)
                                                                                       aaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCAt              <  2:1003562/100‑1 (MQ=255)
                                                                                         taGCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATgg            >  2:3901389/2‑100 (MQ=255)
                                                                                         taGCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATgg            >  1:2463312/2‑100 (MQ=255)
                                                                                         aaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATgg            <  2:3678935/100‑1 (MQ=255)
                                                                                         aaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATgg            <  1:1342656/100‑1 (MQ=255)
                                                                                          aGCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGt           >  1:2474902/1‑100 (MQ=255)
                                                                                           gCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGt           >  1:636166/1‑99 (MQ=255)
                                                                                           gCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTc          >  1:3736277/1‑100 (MQ=255)
                                                                                           gCCTGCCATCTGGCAGGCTAATTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCACGAACAGTTCGCGCAGTTTCTCCAGCATgg            >  1:2696338/1‑98 (MQ=255)
                                                                                           gCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTc          <  1:4321716/100‑1 (MQ=255)
                                                                                            ccTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCa         <  2:1765374/100‑1 (MQ=255)
                                                                                            ccTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCa         <  1:4137154/100‑1 (MQ=255)
                                                                                            ccTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCa         <  2:2862220/100‑1 (MQ=255)
                                                                                             ctcccAGATGGCAGGCTATTTAATAACGGCGCTATTATTGCAACAGCCAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATgg            <  1:3414039/93‑1 (MQ=255)
                                                                                             cTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAg        <  2:693583/100‑1 (MQ=255)
                                                                                               gCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGAc      <  1:1651100/100‑1 (MQ=255)
                                                                                               gCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGAc      >  1:4350368/1‑100 (MQ=255)
                                                                                               gCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGAc      >  2:950477/1‑100 (MQ=255)
                                                                                                ccAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACg     >  2:2820795/1‑100 (MQ=255)
                                                                                                   gATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGtt  <  2:3183868/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTT  >  NC_000913/3722201‑3722399

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