Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,796,101 C→A 100% D243Y (GAC→TAC)  thrS ← threonyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,796,1010CA98.5% 615.1 / 0.6 201D243Y (GAC→TAC) thrSthreonyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (97/101);  total (98/103)
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.17e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 206 total aligned reads displayed.

TGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATG  >  CP009273/1795956‑1796245
                                                                                                                                                 |                                                                                                                                                
  aTACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGt                                                                                                                                              <  1:1402834/148‑1 (MQ=255)
   tACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGtt                                                                                                                                             >  2:2143004/1‑148 (MQ=255)
   tACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGtt                                                                                                                                             >  2:2144906/1‑148 (MQ=255)
   tACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGtt                                                                                                                                             <  2:1216124/148‑1 (MQ=255)
    aCTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGttt                                                                                                                                            >  2:2028455/1‑148 (MQ=255)
    aCTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGttt                                                                                                                                            <  1:161193/148‑1 (MQ=255)
     cTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTc                                                                                                                                           >  2:315805/1‑148 (MQ=255)
      tGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTcg                                                                                                                                          >  2:211680/1‑148 (MQ=255)
       ggTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTcgc                                                                                                                                         >  1:587589/1‑148 (MQ=255)
        gTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCAAATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTcgcg                                                                                                                                        >  1:3136910/1‑148 (MQ=255)
         tACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCgg                                                                                                                                       >  2:3795041/1‑148 (MQ=255)
          aCTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGc                                                                                                                                      >  1:545507/1‑148 (MQ=255)
          aCTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTAAGGTGGTAGCGTTTCGCGGc                                                                                                                                      >  1:3214783/1‑148 (MQ=255)
            tctTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTACAGCATATTGTACAGGTCGAGCTGTTTAACGATTTTACGGTGGTAGCGTTTAGCGGctt                                                                                                                                    >  1:692071/1‑148 (MQ=255)
             ctTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGcttc                                                                                                                                   <  2:105060/148‑1 (MQ=255)
               ttCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGcttctt                                                                                                                                 <  2:3054138/148‑1 (MQ=255)
                  aGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCa                                                                                                                              <  2:12810/148‑1 (MQ=255)
                           cgaacAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCa                                                                                                                     <  2:1757590/148‑1 (MQ=255)
                             aacaaacaCTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAgg                                                                                                                   >  2:2235287/1‑148 (MQ=255)
                               caaacaCTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTa                                                                                                                 <  1:1256580/148‑1 (MQ=255)
                               caaacaCTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTAAGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTa                                                                                                                 >  1:3738224/1‑148 (MQ=255)
                                  acacTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGc                                                                                                              <  2:1567439/148‑1 (MQ=255)
                                  acacTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGc                                                                                                              >  2:3232711/1‑148 (MQ=255)
                                  acacTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGc                                                                                                              <  1:1969922/148‑1 (MQ=255)
                                  acacTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGc                                                                                                              >  2:2528734/1‑148 (MQ=255)
                                  acacTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGc                                                                                                              <  1:2858888/148‑1 (MQ=255)
                                  acacTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGc                                                                                                              <  2:151763/148‑1 (MQ=255)
                                    acTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGt                                                                                                            <  1:1481554/148‑1 (MQ=255)
                                       tCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTaa                                                                                                         >  2:3430253/1‑148 (MQ=255)
                                         cAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGt                                                                                                       >  2:1978129/1‑148 (MQ=255)
                                             tCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCtt                                                                                                   <  2:3770825/148‑1 (MQ=255)
                                             tCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCtt                                                                                                   <  2:1831048/148‑1 (MQ=255)
                                              cACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCttt                                                                                                  <  2:278962/148‑1 (MQ=255)
                                               aCGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCtttt                                                                                                 <  2:1968191/148‑1 (MQ=255)
                                                cGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCttttt                                                                                                >  1:3496089/1‑148 (MQ=255)
                                                  gAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTg                                                                                              <  2:1612074/148‑1 (MQ=255)
                                                   aaGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATTGTACAGGTCGAGCTGTTTTCCGATTTTAAGGTGGTATCGTTTGGCGGCTTGTTCCAGGCGCTGCAGGTAAGCGGTAAGTGGtttttttt                                                                                             >  1:3525297/1‑146 (MQ=255)
                                                    aGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTc                                                                                            <  2:177083/148‑1 (MQ=255)
                                                      aTGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTg                                                                                          >  2:1892803/1‑148 (MQ=255)
                                                      aTGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTg                                                                                          >  2:3786903/1‑148 (MQ=255)
                                                          tCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCa                                                                                      <  2:160757/148‑1 (MQ=255)
                                                                 gtcgtTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGccgta                                                                               <  1:2028072/148‑1 (MQ=255)
                                                                  tcgtTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGccgtac                                                                              >  2:854553/1‑148 (MQ=255)
                                                                  tcgtTGTGCCAGAATACCATACCCGGCGCTTATTACTGCATATGGTACAGGTCGAGCTGTTTAACCATTTTAAGGTGGTATCGTTTCGCGGCTTGTTCCAGGCGCGGCCGGTAATCGTTAAGTGGTTTTTTGTATTCCCAcccccgac                                                                              >  1:1420572/1‑144 (MQ=255)
                                                                   cgtTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGccgtacc                                                                             <  1:1364154/148‑1 (MQ=255)
                                                                     ttGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGccgtaccgt                                                                           <  1:2364519/148‑1 (MQ=255)
                                                                     ttGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGGTAAGTGCTTTTTTGTCTGCCCACGccgtaccgt                                                                           >  2:665004/1‑148 (MQ=255)
                                                                      tgtgCCAGAATACCATACCCGGCGCTTCTTCCTGCCTATGGTACAGGTCGAGCTGTTTAACGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGGTAAGTGCTTTTTTGTGTGCCCACGccgtaccgga                                                                          >  1:177167/1‑146 (MQ=255)
                                                                      tgtgCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGccgtaccgta                                                                          >  2:1106166/1‑148 (MQ=255)
                                                                       gtgCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTaa                                                                         <  2:116596/148‑1 (MQ=255)
                                                                        tgCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTaaa                                                                        <  2:2863143/148‑1 (MQ=255)
                                                                         gCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAAt                                                                       <  2:1480835/148‑1 (MQ=255)
                                                                          ccAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATa                                                                      >  2:2949860/1‑148 (MQ=255)
                                                                            aGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACg                                                                    <  1:2157078/148‑1 (MQ=255)
                                                                             gAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGt                                                                   >  2:1476334/1‑148 (MQ=255)
                                                                             gAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGt                                                                   >  2:2107372/1‑148 (MQ=255)
                                                                                taccataccCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGtttc                                                                >  2:1024079/1‑146 (MQ=255)
                                                                                taccataccCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTAACGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGGTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGc                                                                >  2:2933726/1‑148 (MQ=255)
                                                                                 accataccCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCa                                                               <  2:3163531/148‑1 (MQ=255)
                                                                                   cataccCGGCGCTTCTTCGAGCCTATGGTAAAGGTCGTGCTGTTTACCCATATTACCGTGGTAGCGTTTCGCTGCATCCTCCAGGCGCAGCAGGTAATCGTTAAGTGCTTTTTTGTCTGCCCACTCCGGACCGTAACTAAGTTTCAAc                                                             >  2:3634301/1‑148 (MQ=255)
                                                                                       cccGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACAttt                                                         <  2:3284254/148‑1 (MQ=255)
                                                                                         cGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTttg                                                       <  1:2918015/148‑1 (MQ=255)
                                                                                          ggCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTttgt                                                      >  2:3040574/1‑148 (MQ=255)
                                                                                            cgcTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTttgttg                                                    >  1:3277396/1‑148 (MQ=255)
                                                                                            cgcTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTttgttg                                                    <  2:1924167/148‑1 (MQ=255)
                                                                                                tcttcCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGc                                                >  2:1079217/1‑148 (MQ=255)
                                                                                                tcttcCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGc                                                >  2:432130/1‑148 (MQ=255)
                                                                                                 cttcCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCt                                               >  1:2643062/1‑148 (MQ=255)
                                                                                                   tcCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGt                                             <  1:3795041/148‑1 (MQ=255)
                                                                                                   tcCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGt                                             <  1:1978258/148‑1 (MQ=255)
                                                                                                   tcCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGt                                             <  1:1922436/148‑1 (MQ=255)
                                                                                                      tGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTcgc                                          >  1:105010/1‑148 (MQ=255)
                                                                                                       gCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTcgcc                                         <  1:1160340/148‑1 (MQ=255)
                                                                                                       gCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTcgcc                                         <  1:2689434/148‑1 (MQ=255)
                                                                                                        cATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTcgcca                                        <  2:3767927/148‑1 (MQ=255)
                                                                                                        cATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTcgcca                                        >  2:3635666/1‑148 (MQ=255)
                                                                                                          tatGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTATCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTcgccacg                                      >  2:259823/1‑148 (MQ=255)
                                                                                                          tatGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTcgccacg                                      >  1:1205670/1‑148 (MQ=255)
                                                                                                          tatGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTcgccacg                                      >  1:1864061/1‑148 (MQ=255)
                                                                                                                  aGGTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAg                              <  1:3531816/148‑1 (MQ=255)
                                                                                                                   ggTCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGc                             >  1:2735129/1‑148 (MQ=255)
                                                                                                                     tCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGccc                           <  2:2680115/148‑1 (MQ=255)
                                                                                                                     tCGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGccc                           <  2:545507/148‑1 (MQ=255)
                                                                                                                      cGAGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGcccc                          >  2:2002225/1‑148 (MQ=255)
                                                                                                                        aGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTg                        <  1:3430253/148‑1 (MQ=255)
                                                                                                                        aGCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTg                        >  1:2311321/1‑148 (MQ=255)
                                                                                                                         gCTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGc                       <  2:1937722/148‑1 (MQ=255)
                                                                                                                          cTGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGcc                      >  2:1889424/1‑148 (MQ=255)
                                                                                                                           tGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCg                     <  1:87262/148‑1 (MQ=255)
                                                                                                                           tGTTTACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCg                     >  1:2395197/1‑148 (MQ=255)
                                                                                                                              ttACCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGttt                  >  1:2498093/1‑148 (MQ=255)
                                                                                                                                aCCGATTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTc                <  1:854553/148‑1 (MQ=255)
                                                                                                                                    aTTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTa            <  1:3033775/148‑1 (MQ=255)
                                                                                                                                    aTTTTACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGACGttttattca            >  2:754816/1‑143 (MQ=255)
                                                                                                                                     ttttACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAg           <  1:2243707/148‑1 (MQ=255)
                                                                                                                                     ttttACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAg           >  1:1511470/1‑148 (MQ=255)
                                                                                                                                      tttACGGTGGTAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGt          <  2:3498013/148‑1 (MQ=255)
                                                                                                                                            gtggtAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGaaa    <  1:2235287/148‑1 (MQ=255)
                                                                                                                                            gtggtAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGaaa    <  2:178834/148‑1 (MQ=255)
                                                                                                                                            gtggtAGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGaaa    <  2:177808/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATG  >  CP009273/1795956‑1796245

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