Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,771,583 C→T 5.4% L146L (CTG→TTG)  yfjT → CP4‑57 prophage; putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,771,5830CT 5.4% 239.2 / 2.6 74L146L (CTG→TTG) yfjTCP4‑57 prophage; putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (35/35);  new base T (2/2);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-01

CATTCCCTCAACCCTTGGCTGGCGTAACGGAGATAATGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAGGCCACGTACTTTCGAACGATGAGCGCAGTGG  >  NC_000913/2771435‑2771729
                                                                                                                                                    |                                                                                                                                                  
cATTCCCTCAACCCTTGGCTGGCGTAACGGAGATAATGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTcc                                                                                                                                                    >  3:195434/1‑149 (MQ=255)
cATTCCCTCAACCCTTGGCTGGCGTAACGGAGATAATGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTcc                                                                                                                                                    >  2:377552/1‑149 (MQ=255)
    cccTCAACCCTTGGCTGGCGTAACGGAGATAATGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGc                                                                                                                                                 <  8:83139/148‑1 (MQ=255)
      cTCAACCCTTGGCTGGCGTAACGGAGATAATGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCaa                                                                                                                                              >  1:53243/1‑149 (MQ=255)
         aaCCCTTGGCTGGCGTAACGGAGATAATGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACAt                                                                                                                                           <  4:100047/149‑1 (MQ=255)
         aaCCCTTGGCTGGCGTAACGGAGATAATGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGAACAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAATAGGTCCTGCCAACAt                                                                                                                                           <  2:269805/149‑1 (MQ=255)
           cccTTGGCTGGCGTAACGGAGATAATGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCTAACATTa                                                                                                                                         <  2:22623/149‑1 (MQ=255)
           cccTTGGCTGGCGTAACGGAGATAATGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTa                                                                                                                                         <  4:195434/149‑1 (MQ=255)
            ccTTGGCTGGCGTAACGGAGATAATGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAg                                                                                                                                        >  8:294474/1‑149 (MQ=255)
                     gcGTAACGGAGATAATGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTg                                                                                                                                 <  7:294474/147‑1 (MQ=255)
                          aCGGAGATAATGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAAc                                                                                                                          <  1:135348/149‑1 (MQ=255)
                          aCGGAGATAATGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAAc                                                                                                                          <  8:184496/149‑1 (MQ=255)
                            ggAGATAATGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGa                                                                                                                        >  1:170559/1‑149 (MQ=255)
                            ggAGATAATGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGa                                                                                                                        >  4:149913/1‑149 (MQ=255)
                            ggAGATAATGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGa                                                                                                                        <  3:149913/149‑1 (MQ=255)
                                  aaTGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACg                                                                                                                  >  6:4689/1‑149 (MQ=255)
                                    tGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCTTGCCAACATTAGCGGACTGCACTAACGACTGACGGc                                                                                                                >  5:64899/1‑149 (MQ=255)
                                       cTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGTTCCTGCCAACATTAGCGGACTGCAATAACGACTGACGGCTgc                                                                                                             >  4:291144/1‑149 (MQ=255)
                                       cTGGTTCAGCTGCATCATGCATATCCGTAAGGCCCGGCAGTTCGCCCTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGACTGCACAACAGGGCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTgc                                                                                                             <  8:356454/149‑1 (MQ=255)
                                          gTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTgcgc                                                                                                           <  7:313517/148‑1 (MQ=255)
                                               gCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCt                                                                                                               >  8:392522/1‑139 (MQ=255)
                                               gCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCt                                                                                                               <  7:392522/139‑1 (MQ=255)
                                               gCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTccc                                                                                                     <  6:6086/149‑1 (MQ=255)
                                                   catcatGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCTTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTa                                                                                                   <  6:64899/147‑1 (MQ=255)
                                                      catGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAAc                                                                                                                                             >  7:374730/1‑102 (MQ=255)
                                                      catGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAAc                                                                                                                                             <  8:374730/102‑1 (MQ=255)
                                                        tGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATaa                                                                                             >  2:276543/1‑148 (MQ=255)
                                                        tGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCCTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATaa                                                                                             >  7:11504/1‑148 (MQ=255)
                                                             atCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCaa                                                                                                                                              >  3:304626/1‑94 (MQ=255)
                                                             atCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCaa                                                                                                                                              <  4:304626/94‑1 (MQ=255)
                                                              tCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCt                                                                                      <  5:4689/149‑1 (MQ=255)
                                                               ccGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTc                                                                                                       >  1:338137/1‑131 (MQ=255)
                                                               ccGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTc                                                                                                       <  2:338137/131‑1 (MQ=255)
                                                                 gTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGAtt                                                                                  <  7:180620/149‑1 (MQ=255)
                                                                  tAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGAtt                                                                                  <  1:377552/149‑1 (MQ=255)
                                                                      gCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCt                                                                              >  1:300070/1‑149 (MQ=255)
                                                                       cccGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTc                                                                             <  8:5235/149‑1 (MQ=255)
                                                                         cGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCcg                                                                           >  4:365315/1‑149 (MQ=255)
                                                                           gCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCcggc                                                                         >  4:371639/1‑149 (MQ=255)
                                                                               ttCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTc                                                                     >  5:391902/1‑149 (MQ=255)
                                                                                  gACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGCCGGCGGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCgg                                                                   <  1:276543/148‑1 (MQ=255)
                                                                                   aCGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCgggg                                                                 <  2:53243/149‑1 (MQ=255)
                                                                                      tCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGtt                                                              >  5:9759/1‑149 (MQ=255)
                                                                                        tACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTc                                                            >  5:224454/1‑149 (MQ=255)
                                                                                        tACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTc                                                            >  5:401774/1‑149 (MQ=255)
                                                                                         aCGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCt                                                           >  6:209610/1‑149 (MQ=255)
                                                                                           gTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCttt                                                         <  2:170559/149‑1 (MQ=255)
                                                                                            tCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCtttt                                                        <  5:209610/149‑1 (MQ=255)
                                                                                                tGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGCCGGTCGGGGGTTTCATTTTtct                                                    >  7:208879/1‑149 (MQ=255)
                                                                                                  aCCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTtctc                                                   <  6:391902/148‑1 (MQ=255)
                                                                                                    cAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCTTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCttt                                                         <  1:172038/140‑1 (MQ=255)
                                                                                                    cAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCTTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCttt                                                         >  2:172038/1‑140 (MQ=255)
                                                                                                    cAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCt                                                           <  3:182009/138‑1 (MQ=255)
                                                                                                    cAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCt                                                           >  4:182009/1‑138 (MQ=255)
                                                                                                      gAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTaa                                              >  7:230088/1‑149 (MQ=255)
                                                                                                              gATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCag                                       <  6:224454/148‑1 (MQ=255)
                                                                                                                   tttAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAg                                 <  3:365315/149‑1 (MQ=255)
                                                                                                                   tttAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAg                                 >  5:282197/1‑149 (MQ=255)
                                                                                                                   tttAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCACTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAg                                 >  2:302810/1‑149 (MQ=255)
                                                                                                                        cccTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAGGCCAc                            >  6:183864/1‑149 (MQ=255)
                                                                                                                         ccTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTATCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAGGCCAc                            >  5:187973/1‑148 (MQ=255)
                                                                                                                         ccTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAAtctc                                          <  5:6054/134‑1 (MQ=255)
                                                                                                                         ccTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAAtctc                                          >  6:6054/1‑134 (MQ=255)
                                                                                                                           tGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAGCAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAGGCCACGTa                         >  6:289025/1‑149 (MQ=255)
                                                                                                                             aTGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAGGCCACGTACt                       <  5:183864/149‑1 (MQ=255)
                                                                                                                               gTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAGGCCACGTACttt                     >  5:86591/1‑149 (MQ=255)
                                                                                                                                   ggAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAGGCCACGTACTTTCGaa                 <  8:11504/149‑1 (MQ=255)
                                                                                                                                         cacaAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAGGCCACGTACTTTCGAACGAt             >  5:78054/1‑147 (MQ=255)
                                                                                                                                         cacaAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAGGCCACGTACTTTCGAACGATGa           <  3:371639/149‑1 (MQ=255)
                                                                                                                                          acaAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCcgg                                                                          >  1:314547/1‑85 (MQ=255)
                                                                                                                                          acaAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCcgg                                                                          <  2:314547/85‑1 (MQ=255)
                                                                                                                                           caAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATATCAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAGGCCACGTACttt                     <  5:196133/137‑1 (MQ=255)
                                                                                                                                           caAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATATCAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAGGCCACGTACttt                     >  6:196133/1‑137 (MQ=255)
                                                                                                                                            aaaaGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAGGCCACGTACTTTCGAACGATGAgcg        >  1:375817/1‑149 (MQ=255)
                                                                                                                                                 gTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAGGCCACGTACTTTCGAACGATGAGCGCAGTg   <  6:187973/149‑1 (MQ=255)
                                                                                                                                                  tCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACACCCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAGGCCACGTACTTTCGAACGATGAGCGCAGTgg  >  4:32873/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CATTCCCTCAACCCTTGGCTGGCGTAACGGAGATAATGCCTGGTTCAGCTTCATCATGGATATCCGTAAGGCCCGGCAGTTCGACGTCTACGTCAATGACCAGAAAGTTGGATCGTTTAACCCTGATGTTCGGAATGCACAAAAGGTCCTGCCAACATTAGCGGACTGCACTAACGACTGACGGCTGCGCCTTCCCTACCATAACAACCCTGATTGCCTCCGGCGGTCGGGGGTTTCTTTTTTCTCCTTAATCTCCAGAGGTAGGCCACGTACTTTCGAACGATGAGCGCAGTGG  >  NC_000913/2771435‑2771729

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: