Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,478,527 Δ1 bp intergenic (+36/+23) yjgN → / ← valS hypothetical protein/valine tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,478,5270G.99.5% 947.4 / ‑3.8 204intergenic (+36/+23)yjgN/valShypothetical protein/valine tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (107/96);  total (107/97)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 215 total aligned reads displayed.

AGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCA  >  NC_000913/4478437‑4478614
                                                                                          |                                                                                       
aGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGttt                                                                                    >  1:712338/1‑93 (MQ=255)
   cttcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtca                                                                                 <  2:204261/93‑1 (MQ=255)
    ttcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcat                                                                                <  2:5130580/93‑1 (MQ=255)
    ttcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcat                                                                                >  2:4021424/1‑93 (MQ=255)
     tcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTTGTTTtcatc                                                                               >  2:1353961/1‑93 (MQ=255)
       ttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCAc                                                                             <  2:1136789/93‑1 (MQ=255)
       ttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCAc                                                                             >  2:3572939/1‑93 (MQ=255)
       ttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGGGGATGCTCCCGCCCTCCTTTTCATTAc                                                                             >  1:2393127/1‑93 (MQ=255)
        tAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACt                                                                            <  1:5738158/93‑1 (MQ=255)
         aGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtg                                                                           <  2:4646003/93‑1 (MQ=255)
          gCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgt                                                                          >  1:863512/1‑93 (MQ=255)
           ccAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgtg                                                                         >  1:2082749/1‑93 (MQ=255)
           ccAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgtg                                                                         <  2:2569667/93‑1 (MQ=255)
           ccAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgtg                                                                         <  1:2907489/93‑1 (MQ=255)
             aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtt                                                                       >  2:2862067/1‑93 (MQ=255)
             aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtt                                                                       >  2:4793924/1‑93 (MQ=255)
               tATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtttt                                                                     >  1:331934/1‑93 (MQ=255)
                   tCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGAtt                                                                 >  2:4678918/1‑93 (MQ=255)
                    cACGGGGAGTTATGCCTTCTTTACAATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTa                                                                <  1:187552/93‑1 (MQ=255)
                     agtgggAGTGATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTAc                                                               <  2:307579/90‑1 (MQ=255)
                       ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAg                                                             >  2:933614/1‑93 (MQ=255)
                        gggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                                                            >  1:4889377/1‑93 (MQ=255)
                        gggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                                                            <  2:1659680/93‑1 (MQ=255)
                        gggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGATTTCATCACTGTGTTTTGATTACAgc                                                            >  2:1855598/1‑93 (MQ=255)
                        gagAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                                                            <  2:6149623/91‑1 (MQ=255)
                             ttATGCTGTCTTTACCATTTCTGAAAGAGATAAAAAAGGCAGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCgcggc                                                       <  1:156407/93‑1 (MQ=255)
                             ttATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCgcggc                                                       <  2:5007131/93‑1 (MQ=255)
                             ttATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCgcggc                                                       <  1:3393285/93‑1 (MQ=255)
                              tATTCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCTTGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCtcggcg                                                      >  1:5488390/1‑93 (MQ=255)
                               aTGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGa                                                     <  2:5198805/93‑1 (MQ=255)
                               aTGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGa                                                     <  2:5486905/93‑1 (MQ=255)
                               aTGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGa                                                     >  2:3001829/1‑93 (MQ=255)
                                  ccTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                                                  >  2:2282601/1‑93 (MQ=255)
                                  ccTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                                                  <  1:5067653/93‑1 (MQ=255)
                                   cttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAAc                                                 <  1:1427803/93‑1 (MQ=255)
                                   cttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAAc                                                 <  1:511043/93‑1 (MQ=255)
                                   cttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAAc                                                 >  1:5926883/1‑93 (MQ=255)
                                    ttcttTCCAATTTCTGTAAGAGATAAAAAAGGCGGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACa                                                <  1:2529005/93‑1 (MQ=255)
                                    ttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACa                                                <  2:1475691/93‑1 (MQ=255)
                                     tcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATCACAg                                               >  1:4036712/1‑93 (MQ=255)
                                     tcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAg                                               >  2:2726792/1‑93 (MQ=255)
                                     tcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAg                                               >  1:1023036/1‑93 (MQ=255)
                                      cttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGc                                              >  2:3156092/1‑93 (MQ=255)
                                       tttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGcc                                             <  1:1362040/93‑1 (MQ=255)
                                       tttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGcc                                             >  1:5623897/1‑93 (MQ=255)
                                        ttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCct                                            >  2:2956660/1‑93 (MQ=255)
                                        ttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCct                                            <  1:5273595/93‑1 (MQ=255)
                                        ttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCct                                            >  1:1623699/1‑93 (MQ=255)
                                         tACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAACctg                                           >  1:3621693/1‑93 (MQ=255)
                                          tccATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgc                                          <  1:2156455/92‑1 (MQ=255)
                                          aCCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgc                                          <  2:5234564/93‑1 (MQ=255)
                                          aCCATTTCTGTAAGAGAAAAAAAAGGCCGGAGCATGCCCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgc                                          <  2:4503012/93‑1 (MQ=255)
                                           ccATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTTCAGCGCGGCGATAACAGCctgct                                         >  2:124586/1‑93 (MQ=255)
                                           ccATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgct                                         >  2:3907160/1‑93 (MQ=255)
                                           ccATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgct                                         >  2:3661678/1‑93 (MQ=255)
                                            cATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGTctgctg                                        >  1:694049/1‑93 (MQ=255)
                                            cATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgctg                                        >  2:1600589/1‑93 (MQ=255)
                                              tttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGtt                                      >  1:4837348/1‑93 (MQ=255)
                                               ttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTc                                     <  1:1134837/93‑1 (MQ=255)
                                                  tGTAAGAGATAAAAAGGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAAt                                  <  2:6032746/93‑1 (MQ=255)
                                                  tGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAAt                                  <  2:5669790/93‑1 (MQ=255)
                                                   gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATc                                 >  2:1838073/1‑93 (MQ=255)
                                                    tAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCa                                <  1:3141408/93‑1 (MQ=255)
                                                    tAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCa                                >  1:2633968/1‑93 (MQ=255)
                                                    tAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCa                                >  1:4640686/1‑93 (MQ=255)
                                                      agagATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGt                              <  2:1067669/93‑1 (MQ=255)
                                                      agagATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGt                              >  1:2205469/1‑93 (MQ=255)
                                                      agagATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGt                              >  1:509703/1‑93 (MQ=255)
                                                       gagaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt                             <  1:414123/93‑1 (MQ=255)
                                                         agataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGttt                            <  1:3860766/89‑1 (MQ=255)
                                                         agataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGttt                            <  1:6061065/89‑1 (MQ=255)
                                                          gataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttc                           >  2:181000/4‑93 (MQ=255)
                                                          gataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttc                           <  2:4524378/90‑1 (MQ=255)
                                                          gataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttc                           <  2:2409756/90‑1 (MQ=255)
                                                           ataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcg                          >  1:980952/3‑93 (MQ=255)
                                                            taaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgc                         <  1:5334663/92‑1 (MQ=255)
                                                            taaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgc                         >  1:2897129/2‑93 (MQ=255)
                                                             aaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgct                        >  1:1735372/1‑93 (MQ=255)
                                                                aaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctttc                     >  1:3539578/1‑93 (MQ=255)
                                                                 aaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctttcg                    >  2:5281736/1‑93 (MQ=255)
                                                                    gCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTtcgctttcgctt                 >  2:3027151/1‑93 (MQ=255)
                                                                       cggagcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTcc               >  2:349773/5‑93 (MQ=255)
                                                                        ggagcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCg              <  2:1671645/90‑1 (MQ=255)
                                                                        cgagcatgcTCCGGCGTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCg              <  1:4727307/90‑1 (MQ=255)
                                                                            catgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATa          <  2:3386786/93‑1 (MQ=255)
                                                                             atgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAg         >  1:4506434/1‑93 (MQ=255)
                                                                              tgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGc        >  2:4195359/1‑93 (MQ=255)
                                                                              tgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGc        >  2:1202890/1‑93 (MQ=255)
                                                                               gcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGcc       <  2:3770565/93‑1 (MQ=255)
                                                                                cTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCt      <  1:2698800/93‑1 (MQ=255)
                                                                                cTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCt      >  1:2654653/1‑93 (MQ=255)
                                                                                cTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCt      <  2:1435389/93‑1 (MQ=255)
                                                                                cTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCt      >  2:686419/1‑93 (MQ=255)
                                                                                cTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCt      >  1:145827/1‑93 (MQ=255)
                                                                                 tCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCtt     >  2:1025931/1‑93 (MQ=255)
                                                                                  ccGGCCTTCGTTTTCATCGCTGTGTTTTGATTGCAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTc    <  2:3330214/93‑1 (MQ=255)
                                                                                  ccGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTc    <  1:4249746/93‑1 (MQ=255)
                                                                                   cGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTcc   <  1:3968978/93‑1 (MQ=255)
                                                                                   cGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTcc   <  1:631411/93‑1 (MQ=255)
                                                                                    ggCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCa  >  2:5380096/1‑93 (MQ=255)
                                                                                          |                                                                                       
AGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCA  >  NC_000913/4478437‑4478614

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: