Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,245,143 A→G 7.7% W543R (TGG→CGG)  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,1430AG 7.7% 188.7 / 5.8 65W543R (TGG→CGG) cirAcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter
Reads supporting (aligned to +/- strand):  ref base A (32/28);  new base G (3/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

GAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGA  >  NC_000913/2244995‑2245289
                                                                                                                                                    |                                                                                                                                                  
gAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCa                                                                                                                                                    <  2:721651/149‑1 (MQ=255)
 aaaCATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCAt                                                                                                                                                   <  2:1064554/149‑1 (MQ=255)
   aCATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGGATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCTTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTc                                                                                                                                                 <  1:1031264/149‑1 (MQ=255)
     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:342204/149‑1 (MQ=255)
                                                                ggAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTg                                                                                    <  2:1495187/149‑1 (MQ=255)
                                                                ggAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTg                                                                                    <  2:1209158/149‑1 (MQ=255)
                                                                 gAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTgt                                                                                   >  1:770616/1‑149 (MQ=255)
                                                                          aTCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTaa                                                                          >  2:1221424/1‑149 (MQ=255)
                                                                          aTCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTaa                                                                          <  2:1406517/149‑1 (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:1521897/1‑149 (MQ=255)
                                                                                          ttttCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGccg                                                          >  1:491358/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:1470429/1‑149 (MQ=255)
                                                                                                                                           aCAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCa         >  1:1557174/1‑149 (MQ=255)
                                                                                                                                             aGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCACCAAAGTTTTGGTAGCCCGGTGCAGc       <  2:770616/149‑1 (MQ=255)
                                                                                                                                              aattCCGTTCAACGATGAATGGAATTTTCAGTTCGGTTTCCACGCCCTGTATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTgg                     <  1:123491/132‑1 (MQ=255)
                                                                                                                                              aattCCGTTCAACGATGAATGGAATTTTCAGTTCGGTTTCCACGCCCTGTATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTgg                     >  2:123491/3‑134 (MQ=255)
                                                                                                                                                ttCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGtt    <  2:1512044/149‑1 (MQ=255)
                                                                                                                                                ttCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGtt    <  2:1368960/149‑1 (MQ=255)
                                                                                                                                                  ccATTCATCGTTGAACGGAATTTTCTGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGGGCAGCGTTGa  >  2:889509/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGA  >  NC_000913/2244995‑2245289

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: