Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,466,313 A→G L138P (CTT→CCT)  bfr ← bacterioferritin, iron storage and detoxification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,466,3130AG100.0% 258.9 / NA 72L138P (CTT→CCT) bfrbacterioferritin, iron storage and detoxification protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (38/34);  total (38/34)

AGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCA  >  NC_000913/3466179‑3466447
                                                                                                                                      |                                                                                                                                           
agcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTc                                                                                                                                         <  1:351885/139‑1 (MQ=255)
  caTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCg                                                                                                                                       >  1:115342/1‑139 (MQ=255)
       tAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTcc                                                                                                                                  >  2:163562/1‑139 (MQ=255)
              ggggTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑Gt                                                                                                                      >  1:812470/1‑139 (MQ=255)
                ggTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCg                                                                                                                    >  1:1187618/1‑139 (MQ=255)
                  tATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTAAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGat                                                                                                                  <  2:812470/139‑1 (MQ=255)
                   aTCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGata                                                                                                                 <  2:98316/139‑1 (MQ=255)
                   aTCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGata                                                                                                                 <  2:1187618/139‑1 (MQ=255)
                     cTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGata                                                                                                                 <  1:1355404/137‑1 (MQ=255)
                     cTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGata                                                                                                                 >  2:1355404/1‑137 (MQ=255)
                     cTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATg                                                                                                               <  2:115342/139‑1 (MQ=255)
                     cTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATg                                                                                                               >  2:1406523/1‑139 (MQ=255)
                              aTCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑Gt                                                                                                                      <  1:1375159/123‑1 (MQ=255)
                              aTCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑Gt                                                                                                                      >  2:1375159/1‑123 (MQ=255)
                              aTCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCttcttc                                                                                                      >  2:354303/1‑139 (MQ=255)
                                      ttaACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTc                                                                                                                                         >  1:724269/1‑101 (MQ=255)
                                      ttaACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTc                                                                                                                                         <  2:724269/101‑1 (MQ=255)
                                       taACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCaa                                                                                             >  1:242156/1‑139 (MQ=255)
                                           cGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑Gt                                                                                                                      >  1:1612451/1‑110 (MQ=255)
                                           cGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑Gt                                                                                                                      <  2:1612451/110‑1 (MQ=255)
                                           cGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAAtt                                                                                         >  1:1167590/1‑139 (MQ=255)
                                           cGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAAtt                                                                                         >  2:752044/1‑139 (MQ=255)
                                            gggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAAttt                                                                                        <  2:1167590/139‑1 (MQ=255)
                                             ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCGCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTtcat                                                                                                    <  1:1067476/126‑1 (MQ=255)
                                             ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTtcat                                                                                                    <  2:650201/126‑1 (MQ=255)
                                             ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTc                                                                                       >  1:1351476/1‑139 (MQ=255)
                                             ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGAACAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTtcat                                                                                                    >  2:1067476/1‑126 (MQ=255)
                                             ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCATCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGGTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTtcat                                                                                                    >  1:650201/1‑126 (MQ=255)
                                               aGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTa                                                                                     <  1:1406523/139‑1 (MQ=255)
                                               aGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTa                                                                                     >  2:1299555/1‑139 (MQ=255)
                                                    tctcCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGACAGGTCGATATGGCCTTCTTCATCACGCAAAATTTCTa                                                                                     >  2:1319294/1‑139 (MQ=255)
                                                    tctcCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTatcatc                                                                                <  1:1299555/139‑1 (MQ=255)
                                                             cGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTTCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTc                                                                                                                                         <  2:1053095/78‑1 (MQ=255)
                                                             cGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTc                                                                                                                                         >  1:1053095/1‑78 (MQ=255)
                                                                  cGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGata                                                                                                                 <  1:439279/92‑1 (MQ=255)
                                                                  cGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGata                                                                                                                 >  2:439279/1‑92 (MQ=255)
                                                                   ggCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACg                                                                 <  1:814705/139‑1 (MQ=255)
                                                                    gCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccag                                                                                                                                >  1:429424/1‑80 (MQ=255)
                                                                    gCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccag                                                                                                                                <  2:429424/80‑1 (MQ=255)
                                                                    gCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGt                                                                >  2:355002/1‑139 (MQ=255)
                                                                      tCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTaa                                                              >  2:357708/1‑139 (MQ=255)
                                                                               ttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACg                                                                                                <  2:672161/96‑1 (MQ=255)
                                                                               ttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACg                                                                                                >  1:672161/1‑96 (MQ=255)
                                                                                tcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGaa                                                       <  2:1214497/136‑1 (MQ=255)
                                                                                tcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGaa                                                       >  1:1214497/1‑136 (MQ=255)
                                                                                          gtgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATa                                          <  2:1351476/139‑1 (MQ=255)
                                                                                           tgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagc                                                                                                                               >  2:194636/1‑58 (MQ=255)
                                                                                           tgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagc                                                                                                                               <  1:194636/58‑1 (MQ=255)
                                                                                           tgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATaa                                         >  1:192625/1‑139 (MQ=255)
                                                                                                gCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGaa                                                       <  1:843794/120‑1 (MQ=255)
                                                                                                gCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGaa                                                       >  2:843794/1‑120 (MQ=255)
                                                                                                gCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAAt                                    >  2:806743/1‑139 (MQ=255)
                                                                                                     tAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTc                                                                                       >  2:463987/1‑83 (MQ=255)
                                                                                                     tAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTc                                                                                       <  1:463987/83‑1 (MQ=255)
                                                                                                     tAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCt                               <  2:227058/139‑1 (MQ=255)
                                                                                                                gCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAAttc                    >  1:47678/1‑139 (MQ=255)
                                                                                                                 cccATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAAttct                   <  1:357708/139‑1 (MQ=255)
                                                                                                                     tcttctGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTcgc               <  1:1393062/139‑1 (MQ=255)
                                                                                                                     tcttctGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTcgc               >  2:495348/1‑139 (MQ=255)
                                                                                                                     tcttctGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTcgc               >  2:1393062/1‑139 (MQ=255)
                                                                                                                      cttctGAATCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTcgcg              <  1:806743/139‑1 (MQ=255)
                                                                                                                            aaTCAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATcc        >  2:373441/1‑139 (MQ=255)
                                                                                                                               cAGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGc     >  1:110355/1‑139 (MQ=255)
                                                                                                                                aGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATc         <  1:785273/134‑1 (MQ=255)
                                                                                                                                aGATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATc         >  2:785273/1‑134 (MQ=255)
                                                                                                                                 gATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATc         >  1:1110600/1‑133 (MQ=255)
                                                                                                                                 gATCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATc         <  2:1110600/133‑1 (MQ=255)
                                                                                                                                  aTCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGc                                                                      >  2:255846/1‑71 (MQ=255)
                                                                                                                                  aTCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGc                                                                      <  1:255846/71‑1 (MQ=255)
                                                                                                                                  aTCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCa  <  1:373441/139‑1 (MQ=255)
                                                                                                                                  aTCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCa  <  2:1493643/139‑1 (MQ=255)
                                                                                                                                  aTCAGGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCa  <  1:354303/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                           
AGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCA‑‑‑‑‑GTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCA  >  NC_000913/3466179‑3466447

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: