Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,722,301 A→C 15.8% intergenic (+252/+27) 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,3010AC15.8% 575.5 / 103.8 235intergenic (+252/+27)insK/glySIS150 transposase B/glycine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (91/107);  new base C (20/17);  total (111/124)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-01

Only 100 of 236 total aligned reads displayed.

CTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTT  >  NC_000913/3722202‑3722399
                                                                                                   |                                                                                                  
 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)
        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)
           tcATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAt                                                                                         <  1:592857/100‑1 (MQ=255)
            cATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCtttt                                                                                        <  2:932994/100‑1 (MQ=255)
            cATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAtt                                                                                        >  1:2528723/1‑100 (MQ=255)
            cATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAtt                                                                                        >  1:1883438/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:2784015/100‑1 (MQ=255)
             aTAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAttt                                                                                       >  1:1785934/1‑100 (MQ=255)
             aTAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAttt                                                                                       <  2:2829800/100‑1 (MQ=255)
             aTAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTAttt                                                                                       <  1:2860803/100‑1 (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)
               aaaaTTAATATTTGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCTtttttt                                                                                     <  2:1756300/100‑3 (MQ=38)
               aaaaTTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTtaa                                                                                     <  2:2914411/100‑1 (MQ=255)
               aaaaTTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTtaa                                                                                     <  2:3390726/100‑1 (MQ=255)
                aaaTTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTtaat                                                                                    >  2:1187472/1‑100 (MQ=255)
                   ttAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAAc                                                                                 <  1:2008070/100‑1 (MQ=255)
                    tAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCTttttttatcgc                                                                                <  2:1526983/100‑8 (MQ=37)
                    tAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCATCTGGCAGGCTttttttatcgc                                                                                <  1:2085937/100‑8 (MQ=37)
                         tatGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGtt                                                                           >  2:3613902/1‑100 (MQ=255)
                         tatGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGtt                                                                           >  2:1297551/1‑100 (MQ=255)
                             agaTAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGttatta                                                                       >  2:3364708/1‑100 (MQ=255)
                             agaTAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGttatta                                                                       <  1:2176671/100‑1 (MQ=255)
                             agaTAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGCCGttatta                                                                       <  2:1040589/100‑1 (MQ=255)
                                tAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTg                                                                    <  1:3540053/100‑1 (MQ=255)
                                tAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTg                                                                    <  2:3592819/100‑1 (MQ=255)
                                tAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGCCGttattaaat                                                                    <  2:2309221/100‑4 (MQ=255)
                                tAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGCCGttattaaat                                                                    <  2:3946447/100‑4 (MQ=255)
                                  aTCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCa                                                                  >  2:218147/1‑100 (MQ=255)
                                  aTCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCa                                                                  >  2:4099108/1‑100 (MQ=255)
                                   tCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCaa                                                                 >  2:2676811/1‑100 (MQ=255)
                                    cTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGCCGttattaaatagcc                                                                <  2:2258096/100‑8 (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                                                             <  2:2505462/100‑11 (MQ=38)
                                           cagcagAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGaaa                                                         <  1:4215741/100‑1 (MQ=255)
                                                gaagaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATcc                                                    <  2:1127432/100‑1 (MQ=255)
                                                gaagaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATcc                                                    <  1:65078/100‑1 (MQ=255)
                                                  agaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGc                                                  <  2:3054746/100‑1 (MQ=255)
                                                   gaTTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCa                                                 >  2:1022025/1‑100 (MQ=255)
                                                     ttattaAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAAc                                               >  1:3678913/1‑100 (MQ=255)
                                                       attaAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAAcgc                                             <  2:3876360/100‑1 (MQ=255)
                                                       attaAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAAcgc                                             >  2:2807617/1‑100 (MQ=255)
                                                          aaaGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCa                                          >  2:1677986/1‑100 (MQ=255)
                                                            aGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATGTCCGCAACGCGCAgg                                        >  2:4338280/1‑100 (MQ=255)
                                                            aGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAgg                                        <  2:2039233/100‑1 (MQ=255)
                                                            aGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGACTATTTCATAACGGCGTTAATATTGCAACAGCGAAATATCCACAACGCGCAgg                                        <  1:3073274/100‑1 (MQ=255)
                                                               ttGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAAc                                     <  2:672412/100‑1 (MQ=255)
                                                                 gCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAg                                   <  1:10077/100‑1 (MQ=255)
                                                                  cTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGt                                  >  1:4020763/1‑100 (MQ=255)
                                                                    gTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTc                                >  1:3841924/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)
                                                                         attTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGAGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCa                           >  1:4131189/1‑100 (MQ=255)
                                                                          tttAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAg                          <  2:2012535/100‑1 (MQ=255)
                                                                            tAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGtt                        <  1:975958/100‑1 (MQ=255)
                                                                              gCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAAGATCCGCAACGCGCCGGAACAGTTCACCCAGTTtc                      >  1:691761/1‑100 (MQ=255)
                                                                               cGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTtct                     >  2:3310030/1‑100 (MQ=255)
                                                                                 aTAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTcc                   >  1:44300/1‑100 (MQ=255)
                                                                                   aaaaaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAg                 <  2:1563822/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               >  2:2549649/1‑100 (MQ=255)
                                                                                     aaaaaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCa               >  2:80573/1‑100 (MQ=255)
                                                                                        taGCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATgg            >  2:3901389/2‑100 (MQ=255)
                                                                                        taGCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATgg            >  1:2463312/2‑100 (MQ=255)
                                                                                        aaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATgg            <  1:1342656/100‑1 (MQ=255)
                                                                                        aaGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATgg            <  2:3678935/100‑1 (MQ=255)
                                                                                         aGCCTGCCATCTGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGt           >  1:2474902/1‑100 (MQ=255)
                                                                                         aGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGt           >  2:1653073/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)
                                                                                           ccTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCa         <  2:1768968/100‑1 (MQ=255)
                                                                                            ctcccAGATGGCAGGCTATTTAATAACGGCGCTATTATTGCAACAGCCAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATgg            <  1:3414039/93‑1 (MQ=255)
                                                                                            cTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAg        <  2:693583/100‑1 (MQ=255)
                                                                                              gCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGAc      >  2:950477/1‑100 (MQ=255)
                                                                                              gCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGAc      <  1:1651100/100‑1 (MQ=255)
                                                                                              gCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGAc      >  1:4350368/1‑100 (MQ=255)
                                                                                                  gATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGtt  <  2:3183868/100‑1 (MQ=255)
                                                                                                  gATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGtt  >  1:4284009/1‑100 (MQ=255)
                                                                                                  gATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGAAGTTTCTCCAGCATGGTCAGACGGtt  >  1:881020/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
CTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGGCGTTATTATTGCAACAGCGAAATATCCGCAACGCGCAGGAACAGTTCGCGCAGTTTCTCCAGCATGGTCAGACGGTT  >  NC_000913/3722202‑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: