Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,405,586 C→T 8.1% P51S (CCA→TCA)  accB → biotin carboxyl carrier protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,405,5860CT 8.1% 301.0 / 10.9 99P51S (CCA→TCA) accBbiotin carboxyl carrier protein
Reads supporting (aligned to +/- strand):  ref base C (48/43);  new base T (5/3);  total (53/46)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGATATTCGTAAGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCGGACGCA  >  NC_000913/3405437‑3405732
                                                                                                                                                     |                                                                                                                                                  
tGGATATTCGTAAGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCAc                                                                                                                                                    >  2:898256/1‑150 (MQ=255)
      ttCGTAAGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAatga                                                                                                                                              >  2:1503238/1‑150 (MQ=255)
       tCGTAAGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAatgat                                                                                                                                             >  2:608706/1‑150 (MQ=255)
         gTAAGATTAAAAACCTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATgc                                                                                                                                           <  1:263367/150‑1 (MQ=255)
         gTAAGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATgc                                                                                                                                           <  2:483551/150‑1 (MQ=255)
         gTAAGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATgc                                                                                                                                           <  2:559695/150‑1 (MQ=255)
          tAAGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATgca                                                                                                                                          >  2:374648/1‑150 (MQ=255)
             caTTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTCTCCCTGTGATGCAACAAGCTTACGCTGCATCATTGGTgcagcg                                                                                                                                       <  1:246215/149‑2 (MQ=255)
              aTTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGcagcag                                                                                                                                      <  1:1351196/150‑1 (MQ=255)
               ttAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCATCATTGGTgcagcgta                                                                                                                                     <  1:1566902/150‑4 (MQ=255)
                 aaaaaaCTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGcagcca                                                                                                                                   >  1:648059/1‑150 (MQ=255)
                 aaaaaaCTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGcagcca                                                                                                                                   >  1:1617112/1‑150 (MQ=255)
                  aaaaaCTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATCCAGcagccac                                                                                                                                  >  1:6934/1‑149 (MQ=255)
                   aaaaCTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCATCATTGGTgcagcgtaagct                                                                                                                                 <  2:1157955/150‑8 (MQ=255)
                     aaCTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTc                                                                                                                               <  1:206513/150‑1 (MQ=255)
                       cTGATCTAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTATCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGATGCACCAATGATGCAGCAGCCAGCTCaa                                                                                                                             <  1:736265/150‑1 (MQ=255)
                        tGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAAt                                                                                                                            <  1:1393141/150‑1 (MQ=255)
                        tGATAGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAAt                                                                                                                            >  1:169242/1‑150 (MQ=255)
                         gATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATc                                                                                                                           >  2:371575/1‑150 (MQ=255)
                           tCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACTCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTa                                                                                                                         <  2:1522075/150‑1 (MQ=255)
                              aGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACg                                                                                                                      <  1:538670/150‑1 (MQ=255)
                              aGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACg                                                                                                                      <  1:1239210/150‑1 (MQ=255)
                                  ggTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGc                                                                                                                  <  2:1546457/150‑1 (MQ=255)
                                  ggTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGc                                                                                                                  >  2:873090/1‑150 (MQ=255)
                                      gaagaaTCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATgc                                                                                                                                           >  2:1644360/1‑121 (MQ=255)
                                      gaagaaTCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATgc                                                                                                                                           <  1:1644360/121‑1 (MQ=255)
                                      gaagaaTCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCt                                                                                                              >  1:203249/1‑150 (MQ=255)
                                      gaagaaTCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCt                                                                                                              >  1:573368/1‑150 (MQ=255)
                                        agaaTCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTcc                                                                                                            <  1:1012553/150‑1 (MQ=255)
                                         gaaTCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTccc                                                                                                           >  1:537574/1‑149 (MQ=255)
                                                gCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCg                                                                                                    <  1:1130225/150‑1 (MQ=255)
                                                         aaCTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCa                                                                                           <  1:17181/150‑1 (MQ=255)
                                                           cTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATg                                                                                         >  1:1610500/1‑150 (MQ=255)
                                                            tGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATgg                                                                                        <  1:898256/150‑1 (MQ=255)
                                                            tGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATgg                                                                                        >  1:185009/1‑150 (MQ=255)
                                                             ggAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGa                                                                                       >  1:234929/1‑150 (MQ=255)
                                                             ggAAATTTCTGAAGGCGAAGAGTCAGTACGCATTACCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTACCTTCCATGGa                                                                                       >  1:1126514/1‑150 (MQ=255)
                                                                   ttCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGcc                                                                                 <  1:1270499/150‑1 (MQ=255)
                                                                     cTGTAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCcag                                                                               >  2:448155/1‑150 (MQ=255)
                                                                      tGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCcagc                                                                              >  1:1146165/1‑150 (MQ=255)
                                                                      tGAAGGCGAAGAGTCAGTACGCATTAGACGTGCAGCTCCTGCCGCAAGTTTCACTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAAACTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCcagc                                                                              >  2:1133978/1‑150 (MQ=255)
                                                                       gAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCACCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAg                                                                                     >  1:943484/1‑142 (MQ=255)
                                                                       gAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCACCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAg                                                                                     <  2:943484/142‑1 (MQ=255)
                                                                        aaGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCcagcag                                                                            <  1:120509/150‑1 (MQ=255)
                                                                         cggCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCcagcagc                                                                           <  2:1323118/149‑1 (MQ=255)
                                                                               agagTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGa                                                                     <  2:1537097/150‑1 (MQ=255)
                                                                               agagTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGa                                                                     <  2:234929/150‑1 (MQ=255)
                                                                                gagTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGaa                                                                    <  1:198153/150‑1 (MQ=255)
                                                                                 agTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGaaa                                                                   >  1:817099/1‑150 (MQ=255)
                                                                                  gTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAAt                                                                  >  1:588998/1‑150 (MQ=255)
                                                                                  gTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAAt                                                                  >  1:554929/1‑150 (MQ=255)
                                                                                   tCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATc                                                                 <  1:1506126/150‑1 (MQ=255)
                                                                                         cGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGTAAATCAGTGGt                                                           >  1:1654128/1‑150 (MQ=255)
                                                                                          gCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTc                                                          <  2:1607617/150‑1 (MQ=255)
                                                                                           cATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTca                                                         >  1:149293/1‑150 (MQ=255)
                                                                                                gCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCg                                                    <  1:1485938/150‑1 (MQ=255)
                                                                                                gCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCg                                                    >  1:1487587/1‑150 (MQ=255)
                                                                                                gCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCg                                                    >  2:61935/1‑150 (MQ=255)
                                                                                                gCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCg                                                    >  1:432064/1‑150 (MQ=255)
                                                                                                 ccGTGCAGCTCCTGCCGCAAGTTTACCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGt                                                   >  2:471448/1‑150 (MQ=255)
                                                                                                  cGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTa                                                  >  1:492558/1‑150 (MQ=255)
                                                                                                  cGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTa                                                  >  1:349951/1‑150 (MQ=255)
                                                                                                    ggCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCGCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACg                                                <  2:698846/149‑1 (MQ=255)
                                                                                                      cAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGCAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGtt                                              >  1:1204740/1‑150 (MQ=255)
                                                                                                             tGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGAt                                       <  1:1503238/150‑1 (MQ=255)
                                                                                                              gCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGAtg                                      >  1:1592403/1‑150 (MQ=255)
                                                                                                               ccGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGAtgg                                     >  2:1546377/1‑150 (MQ=255)
                                                                                                                 gCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGAtggtt                                   >  2:136329/1‑150 (MQ=255)
                                                                                                                 gCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGAtggtt                                   <  1:11786/150‑1 (MQ=255)
                                                                                                                  cAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTGCGTTCCCCGAtggttg                                  >  2:360776/1‑150 (MQ=255)
                                                                                                                   aaGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGAtggttgg                                 >  2:1005713/1‑150 (MQ=255)
                                                                                                                     gTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTa                               <  1:471448/150‑1 (MQ=255)
                                                                                                                         cccTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACttt                           >  2:453322/1‑150 (MQ=255)
                                                                                                                          ccTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTc                          <  2:1319256/150‑1 (MQ=255)
                                                                                                                          ccTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTc                          >  1:930536/1‑150 (MQ=255)
                                                                                                                            tgtgATGCAACGAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTa                        >  1:266856/1‑150 (MQ=255)
                                                                                                                             gtgATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTAc                       <  2:641562/150‑1 (MQ=255)
                                                                                                                              tgATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTAcc                      <  1:608706/150‑1 (MQ=255)
                                                                                                                                 tGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCa                   >  1:1043458/1‑150 (MQ=255)
                                                                                                                                  gCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCAc                  >  2:701812/1‑150 (MQ=255)
                                                                                                                                  gCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCAc                  >  2:540317/1‑150 (MQ=255)
                                                                                                                                     acaaGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCAcccc               >  1:479120/1‑150 (MQ=255)
                                                                                                                                        aGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTAc                              >  1:1037422/1‑132 (MQ=255)
                                                                                                                                        aGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTAc                              <  2:1037422/132‑1 (MQ=255)
                                                                                                                                        aGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAg            >  1:284442/1‑150 (MQ=255)
                                                                                                                                        aGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAg            <  1:1133978/150‑1 (MQ=255)
                                                                                                                                         gCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGc           >  2:870426/1‑150 (MQ=255)
                                                                                                                                         ctggctgCTGCATCATTGGTGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGc           >  1:1235382/7‑150 (MQ=255)
                                                                                                                                           ttACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGccc         <  1:458290/150‑1 (MQ=255)
                                                                                                                                           ttACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGccc         <  1:448155/150‑1 (MQ=255)
                                                                                                                                            tACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCGCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCg        <  2:1119997/150‑1 (MQ=255)
                                                                                                                                            tACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCg        <  2:479120/150‑1 (MQ=255)
                                                                                                                                            tACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCg        <  1:162353/150‑1 (MQ=255)
                                                                                                                                             ctgCTGCATCATTGGTGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCgg       >  1:1545060/3‑150 (MQ=255)
                                                                                                                                             ctgCTGCATCATTGGTGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCgg       >  2:920542/3‑150 (MQ=255)
                                                                                                                                             ctgCTGCATCATTGGTGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCgg       >  1:309072/3‑150 (MQ=255)
                                                                                                                                             ctgCTGCATCATTGGTGCAGCAGCCAGCGCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAAGCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCCCCCCAAGCCCgg       >  1:247043/3‑150 (MQ=255)
                                                                                                                                               gCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACATTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCGGAc     <  1:1301483/150‑1 (MQ=255)
                                                                                                                                               gCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACATTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCGGAc     <  1:1295399/150‑1 (MQ=255)
                                                                                                                                                  gCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCGGACGCa  <  1:870426/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGGATATTCGTAAGATTAAAAAACTGATCGAGCTGGTTGAAGAATCAGGCATCTCCGAACTGGAAATTTCTGAAGGCGAAGAGTCAGTACGCATTAGCCGTGCAGCTCCTGCCGCAAGTTTCCCTGTGATGCAACAAGCTTACGCTGCACCAATGATGCAGCAGCCAGCTCAATCTAACGCAGCCGCTCCGGCGACCGTTCCTTCCATGGAAGCGCCAGCAGCAGCGGAAATCAGTGGTCACATCGTACGTTCCCCGATGGTTGGTACTTTCTACCGCACCCCAAGCCCGGACGCA  >  NC_000913/3405437‑3405732

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