Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,708,413 G→T 13.3% P448T (CCG→ACG)  bisC ← biotin sulfoxide reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,708,4130GT13.3% 929.1 / 75.1 275P448T (CCG→ACG) bisCbiotin sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base G (116/122);  new base T (15/22);  total (131/144)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 277 total aligned reads displayed.

ACTTGCTTCATCGGCACCAGATGCTGATTACTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACA  >  CP009273/3708265‑3708561
                                                                                                                                                    |                                                                                                                                                    
  ttGCTTCATCGGCACCAGATGCTGATTACTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGttt                                                                                                                                                   <  1:3255072/150‑1 (MQ=255)
  ttGCTTCATCGGCACCAGATGCTGATTACTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGttt                                                                                                                                                   <  1:3614676/150‑1 (MQ=255)
   tGCTTCATCGGCACCAGATGCTGATTACTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGTtttt                                                                                                                                                  <  2:265351/150‑1 (MQ=255)
   tGCTTCATCGGCACCAGATGCTGATTACTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGtttt                                                                                                                                                  >  2:1933104/1‑150 (MQ=255)
    gCTTCATCGGCACCAGATGCTGATTACTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGttttt                                                                                                                                                 <  2:581171/150‑1 (MQ=255)
      ttCATCGGCACCAGATGCTGATTACTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGTTTTTTGc                                                                                                                                               >  2:2443974/1‑150 (MQ=255)
        cATCGGCACCAGATGCTGATTACTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCa                                                                                                                                             <  2:963038/150‑1 (MQ=255)
          tCGGCACCAGATGCTGATTACTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGTTTTTTGCCAgg                                                                                                                                           >  1:1111300/1‑150 (MQ=255)
          tCGGCACCAGATGCTGATTACTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAgg                                                                                                                                           <  1:2339642/150‑1 (MQ=255)
            ggCACCAGATGCTGATTACTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCa                                                                                                                                         >  2:800616/1‑150 (MQ=255)
             gCACCAGATGCTGATTACTGTACTCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGAGATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCAc                                                                                                                                        <  2:1629098/150‑1 (MQ=255)
                 cAGATGCTGATTACTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGAt                                                                                                                                    >  2:1848953/1‑150 (MQ=255)
                 cAGATGCTGATTACTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCAGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGAt                                                                                                                                    >  2:1848869/1‑150 (MQ=255)
                     tGCTGATTACTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAgg                                                                                                                                <  1:3784305/150‑1 (MQ=255)
                        tGATTACTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGTTTTTTGCCAGGCACGGATCAGGCGa                                                                                                                             <  1:4361706/150‑1 (MQ=255)
                         gATTACTGTAATCACCGGTCATGGTGAGATCATGACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGAt                                                                                                                            <  2:2866728/150‑1 (MQ=255)
                             aCTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGTTTTTTGCCAGGCACGGATCAGGCGATTTGt                                                                                                                        <  2:453708/150‑1 (MQ=255)
                             aCTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGt                                                                                                                        >  2:58755/1‑150 (MQ=255)
                               tGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTAt                                                                                                                      >  1:1884268/1‑150 (MQ=255)
                                  aaTCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGTTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCt                                                                                                                   <  2:345660/150‑1 (MQ=255)
                                  aaTCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCt                                                                                                                   <  1:565633/150‑1 (MQ=255)
                                      aCCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGTTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCtgatg                                                                                                               <  2:3223166/150‑1 (MQ=255)
                                      aCCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCtgatg                                                                                                               <  2:348490/150‑1 (MQ=255)
                                       ccGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCtgatga                                                                                                              <  2:2414871/150‑1 (MQ=255)
                                         ggTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCtgatgatg                                                                                                            <  1:975905/150‑1 (MQ=255)
                                          gTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCtgatgatga                                                                                                           <  1:1780556/150‑1 (MQ=255)
                                            cATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGt                                                                                                         <  1:1035733/150‑1 (MQ=255)
                                             aTGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTa                                                                                                        <  2:1879114/150‑1 (MQ=255)
                                                gTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGTTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAg                                                                                                     <  2:2596622/150‑1 (MQ=255)
                                                  gagaTCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGtt                                                                                                   <  2:2916155/150‑1 (MQ=255)
                                                  gagaTCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGtt                                                                                                   >  1:846550/1‑150 (MQ=255)
                                                  gagaTCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGtt                                                                                                   <  1:3122009/150‑1 (MQ=255)
                                                       cATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCAc                                                                                              >  2:2177032/1‑150 (MQ=255)
                                                       cATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCAc                                                                                              >  2:2672806/1‑150 (MQ=255)
                                                          tACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCAccgc                                                                                           >  2:3973297/1‑150 (MQ=255)
                                                            cGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGTTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCAccgccc                                                                                         <  1:2539213/150‑1 (MQ=255)
                                                              cTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGTTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCAccgcccgc                                                                                       >  2:792552/1‑150 (MQ=255)
                                                               tCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACcgcccgcc                                                                                      <  2:4061333/150‑1 (MQ=255)
                                                               tCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACcgcccgcc                                                                                      <  1:892973/150‑1 (MQ=255)
                                                                 aaaaGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCcca                                                                                    >  2:895498/1‑150 (MQ=255)
                                                                 aaaaGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCcca                                                                                    <  1:598690/150‑1 (MQ=255)
                                                                   aaGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGTTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCccacc                                                                                  >  1:1053120/1‑150 (MQ=255)
                                                                   aaGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCccacc                                                                                  >  1:3997036/1‑150 (MQ=255)
                                                                    agagGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCccacca                                                                                 >  2:1143984/1‑150 (MQ=255)
                                                                      agGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGa                                                                               >  2:3761738/1‑150 (MQ=255)
                                                                        gTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATa                                                                             <  1:3195969/150‑1 (MQ=255)
                                                                         tAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATaa                                                                            <  2:3023180/150‑1 (MQ=255)
                                                                          aGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATAACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATaaa                                                                           >  1:4250447/1‑150 (MQ=255)
                                                                           gTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATaaaa                                                                          <  1:221122/150‑1 (MQ=255)
                                                                             cGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACg                                                                        >  1:355666/1‑150 (MQ=255)
                                                                             cGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACg                                                                        >  1:355870/1‑150 (MQ=255)
                                                                               caggcagAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGaa                                                                      <  2:683205/150‑1 (MQ=255)
                                                                                aggcagAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAat                                                                     >  2:2607758/1‑150 (MQ=255)
                                                                                 ggcagAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAata                                                                    >  1:3834876/1‑150 (MQ=255)
                                                                                  gcagAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAatat                                                                   <  2:2707475/150‑1 (MQ=255)
                                                                                   cagAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATc                                                                  >  2:1263068/1‑150 (MQ=255)
                                                                                   cagAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATc                                                                  <  1:4608323/150‑1 (MQ=255)
                                                                                     gAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCg                                                                >  1:2316518/1‑150 (MQ=255)
                                                                                       aCGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCggg                                                              <  2:440242/150‑1 (MQ=255)
                                                                                       aCGATATCCGAGTGTTTTGCCGCCGCCGTCCATAAGCATTCAGAGATCTCCACCTGCTTCGTTTTTTGCCATGAACGGATCAGGCGATTGTTTTCCTGATGATGAGTAAAGTTGTCACCGCCCGACCAACTGATAAAAATAATATCCggg                                                              <  2:3749466/150‑1 (MQ=255)
                                                                                          atatCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGaaa                                                           >  1:4400109/1‑150 (MQ=255)
                                                                                          atatCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGaaa                                                           >  2:2033997/1‑150 (MQ=255)
                                                                                          atatCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGaaa                                                           >  1:140540/1‑150 (MQ=255)
                                                                                            atCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATg                                                         >  1:3149575/1‑150 (MQ=255)
                                                                                             tCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCTGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGt                                                        <  2:1434845/150‑1 (MQ=255)
                                                                                              ccGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTc                                                       >  2:2257103/1‑150 (MQ=255)
                                                                                              ccGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTc                                                       <  2:4369188/150‑1 (MQ=255)
                                                                                               cgcgTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCg                                                      <  2:1392585/150‑1 (MQ=255)
                                                                                               cgcgTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCg                                                      >  1:3487907/1‑150 (MQ=255)
                                                                                                gcgTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCgg                                                     >  2:3277197/1‑150 (MQ=255)
                                                                                                    gtTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCa                                                 <  1:1532846/150‑1 (MQ=255)
                                                                                                       ttGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATAc                                              <  1:2017688/150‑1 (MQ=255)
                                                                                                       ttGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATAc                                              <  2:1095020/150‑1 (MQ=255)
                                                                                                         gccgccgccGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGTTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCg                                            >  2:2050102/1‑150 (MQ=255)
                                                                                                         gccgccgccGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCg                                            <  2:2887219/150‑1 (MQ=255)
                                                                                                           cgccgccgTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGtt                                          >  2:2087475/1‑150 (MQ=255)
                                                                                                            gccgccgTCCAGTAGCATTCAGCGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTTGTTTCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTtg                                         <  1:1898570/150‑1 (MQ=255)
                                                                                                                 cgTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTg                                    >  1:2892230/1‑150 (MQ=255)
                                                                                                                  gTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGa                                   <  1:2642319/150‑1 (MQ=255)
                                                                                                                   tCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGat                                  <  1:4529664/150‑1 (MQ=255)
                                                                                                                     cAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGatat                                <  2:4033159/150‑1 (MQ=255)
                                                                                                                      aGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATg                               <  1:3719612/150‑1 (MQ=255)
                                                                                                                       gAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATgc                              <  1:3077255/150‑1 (MQ=255)
                                                                                                                          gCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGcc                           <  1:4489846/150‑1 (MQ=255)
                                                                                                                           cATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGcca                          <  1:1125859/150‑1 (MQ=255)
                                                                                                                            aTTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGccac                         >  2:4420433/1‑150 (MQ=255)
                                                                                                                             ttCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGccacc                        >  1:3397014/1‑150 (MQ=255)
                                                                                                                              tCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGccacca                       <  1:502837/150‑1 (MQ=255)
                                                                                                                               cAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAg                      >  1:2767816/1‑150 (MQ=255)
                                                                                                                                agagATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAgg                     >  1:784358/1‑150 (MQ=255)
                                                                                                                                 gagaTCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCATATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAggg                    <  2:2339872/150‑1 (MQ=255)
                                                                                                                                  agaTCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGt                   >  2:777678/1‑150 (MQ=255)
                                                                                                                                      caccacCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTc               <  1:4122619/150‑1 (MQ=255)
                                                                                                                                       accaccAGCTCCGTTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTcc              <  2:721605/150‑1 (MQ=255)
                                                                                                                                       accaccAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTcc              >  1:2674896/1‑150 (MQ=255)
                                                                                                                                       accaccAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTcc              <  1:711406/150‑1 (MQ=255)
                                                                                                                                                tCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTca     <  2:4601511/150‑1 (MQ=255)
                                                                                                                                                  cGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTTAAGTTGGCACCGCCCGCCCAACAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTcaac   <  2:3849022/150‑1 (MQ=255)
                                                                                                                                                  cGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTcaac   <  1:1930263/150‑1 (MQ=255)
                                                                                                                                                  cGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTcaac   <  1:918150/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACTTGCTTCATCGGCACCAGATGCTGATTACTGTAATCACCGGTCATGGTGAGATCATTACGCTCAAAAGAGGTAGTCGCAGGCAGAACGATATCCGCGTGTTTTGCCGCCGCCGTCCAGAAGCATTCAGAGATCACCACCAGCTCCGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACA  >  CP009273/3708265‑3708561

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