Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,515,070 G→A D536D (GAC→GAT fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,515,0700GA98.5% 219.1 / NA 65D536D (GAC→GATfecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (30/34);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGC  >  NC_000913/4514934‑4515206
                                                                                                                                        |                                                                                                                                        
cgtcgtcgtAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGGACGGTGCCAAACAAGCGTTCAGTGTATGCATAAAGATTCCAGCTatca                                                                                                                                       >  2:288737/1‑140 (MQ=255)
 gtcgtcgtAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTatcag                                                                                                                                      <  2:238873/140‑1 (MQ=255)
 gtcgtcgtAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTatcag                                                                                                                                      <  2:182675/140‑1 (MQ=255)
 ctcgtcgtAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTatcag                                                                                                                                      <  2:1002639/139‑1 (MQ=255)
   cgtcgtAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTatcagtc                                                                                                                                    <  2:880571/140‑1 (MQ=255)
   cgtcgtAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTatcagtc                                                                                                                                    <  2:518834/140‑1 (MQ=255)
     tcgtAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTAtcagtcag                                                                                                                                  <  1:25232/140‑1 (MQ=255)
      cgtAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGTTCCAGCTATCAGTCAGGt                                                                                                                                <  2:149046/140‑1 (MQ=255)
         caggggTACCGAGTTCCCAGGTACGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGa                                                                                                                              <  1:465730/137‑1 (MQ=255)
         aGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCGGGTGa                                                                                                                              <  2:800120/140‑1 (MQ=255)
         aGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGa                                                                                                                              <  1:130361/140‑1 (MQ=255)
               tACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGc                                                                                                                        >  2:963880/1‑140 (MQ=255)
                       tCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATCAAGCTTCCAGCTATCAATCAGGCGATAGAGCACGTTctg                                                                                                                >  2:465730/1‑138 (MQ=255)
                                cgcgCTTTTTCCGGTTGAACATGGCCGCTTGTCACCGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAg                                                                                                       <  2:491566/140‑1 (MQ=255)
                                cgcgCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATGAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAg                                                                                                       >  2:313294/1‑140 (MQ=255)
                                cgcgCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAg                                                                                                       <  2:970336/140‑1 (MQ=255)
                                cgcgCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAg                                                                                                       >  1:226526/1‑140 (MQ=255)
                                cgcgCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAATCAGGTGATAGAGCACGTTCAACGCCGGAAg                                                                                                       >  2:379354/1‑140 (MQ=255)
                                  cgcTTTTTCCGGTTCCACATTGTCGCTTTCCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCg                                                                                                     <  1:274781/140‑1 (MQ=255)
                                  cgcTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAg                                                                                                       >  2:146932/1‑138 (MQ=255)
                                  cgcTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAg                                                                                                       <  1:146932/138‑1 (MQ=255)
                                   gcTTTTTCCGGTTCACCATTGCCGCGTTGTACAGCCTAGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCgg                                                                                                    <  1:288737/140‑1 (MQ=255)
                                     tttttCCGGTTCAACATCGCCGCTTTGCACATCCTCGCCTATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTg                                                                                                  <  2:295355/140‑1 (MQ=255)
                                             gTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCg                                                                                                     <  2:722336/129‑1 (MQ=255)
                                             gTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCg                                                                                                     >  1:722336/1‑129 (MQ=255)
                                                        ccGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACttcttc                                                                               >  2:780021/1‑140 (MQ=255)
                                                            tttGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTcgtgc                                                                           <  1:379354/140‑1 (MQ=255)
                                                             ttGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTcgtgcg                                                                          >  1:6439/1‑140 (MQ=255)
                                                                           cAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCg                                                            <  1:963880/140‑1 (MQ=255)
                                                                           cAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCg                                                            <  2:503780/140‑1 (MQ=255)
                                                                           cAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCg                                                            <  2:968956/140‑1 (MQ=255)
                                                                           cAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCg                                                            <  2:226526/140‑1 (MQ=255)
                                                                           cAATCTGGCTGTACTGTACGGTGCCGAACGAGCCATCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCg                                                            <  1:633879/140‑1 (MQ=255)
                                                                           cAATCTGGCTGTACTGCTCGGTGCCGAACGTCCCTTCAGTGTTTGCATGAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCg                                                            <  2:483730/140‑1 (MQ=255)
                                                                                tGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGGTAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCgttgtt                                                       <  2:917483/140‑1 (MQ=255)
                                                                                tGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCgttgtt                                                       <  2:987300/140‑1 (MQ=255)
                                                                                 ggCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTc                                                      <  1:905699/140‑1 (MQ=255)
                                                                                  gCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCt                                                     >  1:232628/1‑140 (MQ=255)
                                                                                     gtactgtacGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCTCTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGt                                                  <  2:660970/140‑1 (MQ=255)
                                                                                            acGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTc                                           >  1:606500/1‑140 (MQ=255)
                                                                                            acGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTc                                           >  1:70700/1‑140 (MQ=255)
                                                                                            acGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTc                                           >  2:69149/1‑140 (MQ=255)
                                                                                               gTGCCGAACGAGCCTCCAGTGTTTGTATATAGATTCCCGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGat                                        <  1:288960/140‑1 (MQ=255)
                                                                                                tGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTTTGATGGTGTTGTTATGGTATGACTCGata                                       >  1:452015/1‑140 (MQ=255)
                                                                                                tGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGTGTTGTTCTGGTATGACTCGata                                       >  1:306544/1‑140 (MQ=255)
                                                                                                tGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGGGTTGTTCTGGTATGACTCGata                                       >  1:336989/1‑140 (MQ=255)
                                                                                                tGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGata                                       >  1:74730/1‑140 (MQ=255)
                                                                                                tGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGata                                       >  1:537471/1‑140 (MQ=255)
                                                                                                tGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGata                                       >  2:940683/1‑140 (MQ=255)
                                                                                                tGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGa                                         <  2:189663/138‑1 (MQ=255)
                                                                                                tGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGa                                         >  1:189663/1‑138 (MQ=255)
                                                                                                 gCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGatat                                      <  2:10575/140‑1 (MQ=255)
                                                                                                  ccGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATg                                     >  1:903707/1‑140 (MQ=255)
                                                                                                   cGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGtt                                   >  2:870560/1‑140 (MQ=255)
                                                                                                     aaCGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTc                                  >  2:55893/1‑140 (MQ=255)
                                                                                                     aaCGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTc                                  >  2:161081/1‑140 (MQ=255)
                                                                                                     aaCGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTc                                  >  2:834349/1‑140 (MQ=255)
                                                                                                      aCGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGa                                <  1:870560/140‑1 (MQ=255)
                                                                                                       cGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGa                                <  2:777617/140‑1 (MQ=255)
                                                                                                          gCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACt                                            >  2:1020165/1‑125 (MQ=255)
                                                                                                          gCCTTCAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACt                                            <  1:1020165/125‑1 (MQ=255)
                                                                                                               cAGTGTTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATa                        <  2:70700/140‑1 (MQ=255)
                                                                                                                 gtgtTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATAcc                      <  2:795328/140‑1 (MQ=255)
                                                                                                                  tgtTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCATGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATa                        >  1:579148/1‑137 (MQ=255)
                                                                                                                  tgtTTGCATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCATGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATa                        <  2:579148/137‑1 (MQ=255)
                                                                                                                        cATAAAGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTg               >  1:776773/1‑140 (MQ=255)
                                                                                                                          tAAAGATTCCAGCTATCGGTCAGGTGATAGAGCACGTTCGACGCCGGAGGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGAt             <  2:306544/140‑1 (MQ=255)
                                                                                                                            aaGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATgg           >  1:598442/1‑140 (MQ=255)
                                                                                                                            aaGATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGGTGTTCTGGGATGGCTCGATATGTTCAAAACTCATACCCGGCGTGATgg           >  2:323206/1‑140 (MQ=255)
                                                                                                                              gATTCCAGCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTc         <  2:776773/140‑1 (MQ=255)
                                                                                                                                     gCTATCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGc  >  2:636670/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                        
CGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAGCCTTGCCAATCTGGCTGTACTGTACGGTGCCGAACGAGCCTTCAGTGTTTGCATAAAGATTCCAGCTGTCAGTCAGGTGATAGAGCACGTTCAACGCCGGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGC  >  NC_000913/4514934‑4515206

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