Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,603,279 A→G S69G (AGT→GGT)  b1516 → putative LACI‑type transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,603,2790AG86.1% 543.1 / 73.3 210S69G (AGT→GGT) b1516putative LACI‑type transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (11/18);  new base G (67/113);  total (79/131)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

Only 100 of 214 total aligned reads displayed.

GCGTGGGATTTTTTACCAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCAGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTGATGGCTTG  >  NC_000913/1603187‑1603368
                                                                                            |                                                                                         
   tGGGATTTTTTACCAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGgtg                                                                                        >  1:3126026/1‑93 (MQ=255)
   tGGGATTTTTTACCAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGgtg                                                                                        >  2:560249/1‑93 (MQ=255)
   tGGGATTTTTTACCAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGgtg                                                                                        >  2:1609571/1‑93 (MQ=255)
    gggATTTTTTACCAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGgtgt                                                                                       <  1:530746/93‑1 (MQ=255)
     ggATTTTTTACCAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGtt                                                                                      <  1:358223/93‑1 (MQ=255)
     ggATTTTTTACCAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGtt                                                                                      <  2:5367680/93‑1 (MQ=255)
       aTTTTTTACCAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTc                                                                                    <  1:1410439/93‑1 (MQ=255)
       aTTTTTTACCAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTc                                                                                    >  1:1693915/1‑93 (MQ=255)
         tttttACCAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTg                                                                                  <  2:852674/93‑1 (MQ=255)
          ttttACCAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTgg                                                                                 <  1:5100466/93‑1 (MQ=255)
             tACCAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCa                                                                              <  1:5648624/93‑1 (MQ=255)
              aCCAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAg                                                                             <  1:3217119/93‑1 (MQ=255)
              aCCAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCAGTGTTTCTGGTCAg                                                                             <  1:1005937/93‑1 (MQ=255)
               ccAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAgg                                                                            <  1:6241975/93‑1 (MQ=255)
               ccAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAgg                                                                            >  2:4822257/1‑93 (MQ=255)
               ccAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAgg                                                                            <  1:497/93‑1 (MQ=255)
                cAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGt                                                                           >  2:5648663/1‑93 (MQ=255)
                cAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGt                                                                           <  1:115415/93‑1 (MQ=255)
                 aGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCAGTGTTTCTGGTCAGGTa                                                                          >  2:5943063/1‑93 (MQ=255)
                    ggTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAg                                                                       <  2:2238405/93‑1 (MQ=255)
                    ggTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAg                                                                       >  1:4584032/1‑93 (MQ=255)
                     gTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGt                                                                      <  1:1043891/93‑1 (MQ=255)
                      tGGCAACGGCGCACACCAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCAGTGTTTCTGGTCAGGTACAGtt                                                                     <  1:5793564/93‑1 (MQ=255)
                      tGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGtt                                                                     <  1:4902748/93‑1 (MQ=255)
                        gCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGa                                                                   <  1:6136135/93‑1 (MQ=255)
                          aaCGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGAtt                                                                 <  1:6281707/93‑1 (MQ=255)
                           aCGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGGTGATta                                                                >  1:3306148/1‑93 (MQ=255)
                            cGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATtaa                                                               <  2:1909408/93‑1 (MQ=255)
                            cGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATtaa                                                               <  2:2770406/93‑1 (MQ=255)
                            cGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCAGTGTTTCTGGTCAGGTACAGTTGATtaa                                                               <  1:3262664/93‑1 (MQ=255)
                              gcgcACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATtaata                                                             >  1:2274304/1‑93 (MQ=255)
                              gcgcACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCAGTGTTTCTGGTCAGGTACAGTTGATtaata                                                             >  1:2976030/1‑93 (MQ=255)
                               cgcACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTGCGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATtaataa                                                            <  2:2518662/93‑1 (MQ=255)
                                gcACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAAc                                                           >  1:2955658/1‑93 (MQ=255)
                                gcACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAAc                                                           <  1:319285/93‑1 (MQ=255)
                                  acaacaAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACtt                                                         <  1:1341086/93‑1 (MQ=255)
                                  acaacaAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACtt                                                         <  2:3559444/93‑1 (MQ=255)
                                   caacaaGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTc                                                        >  1:5842768/1‑93 (MQ=255)
                                   caacaaGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTc                                                        >  2:5804098/1‑93 (MQ=255)
                                   caacaaGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCAGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTc                                                        <  1:4447912/93‑1 (MQ=255)
                                   caacaaGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCAGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTc                                                        <  2:3785185/93‑1 (MQ=255)
                                      caaGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGtc                                                     >  2:3711182/1‑93 (MQ=255)
                                       aaGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGtca                                                    <  1:3993250/93‑1 (MQ=255)
                                        aGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGtcaa                                                   >  2:3044273/1‑93 (MQ=255)
                                         gCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGtcaat                                                  >  1:4199636/1‑93 (MQ=255)
                                          caggTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGtcaatc                                                 <  2:3249268/91‑1 (MQ=255)
                                           gggTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGtcaatca                                                <  1:1739925/93‑1 (MQ=255)
                                           gggTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGtcaatca                                                >  1:4217468/1‑93 (MQ=255)
                                           gggTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGtcaatca                                                <  1:6450011/93‑1 (MQ=255)
                                           gggTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGtcaatca                                                >  1:4801592/1‑93 (MQ=255)
                                            ggTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGtcaatcaa                                               <  2:416526/93‑1 (MQ=255)
                                            ggTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCAGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGtcaatcaa                                               <  2:426361/93‑1 (MQ=255)
                                             gTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAg                                              <  2:6005571/93‑1 (MQ=255)
                                              tAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAgg                                             <  1:4114203/93‑1 (MQ=255)
                                                aaGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGtt                                           >  1:1241425/1‑93 (MQ=255)
                                                 agagCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTta                                          <  1:5169891/93‑1 (MQ=255)
                                                 agagCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTta                                          <  2:1321318/93‑1 (MQ=255)
                                                    gCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATaa                                       <  2:3165283/93‑1 (MQ=255)
                                                       gggCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGc                                    >  2:2806705/1‑93 (MQ=255)
                                                         gCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCa                                  >  1:2757444/1‑93 (MQ=255)
                                                         gCGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCa                                  >  1:2527014/1‑93 (MQ=255)
                                                         gCGTTGATGTGACCTACGACGGGCCGACAGAACCCAGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCa                                  >  2:4552324/1‑93 (MQ=255)
                                                          cGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCAt                                 <  2:3942455/93‑1 (MQ=255)
                                                          cGTTGATGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCAt                                 <  1:6616650/93‑1 (MQ=255)
                                                               aTGTGACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCg                            <  1:4392680/93‑1 (MQ=255)
                                                                tgtgACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGt                           <  2:5033603/93‑1 (MQ=255)
                                                                tgtgACCTACGACGGGCCGACAGAACCCAGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGt                           >  2:163184/1‑93 (MQ=255)
                                                                  ggaacTACGACGGGCCGACAGAACCCGGTCTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGttt                         <  2:1759843/89‑1 (MQ=255)
                                                                   gACCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTc                        <  1:4422569/93‑1 (MQ=255)
                                                                    aCCTACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCt                       <  1:1604507/93‑1 (MQ=255)
                                                                       tACGACGGGCCGTCATATCTCTGTGTACTTGGTCTGGTACAGTTGATTAATTACTTCGTTTATCAAGGTTTTAACGCTACTATCGATTCTGCg                    >  2:6649530/1‑93 (MQ=255)
                                                                       tACGACGGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCg                    >  2:2974477/1‑93 (MQ=255)
                                                                        acgacgGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCgg                   >  2:5619597/1‑93 (MQ=255)
                                                                        acgacgGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCgg                   >  2:1528761/1‑93 (MQ=255)
                                                                        acgacgGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCgg                   <  1:928819/93‑1 (MQ=255)
                                                                         cgacgGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGt                  <  1:3702176/93‑1 (MQ=255)
                                                                          gacgGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGtt                 <  1:2622799/93‑1 (MQ=255)
                                                                           acgGGCCGACAGACCCCAGTCTTTCTGGTCACGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGttt                <  1:1555337/93‑1 (MQ=255)
                                                                            cgGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTc               <  1:6648023/93‑1 (MQ=255)
                                                                            cgGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTc               >  2:2189063/1‑93 (MQ=255)
                                                                            cgGGCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTc               <  2:4039135/93‑1 (MQ=255)
                                                                             gggCCGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCg              >  1:2405825/1‑93 (MQ=255)
                                                                                ccGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCt           <  2:755410/93‑1 (MQ=255)
                                                                                 cGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTg          >  1:3022824/1‑93 (MQ=255)
                                                                                 cGACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTg          <  2:3345836/93‑1 (MQ=255)
                                                                                  gACAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTGa         >  2:4932134/1‑93 (MQ=255)
                                                                                   aCAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTGAt        <  1:3808028/93‑1 (MQ=255)
                                                                                   aCAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTGAt        <  2:156113/93‑1 (MQ=255)
                                                                                    cAGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTGATg       >  2:1433951/1‑93 (MQ=255)
                                                                                     aGAACCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTGATgg      >  1:2074646/1‑93 (MQ=255)
                                                                                      gAACCCAGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTGATGGc     <  2:4493358/93‑1 (MQ=255)
                                                                                       aaCCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTGATGGCt    <  1:3628721/93‑1 (MQ=255)
                                                                                       aaCCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTGATGGCt    <  2:4653896/93‑1 (MQ=255)
                                                                                       aaCCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGATATAACGCCATTATCGTTTCTGCGGTTTCGCCTGATGGCt    >  2:6142778/1‑93 (MQ=255)
                                                                                       aaCCCAGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTGATGGCt    <  1:757230/93‑1 (MQ=255)
                                                                                        aCCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTGATGGCtt   <  2:4388590/93‑1 (MQ=255)
                                                                                        aCCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTGATGGCtt   >  1:4705722/1‑93 (MQ=255)
                                                                                        aCCCGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTGATGGCtt   <  2:3500340/93‑1 (MQ=255)
                                                                                         cccGGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTGATGGCTtg  <  1:5035887/93‑1 (MQ=255)
                                                                                         cccAGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTGATGGCTtg  <  2:373532/93‑1 (MQ=255)
                                                                                            |                                                                                         
GCGTGGGATTTTTTACCAGCGGTGGCAACGGCGCACAACAAGCGGGTAAAGAGCTGGGCGTTGATGTGACCTACGACGGGCCGACAGAACCCAGTGTTTCTGGTCAGGTACAGTTGATTAATAACTTCGTCAATCAAGGTTATAACGCCATTATCGTTTCTGCGGTTTCGCCTGATGGCTTG  >  NC_000913/1603187‑1603368

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