Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 247,885 A→G 5.7% Q83Q (CAA→CAG rayT → RAYT REP element‑mobilizing transposase; TnpA(REP)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913247,8850AG 5.7% 223.5 / 3.5 70Q83Q (CAA→CAGrayTRAYT REP element‑mobilizing transposase; TnpA(REP)
Reads supporting (aligned to +/- strand):  ref base A (38/28);  new base G (2/2);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

ACGCCATTATTAAAGTTAAGCGAGACAGGCCTTTTGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTG  >  NC_000913/247746‑248027
                                                                                                                                           |                                                                                                                                              
aCGCCATTATTAAAGTTAAGCGAGACAGGCCTTTTGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTAccc                                                                                                                                         <  8:83911/147‑1 (MQ=255)
aCGCCATTATTAAAGTTAAGCGAGACAGGCCTTTTGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTAccc                                                                                                                                         <  3:291380/147‑1 (MQ=255)
   ccATTATTAAAGTTAAGCGAGACAGGCCTTTTGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGc                                                                                                                                     >  2:192078/1‑148 (MQ=255)
       tattaAAGTTAAGCGAGACAGGCCTTTTGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTtgt                                                                                                                                 >  2:276037/1‑148 (MQ=255)
        attaAAGTTAAGCGAGACAGGCCTTTTGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTtgtg                                                                                                                                >  5:10118/1‑148 (MQ=255)
          taAAGTTAAGCGAGACAGGCCTTTTGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTAccc                                                                                                                                         >  3:223175/1‑137 (MQ=255)
          taAAGTTAAGCGAGACAGGCCTTTTGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTAccc                                                                                                                                         <  4:223175/137‑1 (MQ=255)
                tAAGCGAGACAGGCCTTTTGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGaa                                                                                                                         >  4:267269/1‑147 (MQ=255)
                    cGAGACAGGCCTTTTGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAata                                                                                                                    <  1:276037/148‑1 (MQ=255)
                       gaCAGGCCTTTTGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCt                                                                                                                 >  5:258779/1‑148 (MQ=255)
                        aCAGGCCTTTTGACATCAACGCCTGTGTCGTGTTGCGAGAGCATCTCCACTGTATCTGGACATTCCCTGAAGGCTATGATGAGTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTg                                                                                                                <  8:241696/148‑1 (MQ=255)
                            gCCTTTTGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGAt                                                                                                                             >  8:209543/1‑131 (MQ=255)
                            gCCTTTTGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGAt                                                                                                                             <  7:209543/131‑1 (MQ=255)
                                tttGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCAc                                                                                                        >  7:53578/1‑148 (MQ=255)
                                tttGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCAc                                                                                                        >  5:45504/1‑148 (MQ=255)
                                 ttGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAat                                                                                                                     >  4:19193/1‑134 (MQ=255)
                                 ttGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAat                                                                                                                     <  3:19193/134‑1 (MQ=255)
                                  tGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGATTGTGGATTGAAAAATATCTGGCAACCACGt                                                                                                      >  5:92339/1‑148 (MQ=255)
                                   gAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGaaa                                                                                                                        <  2:163264/129‑1 (MQ=255)
                                   gAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGaaa                                                                                                                        >  1:163264/1‑129 (MQ=255)
                                     aaTCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGtttt                                                                                                   <  3:267269/148‑1 (MQ=255)
                                     aaTCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGtttt                                                                                                   <  8:53578/148‑1 (MQ=255)
                                      aTCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGttttt                                                                                                  >  4:137490/1‑148 (MQ=255)
                                      aTCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTAGCTGGACTTTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAACTCAATTTACCCATTCTTGTGGAATGAAAAATATCTGGCAACCACCgtttt                                                                                                  >  2:92273/1‑148 (MQ=255)
                                                           ccAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATa                                                                             >  4:182939/1‑148 (MQ=255)
                                                               agCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCa                                                                          >  7:286802/1‑147 (MQ=255)
                                                                    atGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGaaa                                                                                                                        <  7:190565/96‑1 (MQ=255)
                                                                    atGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGaaa                                                                                                                        >  8:190565/1‑96 (MQ=255)
                                                                    atGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAatat                                                                                                                   <  1:158485/101‑1 (MQ=255)
                                                                    atGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAatat                                                                                                                   >  2:158485/1‑101 (MQ=255)
                                                                    atGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGa                                                                      >  6:151550/1‑146 (MQ=255)
                                                                    atGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGa                                                                      <  5:151550/146‑1 (MQ=255)
                                                                    atGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGAt                                                                     >  1:128707/1‑147 (MQ=255)
                                                                     tGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGAtt                                                                    <  3:137490/147‑1 (MQ=255)
                                                                       cACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATtat                                                                  <  1:192078/147‑1 (MQ=255)
                                                                       cACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATta                                                                   >  1:198776/1‑146 (MQ=255)
                                                                        aCTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAg                                                                >  1:230231/1‑148 (MQ=255)
                                                                           gTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGca                                                             >  4:19347/1‑148 (MQ=255)
                                                                           gTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGca                                                             >  5:178616/1‑148 (MQ=255)
                                                                                    cATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGa                                                    >  7:13880/1‑148 (MQ=255)
                                                                                      ttACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGAt                                                   >  3:246765/1‑147 (MQ=255)
                                                                                       tACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGAtta                                                 >  2:71210/1‑148 (MQ=255)
                                                                                        aCCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCa                                                                                <  4:19604/116‑1 (MQ=255)
                                                                                        aCCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCa                                                                                >  3:19604/1‑116 (MQ=255)
                                                                                        aCCTGAAGGCCATCATGATTTTTCCTCGCGCAGCCCGGAAATTAAAAAGCAATATACCCATGCTTGTGCATTGAAAAATATCTGGCAACCACCTTTTTGCCAGAACCCCAACCGCaa                                                                               >  3:23637/1‑117 (MQ=255)
                                                                                         ccTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGTCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATGATGTTGAttata                                               <  2:198776/148‑1 (MQ=255)
                                                                                           tGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCACCCACGTTTTTGGGAGCCCGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGAttatatt                                             >  5:218405/1‑148 (MQ=255)
                                                                                            gAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGaaa                                                                                                                        <  4:3612/72‑1 (MQ=255)
                                                                                            gAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGaaa                                                                                                                        >  3:3612/1‑72 (MQ=255)
                                                                                                   atgatgatTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATaaa                                     >  8:31606/1‑148 (MQ=255)
                                                                                                      atgatTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATAAATcc                                  <  2:128707/148‑1 (MQ=255)
                                                                                                            tttCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATAGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACAAAAGATTATAGGCAGCATGTTGATTATATTTATATAAATCCAGTa                              >  6:178622/1‑146 (MQ=255)
                                                                                                               ccTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATAAATCCAGTAAAGc                          <  4:246765/147‑1 (MQ=255)
                                                                                                               ccTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATAAATCCAGTAAAGc                          <  8:13880/147‑1 (MQ=255)
                                                                                                                             gAAATTAAAAAGCAGTTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGAttatatttat                                          <  4:287162/117‑1 (MQ=255)
                                                                                                                             gAAATTAAAAAGCAGTTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGAttatatttat                                          <  4:291645/117‑1 (MQ=255)
                                                                                                                             gAAATTAAAAAGCAGTTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGAttatatttat                                          >  3:287162/1‑117 (MQ=255)
                                                                                                                             gAAATTAAAAAGCAGTTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGAttatatttat                                          >  3:291645/1‑117 (MQ=255)
                                                                                                                             gAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAg            >  3:271600/1‑147 (MQ=255)
                                                                                                                               aaTTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCaa         >  8:101007/1‑148 (MQ=255)
                                                                                                                               aaTTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCaa         <  7:101007/148‑1 (MQ=255)
                                                                                                                                 ttAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATAAATcc                                  >  5:35241/1‑121 (MQ=255)
                                                                                                                                 ttAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATAAATcc                                  <  6:35241/121‑1 (MQ=255)
                                                                                                                                 ttAAAAAGCAATTTACCCATGCTTGTGGATTCAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCAagt       >  8:304220/1‑148 (MQ=255)
                                                                                                                                   aaaaaGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATAAATCCAGTAAAGCATg                       >  6:273826/1‑130 (MQ=255)
                                                                                                                                   aaaaaGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATAAATCCAGTAAAGCATg                       <  5:273826/130‑1 (MQ=255)
                                                                                                                                   aaaaaGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCAagtg      <  3:19347/147‑1 (MQ=255)
                                                                                                                                   aaaaaGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCAagtg      <  3:182939/147‑1 (MQ=255)
                                                                                                                                     aaaGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCAagtgagt   <  6:178616/148‑1 (MQ=255)
                                                                                                                                      aaGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCAagtgagtg  <  6:92339/148‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
ACGCCATTATTAAAGTTAAGCGAGACAGGCCTTTTGAAATCAACGCCTGGGTCGTTTTGCCAGAGCATATGCACTGTATCTGGACATTACCTGAAGGCGATGATGATTTTTCCTCGCGCTGGCGGGAAATTAAAAAGCAATTTACCCATGCTTGTGGATTGAAAAATATCTGGCAACCACGTTTTTGGGAGCACGCCATCCGCAATACCAAAGATTATAGGCATCATGTTGATTATATTTATATAAATCCAGTAAAGCATGGTTGGGTAAAGCAAGTGAGTG  >  NC_000913/247746‑248027

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