Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,430,948 T→G 5.1% A156A (GCA→GCC tsaC ← tRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein, threonine‑dependent ADP‑forming ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,430,9480TG 5.1% 320.3 / 6.0 109A156A (GCA→GCCtsaCtRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein, threonine‑dependent ADP‑forming ATPase
Reads supporting (aligned to +/- strand):  ref base T (48/54);  new base G (2/4);  total (51/58)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

Only 100 of 112 total aligned reads displayed.

GTGGGCTATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAA  >  NC_000913/3430800‑3431097
                                                                                                                                                    |                                                                                                                                                     
gtgGGCTATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATtgtg                                                                                                                                                      >  2:1217157/1‑150 (MQ=255)
  gggCTATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTgcg                                                                                                                                                    >  1:1018283/1‑150 (MQ=255)
      tATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCACAATTTGGcgcggc                                                                                                                                                <  1:476250/150‑3 (MQ=255)
      tATCGGATTACCAAAAACAGCATAGGTTTCCAATATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGaa                                                                                                                                                <  2:726861/150‑1 (MQ=255)
       aTCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCACAATTTGGcgcggcg                                                                                                                                               <  1:877219/150‑4 (MQ=255)
       aTCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCACAATTTGGcgcggcg                                                                                                                                               <  1:677131/150‑4 (MQ=255)
       aTCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATGTCTGAAGGATTTAAACGCCCCCCCGGTTCACCAGGCACACCCGGGAACGCCGCGCCACATTGTGCGCGacc                                                                                                                                               >  1:348707/1‑148 (MQ=255)
              tACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCt                                                                                                                                        <  2:900700/150‑1 (MQ=255)
               aCCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTa                                                                                                                                       >  1:650785/1‑150 (MQ=255)
                ccAAACACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTGCACCCGTCAGGGCTTCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACATCCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTAc                                                                                                                                      <  2:1656375/150‑1 (MQ=255)
                ccAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCGCAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTAc                                                                                                                                      >  2:1526585/1‑150 (MQ=255)
                ccAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTAc                                                                                                                                      >  2:1198101/1‑150 (MQ=255)
                ccAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTAc                                                                                                                                      <  1:1720125/150‑1 (MQ=255)
                ccAAAAACAGCATAGGTATCCATTATGTTACCCCTGTCGAAACACTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTAc                                                                                                                                      <  1:494205/150‑1 (MQ=255)
                  aaaacCAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTTAAGGATTTAAACGCCCCCCCGTTTCACCTGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTg                                                                                                                                    <  1:1034736/150‑1 (MQ=255)
                  aaaacCAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCTGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTg                                                                                                                                    <  2:1056504/150‑1 (MQ=255)
                  aaaaaCAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTg                                                                                                                                    <  2:1056735/150‑1 (MQ=255)
                  aaaaaCAGCATAGGTTTCAATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGACCTTCGTCTAcag                                                                                                                                    >  1:250868/1‑148 (MQ=255)
                    aaaCAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACACCCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGtt                                                                                                                                  >  2:1757182/1‑150 (MQ=255)
                     aaCAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTc                                                                                                                                 <  2:1758479/150‑1 (MQ=255)
                      aCAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCg                                                                                                                                <  2:1637893/150‑1 (MQ=255)
                         gCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACa                                                                                                                             <  2:2048572/150‑1 (MQ=255)
                         gCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGGTCACCCGGCAGGGCATCGCGGATTTCTGACGGATTTCAACGCCCCCCCGTGTCACCAGCCACACCCGGGAACGCCGCGCCAAATTGAGCGCGAGAGTCGGCGACTGTTCGAcc                                                                                                                             >  1:370389/1‑149 (MQ=255)
                             aGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGt                                                                                                                         <  1:840068/150‑1 (MQ=255)
                              ggTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTg                                                                                                                        >  1:1916342/1‑150 (MQ=255)
                                tttCCATTATGTTACCCCTGGCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCGGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGc                                                                                                                      <  2:1751752/150‑1 (MQ=255)
                                gttCCATTATGTTACCCCTGTCGTAACAGTTCACCCGTCAGGGCATCGCGGATGTCGGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGc                                                                                                                      <  1:1340984/149‑1 (MQ=255)
                                   ccATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACCACCGGGAACGCCGCGCCAAATTGTGCGCGACCTTCGTCTCCTGTTCGACAAGGAGGCAAt                                                                                                                   >  2:1455200/1‑150 (MQ=255)
                                         tGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTc                                                                                                             <  1:866676/150‑1 (MQ=255)
                                              ccccTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGtt                                                                                                        >  1:1497605/1‑150 (MQ=255)
                                              ccccTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGtt                                                                                                        >  2:528887/1‑150 (MQ=255)
                                               cccTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGCACTAAGACTACTGTTCGACAACGTGGCAATCCACTCAAGTTg                                                                                                       >  2:1407977/1‑150 (MQ=255)
                                               cccTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTg                                                                                                       <  1:1503825/150‑1 (MQ=255)
                                               cccTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTg                                                                                                       >  2:1075672/1‑150 (MQ=255)
                                                      taaaCAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTgg                                                                                                <  2:1117974/149‑1 (MQ=255)
                                                         aCAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAg                                                                                             <  2:1732834/150‑1 (MQ=255)
                                                          cAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCATTCCACTCAAGTTGGCACTGGTAGa                                                                                            >  1:1115270/1‑150 (MQ=255)
                                                          cAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGa                                                                                            <  2:1115098/150‑1 (MQ=255)
                                                           aGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGaa                                                                                           <  1:2042704/150‑1 (MQ=255)
                                                            gTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAAGGTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGaaa                                                                                          >  2:1333228/1‑150 (MQ=255)
                                                            gTTCACCCGTCAGGGCATCGCGGATTTATGAAGGATTTAAACGCCCCCCCGTTTCTCCTGGCACAAACGGGAAATCTGCGCCAAATTGTCAGCGAACTTCGTCGAATGTTATACAATGTGTCAATCCACTAAAGTTGGCAGTGGTAGaaa                                                                                          >  2:491591/1‑150 (MQ=255)
                                                               cACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCACCCGTTTCACAGGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCCCTCAAGTTGGCACTGGTAGAAAcaa                                                                                       >  1:1790085/1‑148 (MQ=255)
                                                                   cGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCgg                                                                                   <  1:659088/150‑1 (MQ=255)
                                                                        gggCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTAc                                                                              <  1:742282/150‑1 (MQ=255)
                                                                        gggCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTAc                                                                              <  2:2032392/150‑1 (MQ=255)
                                                                         ggCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTAcc                                                                             <  2:1916341/150‑1 (MQ=255)
                                                                           cATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCAt                                                                           <  2:1830598/150‑1 (MQ=255)
                                                                           cATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCAt                                                                           >  1:1023816/1‑150 (MQ=255)
                                                                             tCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATaa                                                                         >  1:611852/1‑150 (MQ=255)
                                                                                 ggATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCt                                                                     >  1:1264496/1‑150 (MQ=255)
                                                                                    tttCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGc                                                                  <  2:638118/150‑1 (MQ=255)
                                                                                    tttCTGAAGGATTAAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACGGGTAGAAACCAGCGGTTTACCATAAGCCTGGc                                                                  <  2:137432/150‑1 (MQ=255)
                                                                                       cTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGcaca                                                               >  1:767395/1‑150 (MQ=255)
                                                                                       cTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGcaca                                                               <  2:1497604/150‑1 (MQ=255)
                                                                                       cTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAATTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGcaca                                                               >  2:394532/1‑150 (MQ=255)
                                                                                         gAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACaaa                                                             >  2:521796/1‑150 (MQ=255)
                                                                                         gAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACaaa                                                             >  2:244045/1‑150 (MQ=255)
                                                                                         gAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGAAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACaaa                                                             >  2:689198/1‑150 (MQ=255)
                                                                                          aaGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAg                                                            <  1:1759133/150‑1 (MQ=255)
                                                                                          aaGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAg                                                            >  2:550544/1‑150 (MQ=255)
                                                                                           aGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGc                                                           >  2:1883298/1‑150 (MQ=255)
                                                                                           aGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGc                                                           <  1:964406/150‑1 (MQ=255)
                                                                                                 tAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAaccac                                                     >  2:1877430/1‑150 (MQ=255)
                                                                                                 tAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAaccac                                                     >  1:1441312/1‑150 (MQ=255)
                                                                                                     cGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAAc                                                 <  1:887688/150‑1 (MQ=255)
                                                                                                         cccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAg                                                                        >  2:1881388/1‑123 (MQ=255)
                                                                                                         cccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAg                                                                        <  1:1881389/123‑1 (MQ=255)
                                                                                                           cccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATgg                                           <  2:1755523/150‑1 (MQ=255)
                                                                                                             cGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTc                                         <  1:712009/150‑1 (MQ=255)
                                                                                                              ggttCACCAGGCACAACCGGGAACGCCGCGCCAAATTGGGCGCGAACTTCGTCGACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCATGGCACAAAGCAACCACCAACGGATGGTCg                                        <  2:191542/148‑1 (MQ=255)
                                                                                                                ttCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGt                                      <  2:1441311/150‑1 (MQ=255)
                                                                                                                 tCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTg                                     <  1:1593981/150‑1 (MQ=255)
                                                                                                                 tCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAACCAACCACCAACGGATGGTCGGTg                                     >  2:1667401/1‑150 (MQ=255)
                                                                                                                  cACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTgc                                    >  1:1101033/1‑149 (MQ=255)
                                                                                                                  cACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTgc                                    >  1:1093976/1‑149 (MQ=255)
                                                                                                                  cACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGa                                    >  1:1666408/1‑150 (MQ=255)
                                                                                                                    ccAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACt                                  >  1:692579/1‑150 (MQ=255)
                                                                                                                     cAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAACCCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTc                                 >  2:76166/1‑150 (MQ=255)
                                                                                                                      aGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCg                                <  1:1999427/150‑1 (MQ=255)
                                                                                                                       ggCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGt                               >  1:1678590/1‑150 (MQ=255)
                                                                                                                       ggCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGt                               >  1:851317/1‑150 (MQ=255)
                                                                                                                       ggCACAACCGGCAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGt                               >  2:663288/1‑150 (MQ=255)
                                                                                                                          acaGCCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACa                            <  1:887729/150‑1 (MQ=255)
                                                                                                                           caACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTTCCATAAGCCTGGCCCAAAGCAACCACAAACGGATGGACGGTGCCTCGTACAg                           >  1:75861/1‑150 (MQ=255)
                                                                                                                           caACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAg                           <  2:1940232/150‑1 (MQ=255)
                                                                                                                                 ggAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGCACAGCAAGCg                     >  1:1750784/1‑150 (MQ=255)
                                                                                                                                 cgaagttcgcgCACAATTTGGCGCGAACTACGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGAACTGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGAAACGAACAGCAAGCGa                    >  1:1226178/8‑150 (MQ=255)
                                                                                                                                   aaCGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGaa                   <  1:981971/150‑1 (MQ=255)
                                                                                                                                    aCGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAAt                  <  1:1540285/150‑1 (MQ=255)
                                                                                                                                     gttcgcgCACAATTTGGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCTTAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATc                 >  1:631254/4‑150 (MQ=255)
                                                                                                                                     cgccgcGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATc                 <  1:1353494/150‑1 (MQ=255)
                                                                                                                                         gcgcCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAg             >  2:1142616/1‑150 (MQ=255)
                                                                                                                                           gcCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCg           <  1:1210391/150‑1 (MQ=255)
                                                                                                                                           gcCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCg           >  2:1191915/1‑150 (MQ=255)
                                                                                                                                           gcCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCg           <  2:819614/150‑1 (MQ=255)
                                                                                                                                            ccAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCgg          >  1:1586831/1‑150 (MQ=255)
                                                                                                                                              aaaTTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACCGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGcc        <  2:1338277/150‑1 (MQ=255)
                                                                                                                                                aTTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAAGCAAAGCGGCCCg      <  2:169731/150‑1 (MQ=255)
                                                                                                                                                    tgCGCGTACTTCGTCTACTGTTCGACAAGGTGGCAATCCCCTCAAGTTGGCACTGGTAGAAACCCGCGCTTTACCATAACCCTGGCACAAAGCAACCACCTAAGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCac  <  2:69060/150‑2 (MQ=255)
                                                                                                                                                    tgCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCaa  <  2:1614392/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGGGCTATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAA  >  NC_000913/3430800‑3431097

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