Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,418 A→T 7.9% intergenic (‑436/+167) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,4180AT 7.9% 221.4 / 6.8 75intergenic (‑436/+167)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (33/36);  new base T (4/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTG  >  NC_000913/987272‑987566
                                                                                                                                                  |                                                                                                                                                    
aTTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAAt                                                                                                                                                     <  6:184125/148‑1 (MQ=255)
 tttCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGACAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATaa                                                                                                                                                   >  1:37220/1‑149 (MQ=255)
    ccTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGACAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGaa                                                                                                                                                >  3:136575/1‑149 (MQ=255)
     cTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAAt                                                                                                                                               >  2:62689‑M2/1‑148 (MQ=255)
         tatTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAAtgcg                                                                                                                                            <  7:168787‑M2/148‑5 (MQ=255)
           tttAAATTAACTATACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAAtgcgga                                                                                                                                          >  2:128478‑M2/1‑142 (MQ=255)
           tttAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAAtgcggat                                                                                                                                         <  7:206484‑M2/149‑8 (MQ=255)
                 ttAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAAtgcggattagttt                                                                                                                                   >  5:217347‑M2/1‑136 (MQ=255)
                    aCTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTg                                                                                                                                >  6:150853/1‑149 (MQ=255)
                     cTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTgc                                                                                                                               <  8:207839/149‑1 (MQ=255)
                      tAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAATAATGCGGATTAGTTTTTgcg                                                                                                                              >  1:170340/1‑149 (MQ=255)
                        aaCGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGcgcg                                                                                                                            >  4:108325/1‑149 (MQ=255)
                            gAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTaaa                                                                                                                        >  8:127016/1‑149 (MQ=255)
                                ttttGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGtt                                                                                                                    >  2:101635/1‑149 (MQ=255)
                                  ttGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTg                                                                                                                   <  2:3752/148‑1 (MQ=255)
                                           aTGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTaa                                                                                                           <  3:158639/147‑1 (MQ=255)
                                              gAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTc                                                                                                      <  2:37220/149‑1 (MQ=255)
                                               aaaCTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATtgcggattagttttttcttagctaatagcacaattttcatac                                                                                                      <  7:204669‑M2/148‑43 (MQ=255)
                                                   tttATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTaa                                                                                                 >  2:87711/1‑149 (MQ=255)
                                                      aTCGCCCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTGTGCGCGTAAATGTTGCATTGTAATAGTCCGTAACAg                                                                                              >  4:130971/1‑149 (MQ=255)
                                                      aTCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACa                                                                                               <  7:118791/148‑1 (MQ=255)
                                                      aTCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACa                                                                                               <  8:88430/148‑1 (MQ=255)
                                                      aTCGACCTCGCACAAAATCTTCTATCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATtgcggattagttttttcttagctaatagcacaattttcatactattttt                                                                                               >  1:12203‑M2/1‑99 (MQ=255)
                                                       tCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATtgcggattagttttttcttagctaatagcacaattttcatactattttttg                                                                                             >  8:204669‑M2/1‑98 (MQ=255)
                                                       tCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAGTAGTTTTTGAGCGTAAATGGTGCAATGTAATATTCCGTAACAg                                                                                              >  6:189600/1‑148 (MQ=255)
                                                         gACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGga                                                                                            <  5:5615/148‑1 (MQ=255)
                                                         gACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGga                                                                                            <  7:127016/148‑1 (MQ=255)
                                                          aCCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAat                                                                                                          <  2:161362/133‑1 (MQ=255)
                                                          aCCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAat                                                                                                          >  1:161362/1‑133 (MQ=255)
                                                           ccTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGgatga                                                                                         <  7:94154/149‑1 (MQ=255)
                                                               gCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATc                                                                                       <  1:70870/147‑1 (MQ=255)
                                                               gCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATc                                                                                       <  5:136665/147‑1 (MQ=255)
                                                               gCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGt                                                                                     <  8:168787/149‑1 (MQ=255)
                                                               gCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGt                                                                                     >  5:222050/1‑149 (MQ=255)
                                                                cacaAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAata                                                                                                         >  1:178641/1‑128 (MQ=255)
                                                                cacaAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAata                                                                                                         <  2:178641/128‑1 (MQ=255)
                                                                cacaAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGtt                                                                                    <  2:213769/149‑1 (MQ=255)
                                                                  caAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATTTAATATTCCGTAACAGGATGATCGTTAt                                                                                  <  5:77057/149‑1 (MQ=255)
                                                                  caAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTAt                                                                                  >  6:46179/1‑149 (MQ=255)
                                                                      aTCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCg                                                                              <  1:101635/149‑1 (MQ=255)
                                                                              tCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCggg                                                                            >  6:178066/1‑143 (MQ=255)
                                                                              tCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCggg                                                                            <  5:178066/143‑1 (MQ=255)
                                                                                aGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCaaa                                                                    <  8:206484/149‑1 (MQ=255)
                                                                                 gACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCaaa                                                                    >  2:242632/1‑148 (MQ=255)
                                                                                  aCATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATa                                                                  <  5:46179/149‑1 (MQ=255)
                                                                                   cATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGTTTGCAAATa                                                                  >  6:66152/1‑148 (MQ=255)
                                                                                         gAATGACAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGTTTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCt                                                           <  4:136575/149‑1 (MQ=255)
                                                                                          aaTGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCtt                                                          <  5:189600/149‑1 (MQ=255)
                                                                                               cAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACt                                                      >  7:76319/1‑148 (MQ=255)
                                                                                                 aaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaa                                                   <  6:222050/149‑1 (MQ=255)
                                                                                                  aaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaa                                                   >  6:654/1‑148 (MQ=255)
                                                                                                      tAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaacaaac                                              <  3:230422/149‑1 (MQ=255)
                                                                                                         tATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaacaaa                                               >  4:26943/1‑145 (MQ=255)
                                                                                                            gAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATaa                                        >  6:166438/1‑149 (MQ=255)
                                                                                                                 ttGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGACTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGa                                   <  2:46741/149‑1 (MQ=255)
                                                                                                                 ttGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGa                                   <  3:220846/149‑1 (MQ=255)
                                                                                                                   gtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACATATAACTTGAc                                  <  8:76319/148‑1 (MQ=255)
                                                                                                                       tATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGaa                              >  5:95769/1‑148 (MQ=255)
                                                                                                                         ttAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAAAAGGATGATCGTTATGTCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGaaaa                            <  5:150853/148‑1 (MQ=255)
                                                                                                                            gCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCATTGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACtt                                     >  5:202143/1‑136 (MQ=255)
                                                                                                                            gCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACtt                                     <  6:202143/136‑1 (MQ=255)
                                                                                                                             cTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATa                        >  3:190314/1‑148 (MQ=255)
                                                                                                                                aGAAAAAACTAATCCGCAATAAGAATCCGTATTAGTTTTTGCGCGTAAATGTTGAAATGTAATATTCCGTAACAGGATGAGCGTTATGGCGGGTTTGAAAATAGATTGCTTGCATTAACAAAAAGATAACTGGACAGAAAAGATaaaa                     >  6:136364/1‑148 (MQ=255)
                                                                                                                                 gcaaaaaCTAATCCGCATTCTTATTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaacaaac                                              >  1:243685/3‑122 (MQ=38)
                                                                                                                                 gcaaaaaCTAATCCGCATTCTTATTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaacaaac                                              <  2:243685/120‑1 (MQ=38)
                                                                                                                                 gAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaaa                    >  6:104572/1‑148 (MQ=255)
                                                                                                                                  aaaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaaa                    >  7:87708/1‑147 (MQ=255)
                                                                                                                                    aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAaga                <  5:166438/149‑1 (MQ=255)
                                                                                                                                    aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAaga                <  3:108325/149‑1 (MQ=255)
                                                                                                                                     aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAagag               <  6:61319/149‑1 (MQ=255)
                                                                                                                                     aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAagag               >  5:237158/1‑149 (MQ=255)
                                                                                                                                         tAATCCGCATTCTTATTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAg           >  8:88112/1‑149 (MQ=255)
                                                                                                                                              cGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGc              <  7:158845/141‑1 (MQ=255)
                                                                                                                                              cGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGc              >  8:158845/1‑141 (MQ=255)
                                                                                                                                              cGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgc      >  3:1553/1‑149 (MQ=255)
                                                                                                                                                  aTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTg  >  5:29051/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTG  >  NC_000913/987272‑987566

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