Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,915,377 A→G 8.3% intergenic (‑177/+176) 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,3770AG 8.3% 217.3 / 8.9 72intergenic (‑177/+176)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 A (27/39);  new base G (3/3);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

TTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCC  >  NC_000913/3915230‑3915522
                                                                                                                                                   |                                                                                                                                                 
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:588061/1‑149 (MQ=255)
                                                                 ttCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAGTAAACAGGAAAACAg                                                                                 >  2:273639/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                                                                          <  2:19313/124‑1 (MQ=255)
                                                                                                 gcCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGtt                                                                          >  1:19313/1‑124 (MQ=255)
                                                                                                          gACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGTTGTTTTTGCTCTTATTTCGATCAATGAATAAACGGGAAAACAGGGGGTTTTTATATCAGCAGGATCTATGTGAACGCTATTCAg                                        <  1:404768/149‑1 (MQ=255)
                                                                                                           aCCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAgg                                       >  3:177593/1‑149 (MQ=255)
                                                                                                           aCCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAgg                                       >  2:606946/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                                                       <  3:610800/122‑1 (MQ=255)
                                                                                                                      gCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTAt                                                       >  4:610800/1‑122 (MQ=255)
                                                                                                                       ccTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTAt                                            <  2:324857/132‑1 (MQ=255)
                                                                                                                       ccTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTAt                                            >  1:324857/1‑132 (MQ=255)
                                                                                                                        cTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATca                                                                  <  4:151269/109‑1 (MQ=255)
                                                                                                                        cTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATca                                                                  >  3:151269/1‑109 (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                                                                  >  1:83556/4‑98 (MQ=38)
                                                                                                                                   cttaaaTTAAGGTACTGCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATca                                                                  <  2:83556/95‑1 (MQ=39)
                                                                                                                                   cttaaaTTAAGGTACTGCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAag               >  4:224607/4‑149 (MQ=255)
                                                                                                                                   aGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAag               <  3:387492/149‑1 (MQ=255)
                                                                                                                                    ttaaaTTAAGGTACTGCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagc              >  1:139441/3‑149 (MQ=255)
                                                                                                                                    ttaaaTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaag                >  2:570086/3‑146 (MQ=255)
                                                                                                                                    ttaaaTTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCaaaaaaag                <  1:570086/145‑2 (MQ=255)
                                                                                                                                     aaaaTTAAGGCAGTAGCTTAATTTAAGCGTAAGGTCGATGTTTTTGCTCTTATTTCGATCAATGAATAAGCAGGCAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGCACGGGTCACACGCGTAAAAAAAagag             >  4:383564/1‑147 (MQ=255)
                                                                                                                                      aaaTTAAGGTACTGCCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAagcca            <  3:224607/149‑1 (MQ=255)
                                                                                                                                                aGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTcc  >  2:564665/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAATTAAGGCAGTACCTTAATTTAAGCGTCAGGTGGATGTTTTTGCTCTTATTTCGATCAATGAATAAACAGGAAAACAGGGGGGTTTTATATCAGCAGGATCTATGTGAACGCTATTCAGGACGGGTCACACGCGCAAAAAAAAGCCAGCCTGTTTCC  >  NC_000913/3915230‑3915522

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: