Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 894,952 T→A 7.5% intergenic (+56/‑39) potF → / → potG putrescine ABC transporter periplasmic binding protein/putrescine ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913894,9520TA 7.5% 279.6 / 10.3 80intergenic (+56/‑39)potF/potGputrescine ABC transporter periplasmic binding protein/putrescine ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (37/37);  new base A (3/3);  total (40/40)
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.51e-01

GTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGAT  >  NC_000913/894808‑895098
                                                                                                                                                |                                                                                                                                                  
gTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATgg                                                                                                                                                <  4:893172/149‑1 (MQ=255)
    gAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTgcg                                                                                                                                            >  4:524134/1‑149 (MQ=255)
          gTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGca                                                                                                                                      <  1:1454574/149‑1 (MQ=255)
          gTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATAATGGTgcgaa                                                                                                                                          >  4:449816/1‑143 (MQ=255)
          gTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATAATGGTgcgaa                                                                                                                                          <  3:449816/145‑3 (MQ=255)
              aCTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGcacaca                                                                                                                                  >  2:304506/1‑149 (MQ=255)
              aCTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGcacaca                                                                                                                                  <  1:863028/149‑1 (MQ=255)
               ctctGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACAt                                                                                                                                 <  3:524134/149‑1 (MQ=255)
                tctGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACAtt                                                                                                                                <  1:304506/149‑1 (MQ=255)
                tctGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATAATGGTgcgaa                                                                                                                                          <  4:311350/139‑3 (MQ=255)
                tctGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATAATGGTgcgaa                                                                                                                                          >  3:311350/1‑137 (MQ=255)
                    aaaGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAAt                                                                                                                            <  2:52525/149‑1 (MQ=255)
                    aaaGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAAt                                                                                                                            <  2:249984/149‑1 (MQ=255)
                             gATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGagag                                                                                                                   >  4:1171042/1‑149 (MQ=255)
                              aTCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGagc                                                                                                                  >  1:89578/1‑149 (MQ=255)
                                ccGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGagcag                                                                                                                <  1:1012791/149‑1 (MQ=255)
                                  gAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGc                                                                                                                          <  1:127684/137‑1 (MQ=255)
                                  gAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGc                                                                                                                          >  2:127684/1‑137 (MQ=255)
                                  gAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGcc                                                                                                              >  2:49739/1‑149 (MQ=255)
                                     aaTCGACCGTGTGCGCACCTGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGagag                                                                                                                   >  4:1244055/1‑141 (MQ=255)
                                     aaTCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGagag                                                                                                                   <  3:1244055/141‑1 (MQ=255)
                                       tCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGaa                                                                                                         >  4:256977/1‑149 (MQ=255)
                                          aCCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGca                                                                                                                                      <  1:1390942/117‑1 (MQ=255)
                                          aCCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGca                                                                                                                                      >  2:1390942/1‑117 (MQ=255)
                                          aCCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGa                                                                                                      >  2:807839/1‑149 (MQ=255)
                                           ccGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCg                                                                                                                        >  1:1140506/1‑130 (MQ=255)
                                           ccGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCg                                                                                                                        <  2:1140506/130‑1 (MQ=255)
                                           ccGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGAc                                                                                                     <  3:256977/149‑1 (MQ=255)
                                              tgtgCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCt                                                                                                  >  2:1200848/1‑149 (MQ=255)
                                                   gcACCCGCGCGTGGACCAAAGTGAAGTGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGga                                                                                                                      >  2:684732/1‑124 (MQ=255)
                                                   gcACCCGCGCGTGGACCAAAGTGAAGTGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGga                                                                                                                      <  1:684732/124‑1 (MQ=255)
                                                   gcACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATccc                                                                                             >  3:1353761/1‑149 (MQ=255)
                                                       ccGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGc                                                                                         >  2:200666/1‑149 (MQ=255)
                                                          cgcgTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCg                                                                                      <  4:1353761/149‑1 (MQ=255)
                                                             gTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAg                                                                                   <  1:49739/149‑1 (MQ=255)
                                                               ggACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGc                                                                                 <  2:89578/149‑1 (MQ=255)
                                                                gACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCg                                                                                >  4:625091/1‑149 (MQ=255)
                                                                    aaaGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCAGTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGaaaa                                                                            <  3:1171042/149‑1 (MQ=255)
                                                                        tGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGga                                                                                                                      <  4:345250/103‑1 (MQ=255)
                                                                        tGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGga                                                                                                                      >  3:345250/1‑103 (MQ=255)
                                                                         gAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGaa                                                                                                         <  4:1134056/115‑1 (MQ=255)
                                                                         gAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGaa                                                                                                         >  3:1134056/1‑115 (MQ=255)
                                                                          aaGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTa                                                                      >  3:965870/1‑149 (MQ=255)
                                                                               cGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGgcg                                                                 >  4:813456/1‑149 (MQ=255)
                                                                                ggAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGcacaca                                                                                                                                  >  3:735000/1‑83 (MQ=255)
                                                                                ggAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGcacaca                                                                                                                                  <  4:735000/83‑1 (MQ=255)
                                                                                      taatCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTaa                                                                     >  2:1447908/1‑138 (MQ=255)
                                                                                      taatCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTaa                                                                     <  1:1447908/138‑1 (MQ=255)
                                                                                      taatCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAcgc                                                          >  2:1376802/1‑149 (MQ=255)
                                                                                         tCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTa                                                                                                                                                   <  1:1126309/57‑1 (MQ=255)
                                                                                         tCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTa                                                                                                                                                   >  2:1126309/1‑57 (MQ=255)
                                                                                          ccGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCt                                                      >  4:115761/1‑149 (MQ=255)
                                                                                              aGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTa                                                  <  3:813456/149‑1 (MQ=255)
                                                                                              aGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTa                                                  <  3:625091/149‑1 (MQ=255)
                                                                                                  gTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGcacaca                                                                                                                                  <  2:836010/65‑1 (MQ=255)
                                                                                                  gTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGcacaca                                                                                                                                  >  1:836010/1‑65 (MQ=255)
                                                                                                  gTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGaaa                                              >  1:350459/1‑149 (MQ=255)
                                                                                                     gATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGgcgc                                                                <  1:569185/128‑1 (MQ=255)
                                                                                                     gATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGgcgc                                                                >  2:569185/1‑128 (MQ=255)
                                                                                                         ccGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCaa                                       >  3:1139215/1‑149 (MQ=255)
                                                                                                              ggggCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGa                                  <  1:1200848/149‑1 (MQ=255)
                                                                                                               gggCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGAc                                 <  1:807839/149‑1 (MQ=255)
                                                                                                                  cgcACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCaa                              >  3:1046608/1‑149 (MQ=255)
                                                                                                                    caccacACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAAt                            <  4:1046608/149‑1 (MQ=255)
                                                                                                                     accacACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATc                           <  1:92404/149‑1 (MQ=255)
                                                                                                                                 ggCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCa               >  1:1110030/1‑149 (MQ=255)
                                                                                                                                    aaTTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACa            >  1:820514/1‑149 (MQ=255)
                                                                                                                                      ttCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAcgc                                                          <  3:239439/101‑1 (MQ=255)
                                                                                                                                      ttCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTAc                       >  3:883589/1‑136 (MQ=255)
                                                                                                                                      ttCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTAc                       <  4:883589/136‑1 (MQ=255)
                                                                                                                                      ttCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATg          >  3:445721/1‑149 (MQ=255)
                                                                                                                                      ttCGAAACATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAcgc                                                          >  4:239439/1‑101 (MQ=255)
                                                                                                                                       gcGCACCATAATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGaaa                                              >  2:610578/2‑112 (MQ=255)
                                                                                                                                       gcGCACCATAATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGaaa                                              <  1:610578/111‑1 (MQ=255)
                                                                                                                                          cACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGt      <  4:1139215/149‑1 (MQ=255)
                                                                                                                                           aCCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAcgccgc                                                       <  4:1037905/99‑1 (MQ=255)
                                                                                                                                           aCCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAcgccgc                                                       >  3:1037905/1‑99 (MQ=255)
                                                                                                                                           aCCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATcc                                           <  2:1002783/111‑1 (MQ=255)
                                                                                                                                           aCCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATcc                                           >  1:1002783/1‑111 (MQ=255)
                                                                                                                                              aTTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGgat  <  2:350459/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GTGCGAAGCTGTTCACTCTGAAAGTGCAGGATCCGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGAT  >  NC_000913/894808‑895098

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