Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,070,230 T→C 5.9% intergenic (‑50/‑47) insC1 ← / → yoeA IS2 repressor TnpA/CP4‑44 prophage; putative disrupted hemin or colicin receptor;Phage or Prophage Related; interrupted by IS2 and C‑terminal deletion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,070,2300TC 5.9% 203.9 / 3.0 68intergenic (‑50/‑47)insC1/yoeAIS2 repressor TnpA/CP4‑44 prophage; putative disrupted hemin or colicin receptor;Phage or Prophage Related; interrupted by IS2 and C‑terminal deletion
Reads supporting (aligned to +/- strand):  ref base T (33/31);  new base C (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

TCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAG  >  NC_000913/2070126‑2070372
                                                                                                        |                                                                                                                                              
tCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTg                                                                                                    >  8:184883/1‑149 (MQ=11)
  ttttCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGc                                                                                                  >  1:141859/1‑149 (MQ=32)
  ttttCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGc                                                                                                  >  1:47469/1‑149 (MQ=32)
    ttCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCg                                                                                                 >  4:287026/1‑148 (MQ=32)
      ccTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGaa                                                                                               <  6:221059/148‑1 (MQ=34)
             ggTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCAtcctc                                                                                       <  6:75569/149‑1 (MQ=40)
               tACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTccc                                                                                     >  3:239913/1‑149 (MQ=255)
                aCGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCg                                                                                    >  4:188035/1‑149 (MQ=255)
                   ccGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCgtgt                                                                                 >  6:103998/1‑149 (MQ=255)
                    cGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGtt                                                                                >  6:271876/1‑149 (MQ=255)
                      tCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGttt                                                                               >  6:386910/1‑148 (MQ=255)
                        tGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTaa                                                                            >  3:191503/1‑149 (MQ=255)
                            tttCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGAGTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGGGTTTTAATtct                                                                        >  1:367061/1‑149 (MQ=255)
                              tctcCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATtctc                                                                       <  2:7903/148‑1 (MQ=255)
                                tcCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATtc                                                                         <  3:156790/144‑1 (MQ=255)
                                tcCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATtc                                                                         >  4:156790/1‑144 (MQ=255)
                                tcCGGCCCTAAGACATCAATCATCTGTTCTCCAATGAATAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCTTGTTTTAATTCTCAtt                                                                    >  5:59732/1‑149 (MQ=255)
                                 ccGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCAtt                                                                    >  4:35313/1‑148 (MQ=255)
                                  cGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCAtt                                                                    <  6:58278/147‑1 (MQ=255)
                                   ggCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGAt                                                                 <  7:236436/149‑1 (MQ=255)
                                      ccTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGt                                                              >  7:32347/1‑149 (MQ=255)
                                       cTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTc                                                             <  1:6104/149‑1 (MQ=255)
                                         aaGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCag                                                           >  6:54956/1‑149 (MQ=255)
                                           gACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGGCagga                                                         <  2:45331/149‑1 (MQ=255)
                                            aCATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCaggag                                                        <  7:184883/149‑1 (MQ=255)
                                                caatcaTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGGaa                                                    >  1:311749/1‑149 (MQ=255)
                                                 aatcaTATTTTCTCCATTGACTAGTCTATAATCTAGTATTAAGACTATCACTTATTTAATTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAAc                                                   <  6:10754/149‑1 (MQ=21)
                                                   tcatcTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACtt                                                 <  3:188035/149‑1 (MQ=255)
                                                   tcatcTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACt                                                  >  8:353113/1‑148 (MQ=255)
                                                      tcTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTAt                                               >  8:398593/1‑148 (MQ=255)
                                                        tGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAg                                            >  4:148219/1‑149 (MQ=255)
                                                         gTTCTCCAATGACTAGTCTAAAAACTAGTATTTAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGTGAAGCAGCCTCCCGTGTTTTTATTATCATTGATGGTCATGAGGTAACTTATCAgc                                           >  1:122095/1‑149 (MQ=18)
                                                          ttCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATc                                              <  8:139810/145‑1 (MQ=255)
                                                          ttCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAgcg                                          >  3:63791/1‑149 (MQ=255)
                                                          ttCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAgcg                                          <  3:118728/149‑1 (MQ=255)
                                                          ttCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAgc                                           <  7:353113/148‑1 (MQ=255)
                                                          ttCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAg                                            <  7:283986/147‑1 (MQ=255)
                                                               cAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATg                                                                <  6:389465/122‑1 (MQ=255)
                                                               cAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATg                                                                >  5:389465/1‑122 (MQ=255)
                                                               cAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCa                                             <  5:191914/141‑1 (MQ=255)
                                                               cAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCa                                             >  6:191914/1‑141 (MQ=255)
                                                                aaTGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCgg                                    <  2:141859/149‑1 (MQ=255)
                                                                 aTGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGga                                   <  4:191503/149‑1 (MQ=255)
                                                                      tAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATaa                              >  8:378057/1‑149 (MQ=255)
                                                                        gTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAAtt                            <  7:185198/149‑1 (MQ=255)
                                                                         tCTAAAAACTAGTATTAAGACTATCACTTAGTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTa                           >  7:226043/1‑149 (MQ=255)
                                                                           tAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATg                         >  4:347921/1‑149 (MQ=255)
                                                                            aaaaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATgg                        >  6:248160/1‑149 (MQ=255)
                                                                             aaaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTAt                          <  7:398593/146‑1 (MQ=255)
                                                                             aaaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATg                         <  2:345922/147‑1 (MQ=255)
                                                                              aaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGgtg                      >  3:334227/1‑149 (MQ=255)
                                                                               aaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTTGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGTTCAGGAGGTAACTTATCAGCGCGCCGGAGTTAATTATGgtgt                     <  8:245297/149‑1 (MQ=255)
                                                                                 cTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTgg                   <  5:103998/149‑1 (MQ=255)
                                                                                     tATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACt               >  6:119193/1‑149 (MQ=255)
                                                                                       ttAAGACTATCACTTATCTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATg                         >  6:356457/1‑137 (MQ=255)
                                                                                       ttAAGACTATCACTTATCTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATg                         <  3:130831/137‑1 (MQ=255)
                                                                                       ttAAGACTATCACTTATCTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATg                         >  4:130831/1‑137 (MQ=255)
                                                                                       ttAAGACTATCACTTATCTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATg                         <  5:356457/137‑1 (MQ=255)
                                                                                        tAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGtt            >  2:348643/1‑149 (MQ=255)
                                                                                         aaGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGtt            <  4:334227/148‑1 (MQ=255)
                                                                                          aGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGa          <  1:29136/149‑1 (MQ=255)
                                                                                          aGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGa          <  3:347921/149‑1 (MQ=255)
                                                                                           gACTTTCCCTTTTTTATGTGATATTGGTTGTCTGGAGATTCCGGGGGCCAGGCTAGTGGCTAAGCATCCTGCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCTCGCCGGAGATAATTATGGTGTGGGACTGTTGAt         <  1:257238/149‑1 (MQ=38)
                                                                                            aCTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATa        <  5:54956/149‑1 (MQ=255)
                                                                                              tATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGa          >  2:264414/1‑145 (MQ=255)
                                                                                                tCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGAt     <  5:386910/148‑1 (MQ=255)
                                                                                                  aCTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGa   <  5:271876/148‑1 (MQ=255)
                                                                                                  aCTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAAGTATGGTGTGGGACTGTTGATAGATGAg  >  4:236225/1‑149 (MQ=255)
                                                                                                        |                                                                                                                                              
TCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAG  >  NC_000913/2070126‑2070372

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