Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,071,615 T→C 5.7% T48A (ACC→GCC)  yihS ← sulphoquinovose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,071,6150TC 5.7% 719.0 / 20.9 209T48A (ACC→GCC) yihSsulphoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base T (91/106);  new base C (5/7);  total (96/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 210 total aligned reads displayed.

AGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCATTTCATTCCTTTT  >  NC_000913/4071468‑4071763
                                                                                                                                                   |                                                                                                                                                    
aGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTg                                                                                                                                                     <  2:575030/149‑1 (MQ=255)
 ggCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGc                                                                                                                                                    >  2:424800/1‑149 (MQ=255)
  gCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGcc                                                                                                                                                   <  2:615246/149‑1 (MQ=255)
      aCCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTCATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCAtc                                                                                                                                               <  1:322410/149‑1 (MQ=255)
        cAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCAtctc                                                                                                                                             >  2:1523774/1‑149 (MQ=255)
         aGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATctct                                                                                                                                            <  1:56715/149‑1 (MQ=255)
          gccgccATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGCACCCAtctgtg                                                                                                                                           <  2:1578809/149‑4 (MQ=255)
                atatTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGa                                                                                                                                     <  4:863439/149‑1 (MQ=255)
                 tatTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGAt                                                                                                                                    <  2:1289512/149‑1 (MQ=255)
                 tatTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGAt                                                                                                                                    >  3:957653/1‑149 (MQ=255)
                 tatTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGAt                                                                                                                                    <  1:454204/149‑1 (MQ=255)
                     ttttATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCAtct                                                                                                                                              <  1:969105/135‑1 (MQ=255)
                     ttttATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGc                                                                                                                                <  3:1016319/149‑1 (MQ=255)
                     ttttATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGc                                                                                                                                >  2:1749947/1‑149 (MQ=255)
                          tCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCttt                                                                                                                           <  1:243820/149‑1 (MQ=255)
                          tCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGCACCCAtctg                                                                                                                                             <  2:81135/131‑2 (MQ=255)
                          tCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGCACCCAtctg                                                                                                                                             >  1:81135/1‑130 (MQ=255)
                            aCGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCtcttct                                                                                                                                         >  4:1425603/1‑133 (MQ=255)
                            aCGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCtcttct                                                                                                                                         <  3:1425603/133‑1 (MQ=255)
                            aCGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTAt                                                                                                                         >  2:1710338/1‑149 (MQ=255)
                             cGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTAtt                                                                                                                        <  4:163929/149‑1 (MQ=255)
                             cGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTAtt                                                                                                                        <  4:1068963/149‑1 (MQ=255)
                             cGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTAtt                                                                                                                        <  2:1081954/149‑1 (MQ=255)
                              gCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCAtc                                                                                                                                               <  4:1716590/125‑1 (MQ=255)
                              gCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCAtc                                                                                                                                               >  3:1716590/1‑125 (MQ=255)
                               cAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGc                                                                                                                      >  3:56534/1‑149 (MQ=255)
                               cAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGc                                                                                                                      <  3:1429305/149‑1 (MQ=255)
                               cAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGc                                                                                                                      <  4:189026/149‑1 (MQ=255)
                                     gcCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAAcc                                                                                                                >  4:787381/1‑149 (MQ=255)
                                      ccGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCa                                                                                                               <  2:640955/149‑1 (MQ=255)
                                      ccGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCa                                                                                                               >  2:378007/1‑149 (MQ=255)
                                      ccGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCa                                                                                                               <  4:56534/149‑1 (MQ=255)
                                      ccGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGCACCCAtctgtg                                                                                                                                           >  2:952169/1‑118 (MQ=255)
                                       cGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAg                                                                                                              <  1:1710338/149‑1 (MQ=255)
                                             tGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCaaaa                                                                                                        <  1:1519715/149‑1 (MQ=255)
                                                     aTACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGa                                                                                                                                     <  3:545086/112‑1 (MQ=255)
                                                     aTACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGa                                                                                                                                     >  4:545086/1‑112 (MQ=255)
                                                        ccGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGccc                                                                                                                                                  <  1:1830346/96‑1 (MQ=255)
                                                        ccGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGccc                                                                                                                                                  >  2:1830346/1‑96 (MQ=255)
                                                            gATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGc                                                                                               >  1:390916/1‑143 (MQ=255)
                                                             aTCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTc                                                                                                  <  2:1085270/139‑1 (MQ=255)
                                                             aTCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTc                                                                                                  >  1:1085270/1‑139 (MQ=255)
                                                              tCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGa                                                                                       <  1:1549904/149‑1 (MQ=255)
                                                                aaccaacGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCt                                                                                                                    <  2:1040218/118‑1 (MQ=255)
                                                                   caacGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAAtttt                                                                                  >  1:726466/1‑149 (MQ=255)
                                                                   caacGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAAtttt                                                                                  >  3:68095/1‑149 (MQ=255)
                                                                     acGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTa                                                                                <  1:1074374/149‑1 (MQ=255)
                                                                     acGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTa                                                                                <  1:378007/149‑1 (MQ=255)
                                                                      cGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTAc                                                                               <  4:957653/149‑1 (MQ=255)
                                                                        aGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTcttctt                                                                                                                                        <  4:737547/90‑1 (MQ=255)
                                                                        aGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTcttctt                                                                                                                                        >  3:1893427/1‑90 (MQ=255)
                                                                        aGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCa                                                                             <  4:1209319/149‑1 (MQ=255)
                                                                          tAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCaaa                                                                           >  1:923997/1‑149 (MQ=255)
                                                                            aGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACAAAAAt                                                                         >  3:1814647/1‑149 (MQ=255)
                                                                             gcgcCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATc                                                                        >  4:1029495/1‑149 (MQ=255)
                                                                                ccAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCaaa                                                                     <  4:1110803/149‑1 (MQ=255)
                                                                                   ggTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCaaaa                                                                                                        <  1:119748/111‑1 (MQ=255)
                                                                                   ggTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCaaaa                                                                                                        >  2:119748/1‑111 (MQ=255)
                                                                                    gTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCt                                                                                                                    >  2:1846338/1‑98 (MQ=255)
                                                                                    gTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATg                                                                 <  4:68095/149‑1 (MQ=255)
                                                                                        aCCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGc                                                                                                             >  3:761007/1‑101 (MQ=255)
                                                                                        aCCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCaaa                                                                           <  2:1515841/135‑1 (MQ=255)
                                                                                         cccATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGcggtc                                                            <  2:1668429/149‑1 (MQ=255)
                                                                                           cATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAAc                                                                                                       >  4:412102/1‑104 (MQ=255)
                                                                                           cATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAAc                                                                                                       <  3:412102/104‑1 (MQ=255)
                                                                                            aTCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGcggt                                                             >  4:105188/1‑145 (MQ=255)
                                                                                            aTCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGcggt                                                             <  3:105188/145‑1 (MQ=255)
                                                                                              cgccgcAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGttt                                                       <  3:999818/149‑1 (MQ=255)
                                                                                                ccgcAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACAAAAATCAAAGATGCGGTCGGTTTcc                                                     <  4:1814647/149‑1 (MQ=255)
                                                                                                   cagcaACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTAtt                                                                                                                        >  3:848635/1‑79 (MQ=255)
                                                                                                       aaCGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCa                                              >  2:854827/1‑149 (MQ=255)
                                                                                                        aCGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAg                                             >  3:792734/1‑149 (MQ=255)
                                                                                                          ggAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGcc                                           >  2:1455123/1‑149 (MQ=255)
                                                                                                           gAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAAttt                                                                                   >  1:1033823/1‑108 (MQ=255)
                                                                                                            aaTAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCgg                                                                                                <  2:303976/94‑1 (MQ=255)
                                                                                                                aaCGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTc                                                      >  3:1378924/1‑132 (MQ=255)
                                                                                                                aaCGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTc                                                      <  4:1378924/132‑1 (MQ=255)
                                                                                                                       aaCATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGtt                                                        <  2:556515/123‑1 (MQ=255)
                                                                                                                        aCATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCGGCCaaaa                                                                                                        <  3:1070210/74‑1 (MQ=255)
                                                                                                                           tACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAgg                          <  4:826954/149‑1 (MQ=255)
                                                                                                                               aGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCa                                                                                                     >  4:1686675/1‑70 (MQ=255)
                                                                                                                               aGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCa                                                                                                     <  3:1686675/70‑1 (MQ=255)
                                                                                                                                   gTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACg                                                                                         >  1:499159/1‑78 (MQ=255)
                                                                                                                                   gTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACg                                                                                         <  2:499159/78‑1 (MQ=255)
                                                                                                                                    tGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAAc                 >  3:812139/1‑149 (MQ=255)
                                                                                                                                      aTCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCa               >  3:1622183/1‑149 (MQ=255)
                                                                                                                                        gaagaGATGGGCACCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGt                                                  >  3:1286264/5‑112 (MQ=255)
                                                                                                                                        ccACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCAtt             >  1:59217/1‑149 (MQ=255)
                                                                                                                                         cacaGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTaaa                  <  4:1477907/143‑1 (MQ=255)
                                                                                                                                         cacaGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTaaa                  >  3:1477907/1‑143 (MQ=255)
                                                                                                                                         aagaGATGGGCACCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGaa                                                                                      <  3:714451/72‑1 (MQ=255)
                                                                                                                                         aagaGATGGGCACCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGaa                                                                                      >  4:714451/4‑75 (MQ=255)
                                                                                                                                          acaGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGtt                                                 >  1:1637759/1‑111 (MQ=255)
                                                                                                                                          acaGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTtgtg                                 >  4:1785048/1‑127 (MQ=255)
                                                                                                                                          acaGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTtgtg                                 <  3:1785048/127‑1 (MQ=255)
                                                                                                                                          acaGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTtg                                   >  1:1490098/1‑125 (MQ=255)
                                                                                                                                           gaGATGGGCACCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAAtttt                                                                                  <  1:126350/76‑1 (MQ=255)
                                                                                                                                            aGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCATTTCAt         >  1:1765278/1‑149 (MQ=255)
                                                                                                                                             gATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCATTTCAtt        <  3:1192815/149‑1 (MQ=255)
                                                                                                                                                  gTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCATTTCATTCCttt   <  2:726466/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGGAATAAACGTGCAACATACGAGCGGTGATCCACAGATGGGTGCCCATCTCTTCTTTGATTTGCCCTTTATTGCCTAACCAGCCAAAACCAGTCGGCACTACGGAATTTTTACCAAAATCAAAGATGCGGTCGGTTTCCTGTTCCAGCCAACGGTTGTGGCTTAGGGTGTTAAACCATTTCATTCCTTTT  >  NC_000913/4071468‑4071763

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