Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,550,518 T→C 5.2% N260S (AAC→AGC)  gntP ← fructuronate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,550,5180TC 5.2% 3200.6 / 96.7 970N260S (AAC→AGC) gntPfructuronate transporter
Reads supporting (aligned to +/- strand):  ref base T (492/427);  new base C (28/23);  total (520/450)
Fisher's exact test for biased strand distribution p-value = 8.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 1001 total aligned reads displayed.

GATGATCAGAATCACCATTGCAATACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAG  >  NC_000913/4550370‑4550666
                                                                                                                                                    |                                                                                                                                                    
 atgatCAGAATCACCATTGCAATACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGtt                                                                                                                                                     <  4:4004403/149‑1 (MQ=255)
    atCAGAATCACCATTGCAATACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAAc                                                                                                                                                  <  4:3519730/149‑1 (MQ=255)
    atCAGAATCACCATTGCAATACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAAc                                                                                                                                                  <  3:742045/149‑1 (MQ=255)
    atCAGAATCACCATTGCAATACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAAc                                                                                                                                                  >  2:1652229/1‑149 (MQ=255)
    atCAGAATCACCATTGCAATACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAAc                                                                                                                                                  <  2:1868469/149‑1 (MQ=255)
     tCAGAATCACCATTGCAATACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCTCCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGtt                                                                                                                                                     <  2:1133745/145‑1 (MQ=255)
      cAGAATCACCATTGCAATACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCa                                                                                                                                                <  1:6669562/149‑1 (MQ=255)
       aGAATCACCATTGCAATACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGCTAACCAc                                                                                                                                               <  3:1336621/149‑1 (MQ=255)
             aCCATTGCAATACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCca                                                                                                                                         >  3:3519068/1‑149 (MQ=255)
             aCCATTGCAATACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCca                                                                                                                                         <  3:4654725/149‑1 (MQ=255)
              ccATTGCAATACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCcag                                                                                                                                        >  3:5167055/1‑149 (MQ=255)
               cATTGCAATACTCTTCACCGCGCTTTCATAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCcagg                                                                                                                                       >  2:2316810/1‑149 (MQ=255)
                   gCAATACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAggg                                                                                                                                  >  3:66251/1‑149 (MQ=255)
                    cAATACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAggg                                                                                                                                  <  2:513781/149‑1 (MQ=255)
                      aTACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCcag                                                                                                                                        <  1:2809348/141‑1 (MQ=255)
                       tACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAggg                                                                                                                                  <  2:2731833/146‑1 (MQ=255)
                       tACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTa                                                                                                                               >  3:1622674/1‑149 (MQ=255)
                       tACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTa                                                                                                                               >  2:956099/1‑149 (MQ=255)
                          tctTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCt                                                                                                                            >  3:4359840/1‑149 (MQ=255)
                           ctTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCtt                                                                                                                           >  1:1017571/1‑149 (MQ=255)
                             tCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCtttt                                                                                                                         <  1:4010474/149‑1 (MQ=255)
                              cACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTa                                                                                                                        <  2:2111405/149‑1 (MQ=255)
                              cACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTa                                                                                                                        <  3:6075495/149‑1 (MQ=255)
                               aCCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTaa                                                                                                                                                   >  1:4288274/1‑121 (MQ=255)
                               aCCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTAc                                                                                                                       >  2:3871262/1‑149 (MQ=255)
                               aCCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTAc                                                                                                                       <  4:716961/149‑1 (MQ=255)
                                      tttCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCt                                                                                                                            >  1:5140051/1‑137 (MQ=255)
                                       ttCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGa                                                                                                               <  1:5957159/149‑1 (MQ=255)
                                         cAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGCTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGata                                                                                                             <  1:1882356/149‑1 (MQ=255)
                                          aaaaGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGatat                                                                                                            >  4:2986488/1‑149 (MQ=255)
                                          aaaaGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGCTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGatat                                                                                                            >  2:6454258/1‑149 (MQ=255)
                                           aaaGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATAtt                                                                                                           >  3:181498/1‑149 (MQ=255)
                                           aaaGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGCTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATAtt                                                                                                           >  4:5321455/1‑149 (MQ=255)
                                            aaGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTg                                                                                                          <  1:4945407/149‑1 (MQ=255)
                                              gCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGgc                                                                                                        >  1:7200695‑M2/1‑147 (MQ=255)
                                                  tcatcaCCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGgcgatg                                                                                                    >  3:974073‑M2/1‑143 (MQ=255)
                                                    atcaCCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGtt                                                                                                                                                     >  4:131649/1‑98 (MQ=255)
                                                     tcaCCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGt                                                                                                                                >  1:4593708/1‑118 (MQ=255)
                                                     tcaCCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGgcgatggtg                                                                                                 <  4:6268744‑M2/149‑10 (MQ=255)
                                                     tcaCCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGgcgatggtg                                                                                                 <  1:6302487‑M2/149‑10 (MQ=255)
                                                        cccACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGCTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTa                                                                                              >  4:170863/1‑149 (MQ=255)
                                                         ccACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTag                                                                                             <  3:3212497/149‑1 (MQ=255)
                                                         ccACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTAGCTTTTACCAGCCAGATATTGGCGATGGTGGTag                                                                                             >  4:4714366/1‑149 (MQ=255)
                                                         ccACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGCTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTag                                                                                             >  1:5344805/1‑149 (MQ=255)
                                                          cACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTaga                                                                                            <  3:3886202/149‑1 (MQ=255)
                                                           aCTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTagag                                                                                           >  4:236735/1‑149 (MQ=255)
                                                            cTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAggg                                                                                                                                  >  2:462512/1‑109 (MQ=255)
                                                                  gTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATaa                                                                                    >  3:5844911/1‑149 (MQ=255)
                                                                  gTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATaa                                                                                    <  1:6907940/149‑1 (MQ=255)
                                                                     aTGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGa                                                                                 >  1:2588595/1‑149 (MQ=255)
                                                                     aTGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGa                                                                                 >  1:3041610/1‑149 (MQ=255)
                                                                     aTGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGCAGAGATCATAAtaa                                                                                 <  4:5885707/149‑3 (MQ=255)
                                                                      tGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGAt                                                                                <  4:1914713/149‑1 (MQ=255)
                                                                      tGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGAt                                                                                >  4:3309657/1‑149 (MQ=255)
                                                                       gACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTAccagcca                                                                                                                 <  3:5674142/115‑1 (MQ=255)
                                                                       gACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATc                                                                               >  1:1388984/1‑149 (MQ=255)
                                                                           aCGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCg                                                                           <  2:6709309/149‑1 (MQ=255)
                                                                              tGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGa                                                                        <  3:2965817/149‑1 (MQ=255)
                                                                                     ccAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCTATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGgcgatgg                                                                                                   <  3:4093538‑M2/115‑8 (MQ=255)
                                                                                     ccAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCgg                                                                 <  2:2512746/149‑1 (MQ=255)
                                                                                      ctaatagTTCGAATGCCACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGAAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATAATAATGATCGCTGGGATCAGCGGc                                                                >  4:4178900/6‑149 (MQ=255)
                                                                                      cAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGc                                                                >  3:2010869/1‑149 (MQ=255)
                                                                                          gagTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCa                                                            >  2:3255209/1‑149 (MQ=255)
                                                                                             tACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGaa                                                         <  3:6844282/149‑1 (MQ=255)
                                                                                             tACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGaa                                                         >  2:1391334/1‑149 (MQ=255)
                                                                                               cGAATGCAACCACCATCGCAGTAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTTCCAGGCAGGGGAATCTTTTACCAGCCAGATATTGGCGATCGTGATAGAGACCATAATGATCGCCGTGACCAGCGGCACTAGAATc                                                       >  2:4355249/1‑149 (MQ=255)
                                                                                               cGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATc                                                       <  3:3664650/149‑1 (MQ=255)
                                                                                                   tGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGaaa                                                   <  4:5059005/149‑1 (MQ=255)
                                                                                                     cAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAggg                                                                                                                                  >  2:5807678/1‑68 (MQ=255)
                                                                                                      aaCCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGc                                                >  1:1650192/1‑149 (MQ=255)
                                                                                                       accaccATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGaa                                                         <  3:1079787/139‑1 (MQ=255)
                                                                                                           ccaTCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGaaa                                           <  1:5644490/149‑1 (MQ=255)
                                                                                                           ccaTCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGaaa                                           >  2:535043/1‑149 (MQ=255)
                                                                                                           ccaTCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGaaa                                           >  3:6326107/1‑149 (MQ=255)
                                                                                                            caTCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAag                                          >  1:2119338/1‑149 (MQ=255)
                                                                                                             aTCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAaga                                         <  3:2517818/149‑1 (MQ=255)
                                                                                                              tCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAagag                                        <  1:3369263/149‑1 (MQ=255)
                                                                                                              tCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAagag                                        <  2:3520959/149‑1 (MQ=255)
                                                                                                                  aaTAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCa                                    >  4:6060134/1‑149 (MQ=255)
                                                                                                                       aCATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGatta                               >  4:6462052/1‑149 (MQ=255)
                                                                                                                         aTTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGattatt                             >  3:2224600/1‑149 (MQ=255)
                                                                                                                         aTTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGattatt                             <  2:4043223/149‑1 (MQ=255)
                                                                                                                          ttGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTc                            >  1:3791070/1‑149 (MQ=255)
                                                                                                                              aaTCGGCGAGTAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTAGCTGTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCACCGGCACCAGAATCGATACGCCGAAAGAGGGCAGATTATTCatat                        >  4:3712266/1‑149 (MQ=255)
                                                                                                                               aTCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATc                       >  4:2613783/1‑149 (MQ=255)
                                                                                                                                tCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCt                      >  4:1112873/1‑149 (MQ=255)
                                                                                                                                tCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCt                      >  2:5556918/1‑149 (MQ=255)
                                                                                                                                 cGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTa                     >  4:5667134/1‑149 (MQ=255)
                                                                                                                                   gCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTAcc                   >  2:6043415/1‑149 (MQ=255)
                                                                                                                                     gAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCgg                 <  2:7062508/149‑1 (MQ=255)
                                                                                                                                       ggAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCAAGATATGGGCGATGGTGGGAGAGAGCATAATGATCGACGTGATCAGCGGCACCCGAATCCGAACGGCGAAAGAGTGGAGATTATTCATATCTACCGGtt               <  4:6842417/149‑1 (MQ=255)
                                                                                                                                        gAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTAcc                   >  1:233687/1‑144 (MQ=255)
                                                                                                                                        gAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTg              >  1:876543/1‑149 (MQ=255)
                                                                                                                                         aaCCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGa             <  1:6487224/149‑1 (MQ=255)
                                                                                                                                          aCCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGAt            <  2:4882616/149‑1 (MQ=255)
                                                                                                                                          aCCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGAt            <  2:2958842/149‑1 (MQ=255)
                                                                                                                                            cGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCCGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCt          >  3:2690977/1‑149 (MQ=255)
                                                                                                                                             gATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTg         >  2:764693/1‑149 (MQ=255)
                                                                                                                                              aTAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGc        >  3:5362636/1‑149 (MQ=255)
                                                                                                                                                  aGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTc    >  4:2187501/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATGATCAGAATCACCATTGCAATACTCTTCACCGCGCTTTCAAAAGCGTTCATCACCCACTGCATGTCATGACCACGTGCGGTGCCAAAGAGTACGAATGCAACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAG  >  NC_000913/4550370‑4550666

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