Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,436,497 C→A 5.9% G133C (GGC→TGC)  ompN ← outer membrane pore protein N, non‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,436,4970CA 5.9% 219.5 / 3.1 68G133C (GGC→TGC) ompNouter membrane pore protein N, non‑specific
Reads supporting (aligned to +/- strand):  ref base C (33/31);  new base A (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

ATTACTGGCTCCTTCGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTGAAGATTCAGTGTTGTTT  >  NC_000913/1436351‑1436645
                                                                                                                                                  |                                                                                                                                                    
aTTACTGGCTCCTTCGTTGATACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCaa                                                                                                                                                    <  3:158254/149‑1 (MQ=255)
 ttACTGGCTCCTTCGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCaaa                                                                                                                                                   >  5:78598/1‑149 (MQ=255)
     tgacTCCTTCGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAAtt                                                                                                                                                 >  3:76503/4‑147 (MQ=255)
       gCTCCTTCGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCa                                                                                                                                               <  5:303647/147‑1 (MQ=255)
       gCTCCTTCGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAgg                                                                                                                                             <  3:328277/149‑1 (MQ=255)
          ccTTCGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGgcag                                                                                                                                          >  3:259763/1‑149 (MQ=255)
            ttCGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTCCATTGTTATAAGAGTCACAGCCAAATTCAGgcagca                                                                                                                                        >  8:150665/1‑149 (MQ=255)
             tCGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGgcagc                                                                                                                                         <  5:99695/147‑1 (MQ=255)
             tCGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGgcagc                                                                                                                                         >  6:99695/1‑147 (MQ=255)
                  gttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCg                                                                                                                                   <  4:249635/148‑1 (MQ=255)
                     aCCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTc                                                                                                                                <  1:179218/148‑1 (MQ=255)
                     aCCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTc                                                                                                                                >  3:38235/1‑148 (MQ=255)
                      ccTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGAATATCTGTccc                                                                                                                              >  2:222633/1‑148 (MQ=255)
                       cTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAg                                                                                                                             >  3:222216/1‑149 (MQ=255)
                       cTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCAAAATTCAGGCAGCATATCGGTCCAg                                                                                                                             >  1:243002/1‑149 (MQ=255)
                       cTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCAAAATTCAGGCAGCATATCGGTCCAg                                                                                                                             >  1:248123/1‑149 (MQ=255)
                            tACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTc                                                                                                                        <  4:259763/149‑1 (MQ=255)
                              cTGCACAGCGATGTTCAGACCATTTACCAGACCGAAGAATTCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGa                                                                                                                      <  6:262593/149‑1 (MQ=255)
                                  aCCGCGAAGTTCAGACCATTTACAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCTGGCAGCATATCGGTCCAGCCTTCGATGTCg                                                                                                                 <  8:357419/149‑1 (MQ=255)
                                      cgAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTAc                                                                                                              <  6:78598/149‑1 (MQ=255)
                                       gAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGAAGAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACa                                                                                                             >  1:251449/1‑149 (MQ=255)
                                                aCCATTTACCTGGCCGAAGAAATCAGTATTACGAAAAGTAGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATAGGTATAATAGGAAACGCCAAATTCAGGCACAATATCCTTCCAGCCTTCGAAGACGTACATTaaccca                                                                                                     >  8:325424/1‑143 (MQ=255)
                                                aCCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCa                                                                                                     >  5:326867/1‑148 (MQ=255)
                                                   aTTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACAtt                                                                                                           <  6:176938/139‑1 (MQ=255)
                                                   aTTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACAtt                                                                                                           >  5:176938/1‑139 (MQ=255)
                                                       cccAGACCGATGATATGAGTATTACGATAAGTCGCGACGCCATTGGCTCGTCCAGTCTTAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCGTTCGATTTCGTACATTACGCCATAATTAc                                                                                              <  6:312977/147‑1 (MQ=255)
                                                       aCCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTAc                                                                                              <  1:88371/148‑1 (MQ=255)
                                                        ccAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTCGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCTTAATTACGt                                                                                            >  2:67934/1‑149 (MQ=255)
                                                          acaCCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTcc                                                                                          <  5:317249/147‑1 (MQ=255)
                                                              cGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCCGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTaa                                                                                      >  8:19523/1‑149 (MQ=255)
                                                                   aaaTCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATc                                                                                    >  7:172497/1‑146 (MQ=255)
                                                                   aaaTCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATc                                                                                    <  8:172497/146‑1 (MQ=255)
                                                                     aTCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAg                                                                               <  5:286454/149‑1 (MQ=255)
                                                                      tCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAg                                                                               >  1:359882/1‑148 (MQ=255)
                                                                         gTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAAcc                                                                           >  6:357978/1‑149 (MQ=255)
                                                                                gATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCAAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCt                                                                    <  2:243002/149‑1 (MQ=255)
                                                                                gATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCAAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCt                                                                    <  2:248123/149‑1 (MQ=255)
                                                                                      tCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTg                                                                   >  7:319761/1‑144 (MQ=255)
                                                                                      tCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTg                                                                   <  8:319761/144‑1 (MQ=255)
                                                                                      tCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAAt                                                              >  6:302870/1‑149 (MQ=255)
                                                                                          gACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCa                                                          <  6:326867/149‑1 (MQ=255)
                                                                                          gACGCCATTGGCTCGACCAGTCATAAAGTTGGCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCa                                                          <  4:76503/149‑1 (MQ=255)
                                                                                          gACGCCATTGGCTCGACCAGTCATAAAGTCGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCTTAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCa                                                          <  1:67934/149‑1 (MQ=255)
                                                                                                      tCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGc                                                        <  1:88842/139‑1 (MQ=255)
                                                                                                      tCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGc                                                        >  2:88842/1‑139 (MQ=255)
                                                                                                        gACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATg                                              <  2:26386/147‑1 (MQ=255)
                                                                                                            aGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGAc                                        <  1:222633/149‑1 (MQ=255)
                                                                                                              tCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTaa                                                                       >  6:346996/1‑116 (MQ=255)
                                                                                                              tCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACgg                                      <  4:38235/149‑1 (MQ=255)
                                                                                                                aTAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGTAAATGCCAGACGGGt                                    >  2:177124/1‑149 (MQ=255)
                                                                                                                      ttGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCAt                                >  8:43303/1‑147 (MQ=255)
                                                                                                                           tGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCaa                                                 <  6:136535/125‑1 (MQ=255)
                                                                                                                           tGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCaa                                                 >  5:136535/1‑125 (MQ=255)
                                                                                                                                  gTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGttttt                     <  6:123435/146‑1 (MQ=255)
                                                                                                                                  gTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGttttt                     >  5:123435/1‑146 (MQ=255)
                                                                                                                                  gTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGt                         <  3:38471/142‑1 (MQ=255)
                                                                                                                                  gTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGt                         >  4:38471/1‑142 (MQ=255)
                                                                                                                                  gTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTg                   >  1:1379/1‑148 (MQ=255)
                                                                                                                                   tataAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTg                   >  7:261313/1‑147 (MQ=255)
                                                                                                                                   tataAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTGaa                 <  2:359882/149‑1 (MQ=255)
                                                                                                                                     taATAGTCTCCGCCAAATTCATGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACTTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCAt                                <  6:80680/132‑1 (MQ=255)
                                                                                                                                     taAGCGTCAACGCCAAATTCAGGCAGCATATGGGTCCAGCCTTCGTTGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGtt                        <  3:292996/140‑1 (MQ=255)
                                                                                                                                     taAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAATGAACCGAAAACTGCAAATTTCAGCCCGGCAAATGCAAGACGGGACCATGACTGGtt                        >  4:292996/1‑140 (MQ=255)
                                                                                                                                     taAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCAt                                >  5:80680/1‑132 (MQ=255)
                                                                                                                                     taAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTGAAGa               <  2:1379/149‑1 (MQ=255)
                                                                                                                                           tCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTGAAGATTCAgt         >  6:341076/1‑149 (MQ=255)
                                                                                                                                           tCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTGAAGATTCAgt         >  2:148146/1‑149 (MQ=255)
                                                                                                                                             aCCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTGAAGATTCAgtgt       >  2:41850/1‑149 (MQ=255)
                                                                                                                                                 ccAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTGTTTGAAGATTCAGTgttgtt   >  6:327835/1‑149 (MQ=255)
                                                                                                                                                  cAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTGAAGATTCAGTGTTGttt  <  5:302870/149‑1 (MQ=255)
                                                                                                                                                  cAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTGAAGATTCAGTGTTGttt  <  4:29572/149‑1 (MQ=255)
                                                                                                                                                  cAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTGAAGATTCAGTGTTGttt  >  7:264237/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATTACTGGCTCCTTCGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCATTGGCTCGACCAGTCATAAAGTTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTGAAGATTCAGTGTTGTTT  >  NC_000913/1436351‑1436645

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