Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,915,373 C→T 8.3% intergenic (‑173/+180) glmU ← / ← atpC fused N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase/F1 sector of membrane‑bound ATP synthase, epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,915,3730CT 8.3% 213.8 / 9.3 72intergenic (‑173/+180)glmU/atpCfused N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase/F1 sector of membrane‑bound ATP synthase, epsilon subunit
Reads supporting (aligned to +/- strand):  ref base C (26/40);  new base T (3/3);  total (29/43)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.34e-01

CGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCA  >  NC_000913/3915226‑3915512
                                                                                                                                                   |                                                                                                                                           
cGAATTTAAACCGCTTCACCTTGAAAAAACAACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCa                                                                                                                                            <  1:405447/149‑1 (MQ=255)
    ttaaaaCCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTTCAGCAAATTAAGGCAGTAc                                                                                                                                        >  2:404768/4‑149 (MQ=255)
        aaCCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTa                                                                                                                                    >  4:281041/1‑149 (MQ=255)
           cGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAAtt                                                                                                                                 <  3:15311/149‑1 (MQ=255)
               tCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGTACTGCCTTAAttttct                                                                                                                             <  3:522616/149‑4 (MQ=255)
                cACCTTGAAAAAAGTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGc                                                                                                                            <  4:554585/149‑1 (MQ=255)
                cACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGc                                                                                                                            <  3:553212/149‑1 (MQ=255)
                 aCCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGTGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGAGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCg                                                                                                                           >  2:106205/1‑149 (MQ=255)
                   gatGAAAAAAATGTAGTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTc                                                                                                                         >  1:548846/3‑149 (MQ=255)
                   cTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTc                                                                                                                         >  3:353656/1‑149 (MQ=255)
                    ttGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCa                                                                                                                        >  1:152490/1‑149 (MQ=255)
                    ttGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCa                                                                                                                        >  2:15303/1‑149 (MQ=255)
                     tGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAg                                                                                                                       <  3:281041/149‑1 (MQ=255)
                            aCTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATg                                                                                                                >  2:472541/1‑149 (MQ=255)
                                  tttttttCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTg                                                                                                          >  1:339927/1‑149 (MQ=255)
                                   ttttttCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGc                                                                                                         <  2:502780/149‑1 (MQ=255)
                                       ttCATCGTGAAATGGACAGAGGATTTATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCtt                                                                                                     <  1:532120/149‑1 (MQ=255)
                                           tCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTAttt                                                                                                 >  3:280367/1‑149 (MQ=255)
                                           tCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTAttt                                                                                                 <  1:610071/149‑1 (MQ=255)
                                            cGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTc                                                                                                <  1:199731/149‑1 (MQ=255)
                                            cGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTc                                                                                                <  1:472541/149‑1 (MQ=255)
                                             gTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTc                                                                                                <  4:263665/148‑1 (MQ=255)
                                             gTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTc                                                                                                >  3:263665/1‑148 (MQ=255)
                                               gAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGAt                                                                                             <  1:156296/149‑1 (MQ=255)
                                               gAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGAt                                                                                             <  2:548846/149‑1 (MQ=255)
                                                   tGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGATGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAAt                                                                                         <  4:486825/149‑1 (MQ=255)
                                                   tGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAAt                                                                                         <  4:23273/149‑1 (MQ=255)
                                                          agGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATaaa                                                                                  <  1:15303/149‑1 (MQ=255)
                                                          agGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATaaa                                                                                  <  2:315456/149‑1 (MQ=255)
                                                                  ttGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGaaaa                                                                          <  4:353656/149‑1 (MQ=255)
                                                                     ttCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAGTAAACAGGAAAACAg                                                                       >  2:273639/1‑149 (MQ=255)
                                                                     ttCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAGTAAACAGGAAAACAg                                                                       >  2:588061/1‑149 (MQ=255)
                                                                      tCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAGTAAACAGGAAAACAgg                                                                      <  1:588061/149‑1 (MQ=255)
                                                                      tCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAGTAAACAGGAAAACAgg                                                                      <  1:273639/149‑1 (MQ=255)
                                                                       cAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAggg                                                                     <  4:210267/149‑1 (MQ=255)
                                                                        aaTAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAgggg                                                                    >  1:106678/1‑149 (MQ=255)
                                                                         aTAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAggggg                                                                   <  3:390245/149‑1 (MQ=255)
                                                                             cGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGttt                                                               <  1:163463/149‑1 (MQ=255)
                                                                                 ttATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTtata                                                           <  3:11909/149‑1 (MQ=255)
                                                                                  caTCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTatat                                                          <  3:470087/148‑1 (MQ=255)
                                                                                              cGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTa                                              <  3:4316/149‑1 (MQ=255)
                                                                                               gAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTAt                                             <  2:624156/149‑1 (MQ=255)
                                                                                                     gcCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGt                                                                                                                     >  4:492208/1‑71 (MQ=255)
                                                                                                     gcCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGt                                                                                                                     <  3:492208/71‑1 (MQ=255)
                                                                                                     gcCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGtt                                                                >  1:19313/1‑124 (MQ=255)
                                                                                                     gcCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGtt                                                                <  2:19313/124‑1 (MQ=255)
                                                                                                              gACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGTTGTTTTTGCTCTTATTTCGATCAATGAATAAACGGGAAAACAGGGGGTTTTTATATCAGCAGGATCTATGTGAACGCTATTCAg                              <  1:404768/149‑1 (MQ=255)
                                                                                                               aCCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAgg                             >  2:606946/1‑149 (MQ=255)
                                                                                                               aCCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAgg                             >  3:177593/1‑149 (MQ=255)
                                                                                                                     gTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGt                       <  1:475919/149‑1 (MQ=255)
                                                                                                                       cACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTca                     >  4:387492/1‑149 (MQ=255)
                                                                                                                        aCGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTcac                    <  4:232910/149‑1 (MQ=255)
                                                                                                                          gCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTAt                                             >  4:610800/1‑122 (MQ=255)
                                                                                                                          gCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTAt                                             <  3:610800/122‑1 (MQ=255)
                                                                                                                           ccTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTAt                                  <  2:324857/132‑1 (MQ=255)
                                                                                                                           ccTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTAt                                  >  1:324857/1‑132 (MQ=255)
                                                                                                                            cTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATca                                                        >  3:151269/1‑109 (MQ=255)
                                                                                                                            cTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATca                                                        <  4:151269/109‑1 (MQ=255)
                                                                                                                                  tCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaa          >  2:412193/1‑149 (MQ=255)
                                                                                                                                     gCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATcagcag                                                    <  1:470347/104‑1 (MQ=255)
                                                                                                                                     gCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATcagcag                                                    >  2:470347/1‑104 (MQ=255)
                                                                                                                                      ctgAAAATCAAGGCTGTACCTTAATTTAAGCCTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACCGGAAAACAGGGCGGTTTCATATCAGCAGGATCTACGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaaa      <  3:383564/147‑1 (MQ=255)
                                                                                                                                      cAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaaa      <  1:106205/149‑1 (MQ=255)
                                                                                                                                       cttaaaTTAAGGTACTGCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATca                                                        <  2:83556/95‑1 (MQ=39)
                                                                                                                                       cttaaaTTAAGGTACTGCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATca                                                        >  1:83556/4‑98 (MQ=38)
                                                                                                                                       cttaaaTTAAGGTACTGCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAag     >  4:224607/4‑149 (MQ=255)
                                                                                                                                       aGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAag     <  3:387492/149‑1 (MQ=255)
                                                                                                                                        ttaaaTTAAGGTACTGCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagc    >  1:139441/3‑149 (MQ=255)
                                                                                                                                        ttaaaTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaag      <  1:570086/145‑2 (MQ=255)
                                                                                                                                        ttaaaTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaag      >  2:570086/3‑146 (MQ=255)
                                                                                                                                         aaaaTTAAGGCAGTAGCTTAATTTAAGCGTAAGGTCGATGTTTTTGCTCTTATTTCGATCAATGAATAAGCAGGCAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGCACGGGTCACACGCGTAAAAAAAagag   >  4:383564/1‑147 (MQ=255)
                                                                                                                                          aaaTTAAGGTACTGCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagcca  <  3:224607/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
CGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCA  >  NC_000913/3915226‑3915512

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