Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,246,005 G→A 5.5% Y255Y (TAC→TAT cirA ← colicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,246,0050GA 5.5% 237.3 / 2.7 73Y255Y (TAC→TATcirAcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter
Reads supporting (aligned to +/- strand):  ref base G (35/34);  new base A (2/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

GTTTTTGTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCGGTATCGGTGGTCGTTGAGTTTTGCGGGTCATCCTTTTCACGTTTTGCCA  >  NC_000913/2245861‑2246150
                                                                                                                                                |                                                                                                                                                 
ggttttGTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATc                                                                                                                                               <  5:138521/147‑1 (MQ=255)
gTTTTTGTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATc                                                                                                                                               <  7:60517/149‑1 (MQ=255)
 tttttGTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATcc                                                                                                                                              >  4:225288/1‑149 (MQ=255)
          tcGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTaa                                                                                                                                      >  6:212340/1‑148 (MQ=255)
            gACTTTCTCACCGTAGTATTTCAGTGCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTaaaa                                                                                                                                    <  6:52636/148‑1 (MQ=255)
                 tctcACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTg                                                                                                                               <  1:46668/148‑1 (MQ=255)
                     aCCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGAtttt                                                                                                                          >  5:139171/1‑149 (MQ=255)
                      ccGTAGTATTTCAGTTCGCTGGTGCCTTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGAGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGAtttt                                                                                                                          >  5:277881/1‑148 (MQ=255)
                        gtagtaTTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATc                                                                                                                                               <  6:79044/125‑1 (MQ=255)
                        gtagtaTTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATc                                                                                                                                               >  5:79044/1‑125 (MQ=255)
                         tagtaTTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGAtt                                                                                                                      >  8:282750/1‑149 (MQ=255)
                             aTTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGgt                                                                                                                  >  5:15676/1‑149 (MQ=255)
                              tttCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGgt                                                                                                                  <  5:376324/148‑1 (MQ=255)
                                tCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTc                                                                                                               >  1:374301/1‑149 (MQ=255)
                                tCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTc                                                                                                               >  3:243886/1‑149 (MQ=255)
                                      cGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGc                                                                                                          >  7:176013/1‑148 (MQ=255)
                                       gCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATGGTGATTTTGATTCGGTGTCCAGGCa                                                                                                         >  4:193350/1‑148 (MQ=255)
                                         tGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTc                                                                                                               >  2:399814/1‑140 (MQ=255)
                                         tGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTc                                                                                                               <  1:399814/140‑1 (MQ=255)
                                           gTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAAtt                                                                                                     <  5:372388/148‑1 (MQ=255)
                                              ccGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGa                                                                                                  <  6:169938/148‑1 (MQ=255)
                                                 tAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCCGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGGGATTTTGAGTCGGTGTCCAGGCAAATTCGACAtt                                                                                              >  5:300714/1‑149 (MQ=255)
                                                 tAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACAtt                                                                                              <  3:225288/149‑1 (MQ=255)
                                                  aaTCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCATa                                                                                                                                                 <  4:88266/97‑1 (MQ=255)
                                                  aaTCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCATa                                                                                                                                                 <  8:395019/97‑1 (MQ=255)
                                                  aaTCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCATa                                                                                                                                                 >  3:88266/1‑97 (MQ=255)
                                                  aaTCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCATa                                                                                                                                                 >  7:395019/1‑97 (MQ=255)
                                                   aTCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGAtttt                                                                                                                          >  3:48803/1‑119 (MQ=255)
                                                   aTCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGAtttt                                                                                                                          <  4:48803/119‑1 (MQ=255)
                                                       cAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTc                                                                                        >  6:218093/1‑149 (MQ=255)
                                                       cAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTc                                                                                        >  8:80753/1‑149 (MQ=255)
                                                        aaCGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTAGACATTGCCGTcg                                                                                       >  1:161879/1‑149 (MQ=255)
                                                         aCGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTcgc                                                                                      >  1:255731/1‑149 (MQ=255)
                                                             ccATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCt                                                                                  >  1:35444/1‑149 (MQ=255)
                                                             ccATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCt                                                                                  >  1:29471/1‑149 (MQ=255)
                                                                   tGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGaa                                                                            <  7:282750/149‑1 (MQ=255)
                                                                     gCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATc                                                                          <  2:66340/149‑1 (MQ=255)
                                                                       tGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCt                                                                        >  4:349542/1‑149 (MQ=255)
                                                                       tGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCt                                                                        >  6:387632/1‑149 (MQ=255)
                                                                       tGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCt                                                                        >  4:150993/1‑149 (MQ=255)
                                                                        gACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCtt                                                                       >  1:192263/1‑149 (MQ=255)
                                                                         aCGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCTTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCtt                                                                       <  6:139171/148‑1 (MQ=255)
                                                                             agtagtTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATa                                                                  <  5:218093/149‑1 (MQ=255)
                                                                               tagtTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATAcg                                                                >  7:157984/1‑149 (MQ=255)
                                                                                     tGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACAtt                                                                                              >  7:56453/1‑113 (MQ=255)
                                                                                     tGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACAtt                                                                                              <  8:56453/113‑1 (MQ=255)
                                                                                      ggCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGtt                                                         >  5:223001/1‑149 (MQ=255)
                                                                                             cAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCg                                                   <  3:150993/148‑1 (MQ=255)
                                                                                             cAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCg                                                   <  5:387632/148‑1 (MQ=255)
                                                                                              aGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTATGATTCGCTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGTAGAATCCTTCAATACGCGGCGTTTCTCCgg                                                  >  7:181468/1‑148 (MQ=255)
                                                                                                  ggTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCGGTATCg                                             >  1:102776/1‑149 (MQ=255)
                                                                                                  ggTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCGGTATCg                                             <  5:212340/149‑1 (MQ=255)
                                                                                                   gTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCGGTATCgg                                            >  4:362881/1‑149 (MQ=255)
                                                                                                    tttttGTCCAGCGAGTCGGAATCACGATCCTGACGTTCGAAACCGTATCCGGCAGTACAATCGTGATTATGATTCGTTGTCCAGGCAAATTCGACATTGCAGTCGCTTCTGGAGAATCCTTCAATACTCGTCGTTTATCCGGTATCggt                                           >  6:378754/1‑149 (MQ=255)
                                                                                                     ttttGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCGGTATCggtg                                          <  3:193350/149‑1 (MQ=255)
                                                                                                         gTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCGGTATCGGTGGTc                                       <  8:176013/148‑1 (MQ=255)
                                                                                                            cAGCGAGTCGGAATCACGATCCTGACTGTCGAAACCGTATCCGTGAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAATTTCGTCATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCTGCGTTTCTCCGGTATCGGTGGTCGtt                                    <  2:374301/148‑1 (MQ=255)
                                                                                                            cAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCGGTATCGGTGGTCGtt                                    <  4:243886/148‑1 (MQ=255)
                                                                                                                   tCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCGGTATCGGTGGTCGTTGAGTTTTg                            <  5:66610/149‑1 (MQ=255)
                                                                                                                     ggAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCGGTATCGGTGGTCGTTGAGTTTTGc                           <  2:192263/148‑1 (MQ=255)
                                                                                                                     ggAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCGGTATCGGTGGTCGTTGAGTTTTGc                           <  6:15676/148‑1 (MQ=255)
                                                                                                                      taaTCACGATCCTGACGGTCGAAAACGTATCAGGCAGTAAAATCTTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGTCGTTTCTCCTGTATCGGTGGTCGTTGAGTTTTGCgg                         <  6:277881/148‑1 (MQ=255)
                                                                                                                       aaTCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCGGTATCGGTGGTCGTTGAGTTTTGCggg                        <  2:102776/149‑1 (MQ=255)
                                                                                                                         tCACGATCCTGACGGTCGATACCGTATCCGGCAGTAAAATCGTGATTTTGAGTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCGGTATCGGTGGTCGTTGAGTTTTGCGGGt                       >  1:182233/1‑148 (MQ=255)
                                                                                                                         tCACCATCCGCACTGTCGAACCCGTATCCGCCAGTAAACTCGTGCTTTTGATTCGTAGTCCAGGCAAATTCGACATAGCCGTCGCGGCTGGAGAATCCTTCAATCCGCGTCGTTTCTCCGGTATCGGTGGTCGTTGAGTTTTGCggg                        <  4:109839/147‑1 (MQ=255)
                                                                                                                           aCGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGaa                                                                            <  5:28490/93‑1 (MQ=255)
                                                                                                                           aCGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGaa                                                                            >  6:28490/1‑93 (MQ=255)
                                                                                                                           aCGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCGGTATCGGTGGTCGTTGAGTTTTGCGGGTCa                     <  3:349542/148‑1 (MQ=255)
                                                                                                                               tCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCGGTATCggtgg                                         <  1:119241/124‑1 (MQ=255)
                                                                                                                               tCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCGGTATCggtgg                                         >  2:119241/1‑124 (MQ=255)
                                                                                                                                             aCCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCGGTATCGGTGGTCGTTGAGTTTTGCGGGTCATCCTTTTCACGTTTTGCCa  <  3:362881/149‑1 (MQ=255)
                                                                                                                                               cGTATCGGGCAGTAATATCGTGATTTTGATTCGGTGTCCAGGCATATTCGAGATTGCCGTCGCGGCTGGAGTATGCTTCAATACGCGGCGTTTCTCCGGTATCGGTGGTCGTTGAGTTTTGCGGGTCTTCCTTTTCa            <  8:250243/137‑1 (MQ=255)
                                                                                                                                               cGTATCCGGCAGTAAAATCTTGATTTTGATTCGGTGTCCACGCAAATTCGACATTGCCGTCTCGGCTGGAGAATCCTTCAATACGCGGCGTTTATCCGGTATCGGTGGTCGTTGAGTATTGCGGGTCATCCTTTTCa            >  7:250243/1‑137 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GTTTTTGTTCTCGACTTTCTCACCGTAGTATTTCAGTTCGCTGGTGCCGTAATCCCAACGCCCATTATGGCTGACGGAGTAGTTCTGGCGTTCCAGGCGGTTTTTGTCCAGCGAGTCGGAATCACGATCCTGACGGTCGAAACCGTATCCGGCAGTAAAATCGTGATTTTGATTCGGTGTCCAGGCAAATTCGACATTGCCGTCGCGGCTGGAGAATCCTTCAATACGCGGCGTTTCTCCGGTATCGGTGGTCGTTGAGTTTTGCGGGTCATCCTTTTCACGTTTTGCCA  >  NC_000913/2245861‑2246150

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: