Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,708,823 T→C R128R (CGT→CGC nadB → quinolinate synthetase B protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,708,8230TC97.6% 505.0 / 0.7 166R128R (CGT→CGCnadBquinolinate synthetase B protein
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (91/71);  total (92/74)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 178 total aligned reads displayed.

CGACCAGGGGGTGTTGTTTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGTATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCGAATATTCGCGTGCTGGAGC  >  NC_000913/2708733‑2708914
                                                                                          |                                                                                           
  aCCAGGGGGTGTTGTTTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCAttc                                                                                         <  1:3016358/93‑1 (MQ=255)
   ccAGGGGGTGTTGTTTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCAttct                                                                                        <  1:3944688/93‑1 (MQ=255)
    cAGGGGGTGTTGTTTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCAttctt                                                                                       <  1:3309911/93‑1 (MQ=255)
    cAGGGGGTGTTGTTTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCAttctt                                                                                       >  1:4663901/1‑93 (MQ=255)
    cAGGGGGTGTTGTTTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCAttctt                                                                                       >  2:2464220/1‑93 (MQ=255)
    cAGGGGGTGTTGTTTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCAttctt                                                                                       >  2:1599362/1‑93 (MQ=255)
     aGGGGGTGTTGTTTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCAttcttc                                                                                      <  1:2668347/93‑1 (MQ=255)
      gggggTGTTGTTTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCa                                                                                     >  2:2632179/1‑93 (MQ=255)
       ggggTGTTGTTTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCAt                                                                                    >  2:124196/1‑93 (MQ=255)
       ggggTGTTGTTTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCAt                                                                                    >  2:2170413/1‑93 (MQ=255)
        gggTGTTGTTTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCAtc                                                                                   >  1:1569270/1‑92 (MQ=255)
         ggTGTTGTTTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATgc                                                                                  <  1:1039053/93‑1 (MQ=255)
          gtgtTGTTTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATgcc                                                                                 >  1:2603313/1‑93 (MQ=255)
          gtgtTGTTTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATgcc                                                                                 <  2:2051268/93‑1 (MQ=255)
            gttgttTGATACCCACATTCAACCGAATGGTGAAGAAAGTTATCATCTGACCCGTGAATGTGGACATAGTCCTCGTCGCATTCTTCTTgccgc                                                                               >  2:1994288/1‑93 (MQ=255)
            gttgttTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATgccgc                                                                               <  2:2190476/93‑1 (MQ=255)
            gttgttTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGACACCGTCGCATTCTTCCTgccgc                                                                               >  1:3846536/1‑93 (MQ=255)
               ttttGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGa                                                                            <  2:3944939/92‑1 (MQ=255)
               gttTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGccgccgt                                                                            >  1:3930311/1‑92 (MQ=255)
               gttTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGa                                                                            <  1:962265/93‑1 (MQ=255)
               gttTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGa                                                                            >  2:3697531/1‑93 (MQ=255)
                 ttGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACg                                                                          >  2:1182589/1‑93 (MQ=255)
                    aTACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCa                                                                       >  1:2398537/1‑93 (MQ=255)
                     tACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCAc                                                                      <  2:2800755/93‑1 (MQ=255)
                      aCCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCAcc                                                                     <  1:4165245/93‑1 (MQ=255)
                       cccACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCg                                                                    >  2:2217752/1‑93 (MQ=255)
                       cccACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCg                                                                    <  2:1231614/93‑1 (MQ=255)
                       cccACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCg                                                                    <  1:3206261/93‑1 (MQ=255)
                       cccACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCg                                                                    >  2:3118690/1‑93 (MQ=255)
                        ccACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCgg                                                                   <  1:4547679/93‑1 (MQ=255)
                        ccACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCgg                                                                   >  2:884693/1‑93 (MQ=255)
                           caTTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTag                                                                <  2:4642881/93‑1 (MQ=255)
                           caTTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTag                                                                <  2:4323447/93‑1 (MQ=255)
                            aTTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTaga                                                               >  2:1887327/1‑93 (MQ=255)
                            aTTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTaga                                                               <  1:2225470/93‑1 (MQ=255)
                            aTTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTaga                                                               <  1:1729481/93‑1 (MQ=255)
                            aTTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTaga                                                               >  2:3479075/1‑93 (MQ=255)
                             ttCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTagag                                                              >  1:279099/1‑93 (MQ=255)
                             ttCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTagag                                                              >  2:1995400/1‑93 (MQ=255)
                              tCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAgaga                                                             <  1:4219841/93‑1 (MQ=255)
                                aaCCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAg                                                           <  2:50955/93‑1 (MQ=255)
                                aaCCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAg                                                           >  2:4145402/1‑93 (MQ=255)
                                 aCCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGt                                                          <  2:2850208/93‑1 (MQ=255)
                                 aCCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGt                                                          >  2:460096/1‑93 (MQ=255)
                                 aCCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGt                                                          <  1:4012042/93‑1 (MQ=255)
                                  ccGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTa                                                         >  2:1815520/1‑93 (MQ=255)
                                   cGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAg                                                        <  1:726863/93‑1 (MQ=255)
                                    gAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGa                                                       >  2:3355295/1‑93 (MQ=255)
                                      aTGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCTCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGaaa                                                     <  1:2500883/93‑1 (MQ=255)
                                      aTGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGaaa                                                     >  1:2456594/1‑93 (MQ=255)
                                       tGGCGAATAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCGGACGCCACCGGTAGAGAAGTAGAAAc                                                    <  1:1932667/93‑1 (MQ=255)
                                       tGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGCTAGAGAAGTAGaaaa                                                    >  1:4526979/1‑92 (MQ=255)
                                        ggCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAAcc                                                   <  1:1084498/93‑1 (MQ=255)
                                         gCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGTATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCa                                                  <  1:237507/93‑1 (MQ=255)
                                          cGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCAc                                                 <  1:3331755/93‑1 (MQ=255)
                                          cGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCAc                                                 >  2:1279219/1‑93 (MQ=255)
                                          cGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCAc                                                 >  1:151164/1‑93 (MQ=255)
                                          cGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCAc                                                 <  1:2184530/93‑1 (MQ=255)
                                           gaagaaAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACg                                                >  2:2715471/1‑93 (MQ=255)
                                            aagaaAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCCTTCTTCATGCCGCCGACGCCACCGGTAGAGAAGCACAAACCAAGc                                               >  2:4012318/1‑93 (MQ=255)
                                            aagaaACGTACCATCTGACCCCTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGc                                               <  1:843313/93‑1 (MQ=255)
                                             agaaAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCt                                              <  2:3589744/93‑1 (MQ=255)
                                             agaaAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCt                                              >  1:1802866/1‑93 (MQ=255)
                                              gaaAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTg                                             >  1:4114038/1‑93 (MQ=255)
                                               aaaTTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTgg                                            <  2:1750777/93‑1 (MQ=255)
                                                  gTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGa                                         <  2:657050/93‑1 (MQ=255)
                                                   ttACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAg                                        >  1:290967/1‑93 (MQ=255)
                                                    tACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGc                                       >  2:1367274/1‑93 (MQ=255)
                                                    tACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGc                                       >  2:4096381/1‑93 (MQ=255)
                                                     aCCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCa                                      >  2:2479628/1‑93 (MQ=255)
                                                     aCCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGAAGGGAAGTAGAAAACACGCTGGTGAGCa                                      >  2:180960/1‑93 (MQ=255)
                                                      ccATCTGACCGGAGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCaa                                     <  2:3134756/93‑1 (MQ=255)
                                                       cATCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAg                                    <  2:303937/93‑1 (MQ=255)
                                                        aTCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAgg                                   <  1:977632/93‑1 (MQ=255)
                                                        aTCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAgg                                   <  2:2851644/93‑1 (MQ=255)
                                                         tCTGACCCGTGAAGGTGGACATAGTCCCCGTCTCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGgc                                  <  1:3108279/93‑1 (MQ=255)
                                                         tCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGgc                                  >  1:4634072/1‑93 (MQ=255)
                                                         tCTGACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGgc                                  >  2:2959196/1‑93 (MQ=255)
                                                            gACCCGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCt                               >  2:1934569/1‑93 (MQ=255)
                                                              cccGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGa                             >  1:144052/1‑93 (MQ=255)
                                                               ccGTGAAGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTTAGCAAGGCGCTGac                            >  1:2805681/1‑92 (MQ=255)
                                                                    aaNGTGNACATAGTCACCGTCGCATTCTTCATGCCGCCGACCCCACCGGTAGAGAAGTAGAACCCACGCTGGTGAGCAAGGCGCTGAACCATc                       <  1:1358180/93‑1 (MQ=255)
                                                                     aGGTGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATcc                      <  2:1323639/93‑1 (MQ=255)
                                                                      ggTGGACATAGTCCCCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCg                     <  1:3442176/93‑1 (MQ=255)
                                                                        tGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCGaa                   >  1:3406361/1‑93 (MQ=255)
                                                                        tGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCGaa                   <  1:4049682/93‑1 (MQ=255)
                                                                        tGGACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCGaa                   <  2:2232639/93‑1 (MQ=255)
                                                                          gACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCGAata                 >  1:3612378/1‑93 (MQ=255)
                                                                          gACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCGAata                 >  2:3233340/1‑93 (MQ=255)
                                                                          gACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCGAata                 >  1:3213502/1‑93 (MQ=255)
                                                                          gACATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCGAata                 <  2:15113/93‑1 (MQ=255)
                                                                            cATAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCGAATAtt               >  2:2497557/1‑93 (MQ=255)
                                                                             aTAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCGAATATTc              <  1:4121359/93‑1 (MQ=255)
                                                                             aTAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAAACACGCTGGTGAGCAAGGCGCTGAACCATCCGAATATTc              >  2:525261/1‑93 (MQ=255)
                                                                              tAGTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCGAATATTcg             >  2:2272078/1‑93 (MQ=255)
                                                                                gTCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCGAATATTtgcg           >  1:3260886/1‑93 (MQ=255)
                                                                                 tCACCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCGAATATTCGCGt          >  2:1380203/1‑93 (MQ=255)
                                                                                 gccccGTCGTATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCGAATATTCGCGt          <  1:880225/90‑1 (MQ=255)
                                                                                   aCCGTCGCATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCGAATATTCGCGTGc        <  1:601213/93‑1 (MQ=255)
                                                                                       tcgcATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCGAATATTCGCGTGCTGGa    <  2:833970/93‑1 (MQ=255)
                                                                                          |                                                                                           
CGACCAGGGGGTGTTGTTTGATACCCACATTCAACCGAATGGCGAAGAAAGTTACCATCTGACCCGTGAAGGTGGACATAGTCACCGTCGTATTCTTCATGCCGCCGACGCCACCGGTAGAGAAGTAGAAACCACGCTGGTGAGCAAGGCGCTGAACCATCCGAATATTCGCGTGCTGGAGC  >  NC_000913/2708733‑2708914

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