Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,245,156 C→T 7.1% P538P (CCG→CCA)  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,1560CT 7.1% 215.1 / 4.8 71P538P (CCG→CCA) cirAcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter
Reads supporting (aligned to +/- strand):  ref base C (39/26);  new base T (3/2);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

ACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGC  >  NC_000913/2245008‑2245304
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                                                               >  1:495305/1‑149 (MQ=255)
       ttCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAAttt                                                                                                                                               >  2:429493/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:1495187/149‑1 (MQ=255)
                                                   ggAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTg                                                                                                   <  2:342204/149‑1 (MQ=255)
                                                   ggAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTg                                                                                                   <  2:1209158/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:516023/1‑149 (MQ=255)
                                                                                              aTCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGtct                                                        >  1:1512044/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                                              >  1:1066842/1‑149 (MQ=255)
                                                                                                        tCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACcaac                                              >  2:1121695/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:1368960/149‑1 (MQ=255)
                                                                                                                                   ttCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGtt                   <  2:1512044/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)
                                                                                                                                              gTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACg        >  2:621317/1‑149 (MQ=255)
                                                                                                                                               ttGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGt       >  1:654924/1‑149 (MQ=255)
                                                                                                                                               ttGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGt       >  2:753961/1‑149 (MQ=255)
                                                                                                                                               ttGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGt       >  2:879014/1‑149 (MQ=255)
                                                                                                                                                 gAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTAc     >  1:1035178/1‑149 (MQ=255)
                                                                                                                                                  aaCGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGa          <  1:1478871/143‑1 (MQ=255)
                                                                                                                                                  aaCGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGa          >  2:1478871/1‑143 (MQ=255)
                                                                                                                                                  aaCGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACg    >  2:762640/1‑149 (MQ=255)
                                                                                                                                                    cGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGgc  >  1:1318900/1‑149 (MQ=255)
                                                                                                                                                    cGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGATAAATTCCGGTATGCGCCGTCCGTTAGAGCCCGTCTCAAAACCAACAAAGATTATGTAGACCGGTGCAGCGTTGACTTCAGACGTAAgga  >  1:275099/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGC  >  NC_000913/2245008‑2245304

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: ‑