Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,247,033 T→C 7.7% intergenic (‑264/+30) cirA ← / ← lysP colicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,247,0330TC 7.7% 187.9 / 4.2 63intergenic (‑264/+30)cirA/lysPcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base C (2/3);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

CCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATAC  >  NC_000913/2246885‑2247174
                                                                                                                                                    |                                                                                                                                             
cccAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGtt                                                                                                                                               <  2:122306/149‑1 (MQ=255)
cccAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGtt                                                                                                                                               <  3:12200/149‑1 (MQ=255)
cccAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGtt                                                                                                                                               >  3:126983/1‑149 (MQ=255)
      aaGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAgg                                                                                                                                         <  1:178576/149‑1 (MQ=255)
      aaGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAgg                                                                                                                                         <  1:411286/149‑1 (MQ=255)
           ggAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCAACCGAGAGGGCttt                                                                                                                                    <  4:157675/149‑2 (MQ=255)
           ggAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCAACCGAGAGGGCttt                                                                                                                                    <  2:46264/149‑2 (MQ=255)
                   ctctGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaaggaa                                                                                                                            >  2:274808/1‑149 (MQ=255)
                    tctGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAg                                                                                                                                   >  3:355180/1‑141 (MQ=255)
                    tctGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAg                                                                                                                                   <  4:355180/141‑1 (MQ=255)
                    tctGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaaggaag                                                                                                                           >  1:532647/1‑149 (MQ=255)
                    tctGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaaggaag                                                                                                                           >  4:426008/1‑149 (MQ=255)
                     ctGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaaggaagg                                                                                                                          <  3:304910/149‑1 (MQ=255)
                       gATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAgg                                                                                                                                         <  3:32906/132‑1 (MQ=255)
                       gATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAgg                                                                                                                                         >  4:32906/1‑132 (MQ=255)
                       gATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGgga                                                                                                                        <  2:82808/149‑1 (MQ=255)
                       gATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGgga                                                                                                                        <  2:81916/149‑1 (MQ=255)
                             tCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAAc                                                                                                                  >  2:48909/1‑149 (MQ=255)
                                   ttttttGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTAttt                                                                                                            >  4:207838/1‑149 (MQ=255)
                                      tttGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCtt                                                                                                         <  1:568202/149‑1 (MQ=255)
                                                 cgTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCg                                                                                              <  2:543617/149‑1 (MQ=255)
                                                  gTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCgg                                                                                             >  2:489940/1‑149 (MQ=255)
                                                     aCGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGaa                                                                                          <  3:207838/149‑1 (MQ=255)
                                                     aCGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGaa                                                                                          <  1:465671/149‑1 (MQ=255)
                                                       gATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAGAAAGCCCTCTCGGTTga                                                                                                                                             >  4:421252/1‑96 (MQ=255)
                                                       gATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTga                                                                                                                                             <  3:421252/96‑1 (MQ=255)
                                                       gATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACt                                                                                        >  2:590142/1‑149 (MQ=255)
                                                             ggaaaCGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCAttt                                                                                  >  4:401744/3‑149 (MQ=255)
                                                             gCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCAttt                                                                                  >  2:503591/1‑149 (MQ=255)
                                                              cAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTc                                                                                 >  1:33108/1‑149 (MQ=255)
                                                                aaCGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGc                                                                               <  2:70254/149‑1 (MQ=255)
                                                                aaCGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGc                                                                               >  3:116309/1‑149 (MQ=255)
                                                                aaCGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGc                                                                               >  3:310380/1‑149 (MQ=255)
                                                                         aaTTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTa                                                                      <  2:28935/149‑1 (MQ=255)
                                                                             aTCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGaa                                                                  <  4:495452/149‑1 (MQ=255)
                                                                             aTCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGaa                                                                  >  4:439534/1‑149 (MQ=255)
                                                                                  tCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAg                                                             >  2:305600/1‑149 (MQ=255)
                                                                                  tCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAg                                                             >  2:23262/1‑149 (MQ=255)
                                                                                    cAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGtt                                                           >  2:380143/1‑149 (MQ=255)
                                                                                     aaaTGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTc                                                          >  3:312488/1‑149 (MQ=255)
                                                                                     aaaTGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTc                                                          >  3:584477/1‑149 (MQ=255)
                                                                                      aaTGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTcc                                                         >  1:346263/1‑149 (MQ=255)
                                                                                       aTGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCAGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCt                                                        <  1:305600/149‑1 (MQ=255)
                                                                                        tGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGACGTGAGTTCCtt                                                       <  3:378331/149‑1 (MQ=255)
                                                                                          ttAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCTAGGAAGGGAGGAAACTTATTTCTTACGGTTCTGCGGGAACTTCATTCCGTTTTACCGTACGAGGTGAGTTCCttat                                                     <  4:379172/149‑3 (MQ=255)
                                                                                          ttAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTa                                                     >  1:198048/1‑149 (MQ=255)
                                                                                           tAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTaa                                                    <  2:33108/149‑1 (MQ=255)
                                                                                                tttATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGc                                               <  2:198048/149‑1 (MQ=255)
                                                                                                     atCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTa                                          <  3:401744/149‑1 (MQ=255)
                                                                                                      tCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAATGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGa                                                              >  3:471800/1‑128 (MQ=255)
                                                                                                      tCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAATGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGa                                                              <  4:471800/128‑1 (MQ=255)
                                                                                                         gCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTaa                                                    <  4:57648/135‑1 (MQ=255)
                                                                                                         gCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTaa                                                    >  3:57648/1‑135 (MQ=255)
                                                                                                                       ggAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAg                        <  2:109521/149‑1 (MQ=255)
                                                                                                                             ggAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAg                  >  3:490027/1‑149 (MQ=255)
                                                                                                                                aTTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCgg               >  2:2415/1‑149 (MQ=255)
                                                                                                                                     aaaaGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATAc          >  3:326912/1‑149 (MQ=255)
                                                                                                                                      taaGCCCTCTCAACCGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTa                                                     >  4:540226/2‑105 (MQ=255)
                                                                                                                                      taaGCCCTCTCAACCGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATAcc         >  2:326876/2‑149 (MQ=255)
                                                                                                                                      taaGCCCTCTCAACCGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAATTCCTTTa                                                     <  3:540226/104‑1 (MQ=255)
                                                                                                                                      aaaGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATAcc         <  2:91430/149‑1 (MQ=255)
                                                                                                                                      aaaGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATAcc         >  3:355519/1‑149 (MQ=255)
                                                                                                                                            ctctCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAAtata   <  1:274808/149‑1 (MQ=255)
                                                                                                                                             tctcGGTTGATAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATAc  <  2:539837/149‑1 (MQ=255)
                                                                                                                                             tctcGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATAc  >  4:16645/1‑149 (MQ=255)
                                                                                                                                             tctcGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATAc  >  2:334360/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
CCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATAC  >  NC_000913/2246885‑2247174

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