Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,412,644 A→C 5.8% intergenic (‑50/+33) yfbU ← / ← yfbV UPF0304 family protein/UPF0208 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,412,6440AC 5.8% 308.3 / 8.0 86intergenic (‑50/+33)yfbU/yfbVUPF0304 family protein/UPF0208 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (40/41);  new base C (2/3);  total (42/44)
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.94e-01

AACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGC  >  NC_000913/2412510‑2412755
                                                                                                                                      |                                                                                                               
aaCGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCcttctt                                                                                                        <  3:380889/144‑1 (MQ=255)
aaCGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCcttctt                                                                                                        >  4:380889/1‑144 (MQ=255)
    ttCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGc                                                                                                >  2:462660/1‑148 (MQ=255)
     tCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTa                                                                                                                 <  3:851208/130‑1 (MQ=255)
     tCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTa                                                                                                                 >  4:851208/1‑130 (MQ=255)
         cATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCcttctt                                                                                                        <  4:607498/135‑1 (MQ=255)
         cATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCcttctt                                                                                                        >  3:607498/1‑135 (MQ=255)
         cATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCATACTCAGCCTTCttt                                                                                                       <  2:586553/136‑1 (MQ=255)
         cATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCATACTCAGCCTTCttt                                                                                                       >  1:586553/1‑136 (MQ=255)
                         cATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATg                                                                                                               >  3:867086/1‑112 (MQ=255)
                         cATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATg                                                                                                               <  4:867086/112‑1 (MQ=255)
                               catcatCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACaaa                                                                          <  1:123974/143‑1 (MQ=255)
                               catcatCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACaaa                                                                          >  2:123974/1‑143 (MQ=255)
                               catcatCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAatcatc                                                                     >  1:780380/1‑148 (MQ=255)
                                 tcatcTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCctt                                                                                                           >  3:902086/1‑108 (MQ=255)
                                 tcatcTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCctt                                                                                                           <  4:902086/108‑1 (MQ=255)
                                  catcTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCATACTCAGCCTTCttttt                                                                                                     >  2:82687/1‑112 (MQ=38)
                                  catcTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCATACTCAGCCTTCttttt                                                                                                     <  1:82687/113‑2 (MQ=39)
                                       tGTACTTATTAAACAAAATAAGACGTTGAGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCtttt                                                                                                      <  3:1077282/107‑1 (MQ=255)
                                       tGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCtttt                                                                                                      >  4:1077282/1‑107 (MQ=255)
                                       tGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAg                                                                                        >  3:529815/1‑121 (MQ=255)
                                       tGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAg                                                                                        <  4:529815/121‑1 (MQ=255)
                                        gTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGt                                                            <  3:340199/148‑1 (MQ=255)
                                               ttAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTtat                                                                                                    >  4:326880/1‑101 (MQ=255)
                                               ttAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTtat                                                                                                    <  3:326880/101‑1 (MQ=255)
                                                   aCAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCcttct                                                                                                         >  2:804413/1‑92 (MQ=255)
                                                   aCAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCcttct                                                                                                         <  1:804413/92‑1 (MQ=255)
                                                   aCAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCATACTCAGCCTTCttt                                                                                                       <  1:428557/94‑1 (MQ=38)
                                                    cAAAATAAGATGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACaa                                                                           <  1:850524/121‑1 (MQ=255)
                                                    cAAAATAAGATGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACaa                                                                           >  2:850524/1‑121 (MQ=255)
                                                    cAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTAt                                                                                                                >  1:834420/1‑84 (MQ=255)
                                                    cAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTAt                                                                                                                <  2:834420/84‑1 (MQ=255)
                                                      aaaTAAGACGTTGCGCGTTGGTCATTTCCATTGTTTACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCtttt                                                                                                      <  3:465975/92‑1 (MQ=255)
                                                      aaaTAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCtttt                                                                                                      >  4:465975/1‑92 (MQ=255)
                                                      aaaTAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCa                                                                                         >  2:367421/1‑105 (MQ=255)
                                                      aaaTAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCa                                                                                         <  1:367421/105‑1 (MQ=255)
                                                           aGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAgg                                                                                       <  1:825932/102‑1 (MQ=255)
                                                           aGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAgg                                                                                       >  2:825932/1‑102 (MQ=255)
                                                            gACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCa                                                                                         >  4:138067/1‑99 (MQ=255)
                                                            gACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCa                                                                                         <  3:138067/99‑1 (MQ=255)
                                                                        tGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCg                            <  1:462660/148‑1 (MQ=255)
                                                                        tGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCg                            <  2:780380/148‑1 (MQ=255)
                                                                           tCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCaa                                                                   <  1:567535/106‑1 (MQ=255)
                                                                           tCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCaa                                                                   >  2:567535/1‑106 (MQ=255)
                                                                              tttCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTtat                                                                                                    >  1:1008046/1‑70 (MQ=255)
                                                                              tttCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTtat                                                                                                    <  2:1008046/70‑1 (MQ=255)
                                                                              tttCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAATCGTGTCAg                      >  2:917237/1‑148 (MQ=255)
                                                                              tttCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAg                      >  3:724922/1‑148 (MQ=255)
                                                                                    ttgttgACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGttt                                                          >  2:769809/1‑106 (MQ=255)
                                                                                    ttgttgACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGttt                                                          <  1:769809/106‑1 (MQ=255)
                                                                                     tgttgACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGt                                                                                    <  4:444010/79‑1 (MQ=255)
                                                                                     tgttgACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGt                                                                                    >  3:444010/1‑79 (MQ=255)
                                                                                        tgACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCttt                                                                                                       <  3:836085/57‑1 (MQ=255)
                                                                                        tgACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCttt                                                                                                       >  4:836085/1‑57 (MQ=255)
                                                                                        tgACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGaa                                          >  3:273691/1‑118 (MQ=255)
                                                                                        tgACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGaa                                          <  4:273691/118‑1 (MQ=255)
                                                                                             ccTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGtttt                                                         >  2:624586/1‑98 (MQ=255)
                                                                                             ccTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGtttt                                                         <  1:624586/98‑1 (MQ=255)
                                                                                                 tATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAgg   <  4:724922/148‑1 (MQ=255)
                                                                                                  aTCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTACTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGc  <  1:917237/148‑1 (MQ=255)
                                                                                                   tCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCa                                                    <  2:192783/97‑1 (MQ=255)
                                                                                                   tCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAg                              <  1:452996/119‑1 (MQ=255)
                                                                                                   tCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAg                              >  2:452996/1‑119 (MQ=255)
                                                                                                   tCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTACCCa                                                    >  1:192783/1‑97 (MQ=255)
                                                                                                       tctACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGt                                                            <  1:946101/85‑1 (MQ=255)
                                                                                                       tctACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGt                                                            >  2:946101/1‑85 (MQ=255)
                                                                                                       tctACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttgt                                               <  2:224637/98‑1 (MQ=255)
                                                                                                       tctACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttgt                                               >  1:224637/1‑98 (MQ=255)
                                                                                                        ctactaCGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAg                                                             <  1:1063800/83‑1 (MQ=255)
                                                                                                        ctactaCGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAg                                                             >  2:1063800/1‑83 (MQ=255)
                                                                                                               ggTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTa                                                        <  3:819823/81‑1 (MQ=255)
                                                                                                               agTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTCTTa                                                        >  4:819823/2‑81 (MQ=255)
                                                                                                                 tGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCaa                                                                   <  3:461545/68‑1 (MQ=255)
                                                                                                                 tGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCaa                                                                   >  4:461545/1‑68 (MQ=255)
                                                                                                                    aaaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCgtgt                         >  2:631162/1‑107 (MQ=255)
                                                                                                                    aaaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCgtgt                         <  1:631162/107‑1 (MQ=255)
                                                                                                                          agaagGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttgtt                                              >  2:423682/1‑80 (MQ=255)
                                                                                                                          agaagGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttgtt                                              <  1:423682/80‑1 (MQ=255)
                                                                                                                           gaagGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGcc              >  3:366255/1‑111 (MQ=255)
                                                                                                                           gaagGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGcc              <  4:366255/111‑1 (MQ=255)
                                                                                                                            aagGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgt                                                  >  2:353014/1‑74 (MQ=255)
                                                                                                                            aagGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgt                                                  <  1:353014/74‑1 (MQ=255)
                                                                                                                                cggAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGa                                           >  2:425756/3‑77 (MQ=255)
                                                                                                                                cTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGa                                           <  1:425756/77‑1 (MQ=255)
                                                                                                                                 tGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGc  <  3:988467/117‑1 (MQ=255)
                                                                                                                                 tGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGc  >  4:988467/1‑117 (MQ=255)
                                                                                                                                  gAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAAtt                                                                               >  2:1120522/1‑39 (MQ=255)
                                                                                                                                  gAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAAtt                                                                               <  1:1120522/39‑1 (MQ=255)
                                                                                                                                      |                                                                                                               
AACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGC  >  NC_000913/2412510‑2412755

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: