Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,245,148 T→A 8.9% D541V (GAT→GTT)  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,245,1480TA 8.9% 185.2 / 6.4 64D541V (GAT→GTT) cirAcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base A (4/2);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAG  >  NC_000913/2245000‑2245295
                                                                                                                                                    |                                                                                                                                                   
aTAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCCGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCAt                                                                                                                                                     >  1:1035499/1‑149 (MQ=255)
        aCCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAGCGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAAc                                                                                                                                             <  2:1226293/149‑1 (MQ=255)
         ccAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGCCCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACg                                                                                                                                            >  1:467928/1‑149 (MQ=255)
            gTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGaa                                                                                                                                         <  2:556528/149‑1 (MQ=255)
             tCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAAt                                                                                                                                        >  1:1209158/1‑149 (MQ=255)
               ttCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAAttt                                                                                                                                      >  2:429493/1‑149 (MQ=255)
               ttCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAAttt                                                                                                                                      >  1:495305/1‑149 (MQ=255)
               ttCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAAttt                                                                                                                                      >  2:1545989/1‑149 (MQ=255)
                 ccagcgccagcgGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCGTTCAACGATGAATGGaaactgt                                                                                                                                    <  1:926539/149‑6 (MQ=255)
                      gccagcgGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTc                                                                                                                               <  2:549699/149‑1 (MQ=255)
                       ccagcgGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCAGTCATCGGTGAACGGAATTTTCAGTTCg                                                                                                                              >  2:1260651/1‑149 (MQ=255)
                         agcgGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGATTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGt                                                                                                                            <  1:1545522/149‑1 (MQ=255)
                              tttCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCa                                                                                                                       >  1:1406517/1‑149 (MQ=255)
                                tccagtccagCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACg                                                                                                                     >  1:1495187/1‑149 (MQ=255)
                                  cagtccagCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGcc                                                                                                                   >  1:12663/1‑149 (MQ=255)
                                        agCGTTCCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCAACGTTGAACGGAATTTTCAGTTCGTTTTCCACGCCCTGatt                                                                                                             >  1:515615/1‑147 (MQ=255)
                                        agCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAAt                                                                                                             <  1:185086/149‑1 (MQ=255)
                                            tACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGCACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGa                                                                                                         >  2:1165033/1‑149 (MQ=255)
                                            tACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGa                                                                                                         <  2:199113/149‑1 (MQ=255)
                                                 ttAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCttt                                                                                                    <  1:672097/149‑1 (MQ=255)
                                                     cAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTa                                                                                                <  1:849516/149‑1 (MQ=255)
                                                       ctATGGAACGGCAGATCGGCTAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTTGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAAc                                                                                              <  2:139759/148‑1 (MQ=255)
                                                          tGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGtt                                                                                           <  1:664613/149‑1 (MQ=255)
                                                           gtaaCGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTg                                                                                          <  2:129625/147‑1 (MQ=255)
                                                           ggAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTg                                                                                          <  2:1209158/149‑1 (MQ=255)
                                                           ggAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTg                                                                                          <  2:342204/149‑1 (MQ=255)
                                                           ggAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTg                                                                                          <  2:1495187/149‑1 (MQ=255)
                                                            gAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTgt                                                                                         >  1:770616/1‑149 (MQ=255)
                                                                     aTCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTaa                                                                                <  2:1406517/149‑1 (MQ=255)
                                                                     aTCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTaa                                                                                >  2:1221424/1‑149 (MQ=255)
                                                                     aTCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCAAGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTaa                                                                                >  2:1109823/1‑149 (MQ=255)
                                                                        ggATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATa                                                                             >  2:1505750/1‑149 (MQ=255)
                                                                          aTAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAcc                                                                           <  2:1287956/149‑1 (MQ=255)
                                                                           tAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCg                                                                          <  2:578960/149‑1 (MQ=255)
                                                                            aGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGTTAAATACCgg                                                                         <  2:1358562/149‑1 (MQ=255)
                                                                            aGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCgg                                                                         >  1:113690/1‑149 (MQ=255)
                                                                                     ttttCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGccg                                                                >  1:491358/1‑149 (MQ=255)
                                                                                     ttttCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGccg                                                                >  1:1521897/1‑149 (MQ=255)
                                                                                           ccGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGtt                                                          >  2:588849/1‑149 (MQ=255)
                                                                                                 cTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGcc                                                    >  1:187826/1‑149 (MQ=255)
                                                                                                      aTCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGtct                                               >  1:1512044/1‑149 (MQ=255)
                                                                                                      aTCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGtct                                               >  1:516023/1‑149 (MQ=255)
                                                                                                       tCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGtctc                                              >  1:346821/1‑149 (MQ=255)
                                                                                                           gACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCaaaa                                          <  1:1191567/149‑1 (MQ=255)
                                                                                                            aCCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAAc                                         <  2:541195/149‑1 (MQ=255)
                                                                                                              cATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACca                                       >  2:1441602/1‑149 (MQ=255)
                                                                                                                tCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACcaac                                     >  2:1121695/1‑149 (MQ=255)
                                                                                                                tCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACcaac                                     >  1:1066842/1‑149 (MQ=255)
                                                                                                                      tAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGtt                               >  1:815806/1‑149 (MQ=255)
                                                                                                                       aGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTGCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGttt                              <  1:976589/149‑1 (MQ=255)
                                                                                                                        ggTGTAGTTGATCGACAGTTTCCATTCATCGTTGAAAGGAATTTTCAGTTCGGTTTCCACGCGCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGTTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGtttt                             <  2:467928/149‑1 (MQ=255)
                                                                                                                                   tCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGt                  <  1:429493/149‑1 (MQ=255)
                                                                                                                                     tgaaaattCCGTTCAACGATGAATGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGc                >  2:733893/7‑149 (MQ=255)
                                                                                                                                     gACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGc                >  2:1387369/1‑149 (MQ=255)
                                                                                                                                      gaaaattCCGTTCAACAATGAATGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCa               >  1:227533/6‑149 (MQ=255)
                                                                                                                                      aCAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCa               >  1:1557174/1‑149 (MQ=255)
                                                                                                                                      aCAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCa               >  1:1470429/1‑149 (MQ=255)
                                                                                                                                        aGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCACCAAAGTTTTGGTAGCCCGGTGCAGc             <  2:770616/149‑1 (MQ=255)
                                                                                                                                         aattCCGTTCAACGATGAATGGAATTTTCAGTTCGGTTTCCACGCCCTGTATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTgg                           >  2:123491/3‑134 (MQ=255)
                                                                                                                                         aattCCGTTCAACGATGAATGGAATTTTCAGTTCGGTTTCCACGCCCTGTATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTgg                           <  1:123491/132‑1 (MQ=255)
                                                                                                                                           ttCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGtt          <  2:1512044/149‑1 (MQ=255)
                                                                                                                                           ttCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGtt          <  2:1368960/149‑1 (MQ=255)
                                                                                                                                             ccATTCATCGTTGAACGGAATTTTCTGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGGGCAGCGTTGa        >  2:889509/1‑149 (MQ=255)
                                                                                                                                                ttCATCGTTGTACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGt     <  1:1545989/149‑1 (MQ=255)
                                                                                                                                                   aTCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAg  >  1:707477/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAG  >  NC_000913/2245000‑2245295

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