Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,015,399 Δ1 bp 5.7% coding (240/579 nt) mocA → CTP:molybdopterin cytidylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,015,3990A. 5.7% 232.1 / 6.3 70coding (240/579 nt)mocACTP:molybdopterin cytidylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (36/30);  new base . (2/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CTTGATACAAGTATCAAAAATGCGTTGCAGTTTTGTAGCCGAATTATTTTAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATC  >  NC_000913/3015253‑3015546
                                                                                                                                                  |                                                                                                                                                   
cTTGATACAAGTATCAAAAATGCGTTGCAGTTTTGTAGCCGAATTATTTTAGTCACCGGCTATCGTGGTAATCAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTaa                                                                                                                                                    <  8:184954/148‑1 (MQ=255)
    aTACAAGTATCAAAAATGCGTTGCAGTTTTGTAGCCGAATTATTTTAGTCACTGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATcc                                                                                                                                               >  8:208968/1‑149 (MQ=255)
            aTCAAAAATGCGTTGCAGTTTTGTAGCCGAATTATTTTAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGc                                                                                                                                       >  3:112442/1‑149 (MQ=255)
                       gTTGCAGTTTTGTAGCAGAATTCTTTTAGGCACAGGATATCGTGGTAATGAACTGCAAGAACGTTATGAGAACCAGATCAATATCAGTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCg                                                                                                                            <  7:173507/149‑1 (MQ=255)
                          gCAGTTTTGTAGCCGAATTATTTTAGTCACCGGCTATCGTGGTTATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAg                                                                                                                         >  4:158252/1‑149 (MQ=255)
                           cAGTTTTGTAGCCGAATTATTTTAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGt                                                                                                                        <  3:59400/149‑1 (MQ=255)
                            aGTTTTGTAGCCGATTTATTTTAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTg                                                                                                                       >  6:103009/1‑149 (MQ=255)
                              ttttGTAGCCGAATTATTTTAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGc                                                                                                                      >  8:198794/1‑148 (MQ=255)
                                 tGTAGCCGAATTATTTTAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCaaa                                                                                                                   <  3:232056/148‑1 (MQ=255)
                                     gCCGAATTATTTTAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGTCCGCAGTACCCGCAGTGCAAACAGaa                                                                                                              >  1:166463/1‑149 (MQ=255)
                                      ccGAATTATTTTAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAAc                                                                                                             <  7:198794/149‑1 (MQ=255)
                                       cGAATTATTTTAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACa                                                                                                            >  8:4366/1‑149 (MQ=255)
                                       cGAATTATTTTAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACAt                                                                                                           >  8:99666/1‑149 (MQ=255)
                                       cGAATTATTTTAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACAt                                                                                                           >  8:103061/1‑149 (MQ=255)
                                         aaTTATTTTAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATGACTATTATTCACAGCCCCGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACAtt                                                                                                          >  6:179710/1‑149 (MQ=255)
                                           ttattTTAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGt                                                                                                        >  3:147346/1‑149 (MQ=255)
                                             attTTAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGtt                                                                                                       <  1:33994/148‑1 (MQ=255)
                                              ttttAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGtttt                                                                                                     <  3:82465/149‑1 (MQ=255)
                                              ttttAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACAAGAGCAATATCACTATTATTCACAACCAAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGtttt                                                                                                     >  5:69775/1‑149 (MQ=255)
                                              ttttAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCAAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGtttt                                                                                                     >  5:192891/1‑149 (MQ=255)
                                                   gTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCATACAGAACATTGTTTTCTCAc                                                                                                >  1:147824/1‑149 (MQ=255)
                                                   gTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCAc                                                                                                <  1:154963/149‑1 (MQ=255)
                                                      aCCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCa                                                                                             >  7:465/1‑149 (MQ=255)
                                                          gCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGTGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGt                                                                                         <  2:147507/149‑1 (MQ=255)
                                                             aTCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGa                                                                                       >  4:115646/1‑148 (MQ=255)
                                                              tCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACa                                                                                     >  8:200762/1‑149 (MQ=255)
                                                                gTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACAt                                                                                    <  4:111067/148‑1 (MQ=255)
                                                                gTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATAATGCGCAGGGTTTACTAACCTCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATg                                                                                   <  6:203050/149‑1 (MQ=255)
                                                                   gTAATGAACTGCACGAACGTTATGCAAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCa                                                                                <  6:192891/149‑1 (MQ=255)
                                                                    tAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAAc                                                                              <  7:103061/149‑1 (MQ=255)
                                                                    tAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAAc                                                                              <  7:99666/149‑1 (MQ=255)
                                                                         aaCTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCt                                                                           <  8:85399/148‑1 (MQ=255)
                                                                              cACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCAt                                                                     >  6:39992/1‑149 (MQ=255)
                                                                              cACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCAt                                                                     >  8:140359/1‑149 (MQ=255)
                                                                              cACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCAt                                                                     >  6:214576/1‑149 (MQ=255)
                                                                               aCGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATc                                                                    >  4:224830/1‑149 (MQ=255)
                                                                               aCGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATc                                                                    >  3:106847/1‑149 (MQ=255)
                                                                                      ttATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATAttt                                                             >  6:106008/1‑149 (MQ=255)
                                                                                             aaCCAGAGCAAGATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGaa                                                       >  2:3312/1‑148 (MQ=255)
                                                                                              aCCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGaaa                                                      >  2:211451/1‑148 (MQ=255)
                                                                                              aCCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCACGGTTTACTAACCTCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCCCCATCGATATTTTTAGaaaa                                                     <  6:86955/149‑1 (MQ=255)
                                                                                                  gagCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACATAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAAtgt                                                 >  6:120492/1‑147 (MQ=255)
                                                                                                   agCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTg                                                <  2:147824/149‑1 (MQ=255)
                                                                                                     cAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGt                                              <  4:112442/149‑1 (MQ=255)
                                                                                                        tatCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGa                                                                                       <  2:76315/105‑1 (MQ=255)
                                                                                                        tatCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGa                                                                                       >  1:76315/1‑105 (MQ=255)
                                                                                                           cACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCAACCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTAc                                        >  1:199171/1‑149 (MQ=255)
                                                                                                            aCTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTAc                                        <  3:115646/148‑1 (MQ=255)
                                                                                                              tattatTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGa                                      >  4:139317/1‑148 (MQ=255)
                                                                                                              tattatTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGa                                      >  1:14654/1‑148 (MQ=255)
                                                                                                                  attCCCAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCCCCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAatga                                 >  3:130691/1‑149 (MQ=255)
                                                                                                                     cacaACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAatgatg                               <  2:166463/148‑1 (MQ=255)
                                                                                                                      acaACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGgc                             >  8:101500/1‑149 (MQ=255)
                                                                                                                       caACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGgcg                            <  7:208968/149‑1 (MQ=255)
                                                                                                                          cccAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCaa                         >  4:108326/1‑149 (MQ=255)
                                                                                                                            cAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATa                       <  2:14654/149‑1 (MQ=255)
                                                                                                                              gATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACt                     >  4:104616/1‑149 (MQ=255)
                                                                                                                                ttATGCGCAGGGTTTACTAACATCCGTCCAGGCCGCAGTACCCGCAGGGCAAACAGAACATTGGTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGc                   <  8:175019/149‑1 (MQ=255)
                                                                                                                                ttATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGc                   <  3:108326/149‑1 (MQ=255)
                                                                                                                                   tGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCa                 <  7:200762/148‑1 (MQ=255)
                                                                                                                                     cgcAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACt               <  3:158252/148‑1 (MQ=255)
                                                                                                                                       cAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCa            <  5:103009/149‑1 (MQ=255)
                                                                                                                                        aGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCAt           <  3:104616/149‑1 (MQ=255)
                                                                                                                                         gggTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATa          <  7:4366/149‑1 (MQ=255)
                                                                                                                                          ggTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATaa         >  3:94550/1‑149 (MQ=255)
                                                                                                                                          ggTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATaa         >  3:188079/1‑149 (MQ=255)
                                                                                                                                           gTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAAt        >  1:139160/1‑149 (MQ=255)
                                                                                                                                            tttACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATg       <  2:139160/149‑1 (MQ=255)
                                                                                                                                            tttACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATg       <  8:465/149‑1 (MQ=255)
                                                                                                                                            tttACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATg       <  5:39992/149‑1 (MQ=255)
                                                                                                                                                 tAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATc  >  7:117575/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTTGATACAAGTATCAAAAATGCGTTGCAGTTTTGTAGCCGAATTATTTTAGTCACCGGCTATCGTGGTAATGAACTGCACGAACGTTATGCGAACCAGAGCAATATCACTATTATTCACAACCCAGATTATGCGCAGGGTTTACTAACATCCGTCAAGGCCGCAGTACCCGCAGTGCAAACAGAACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATC  >  NC_000913/3015253‑3015546

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: