Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 328,379 G→T 5.7% L119M (CTG→ATG) ‡ betB ← betaine aldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913328,3790GT 5.7% 299.5 / 4.0 87L119M (CTG→ATG) ‡betBbetaine aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (43/39);  new base T (2/3);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCG  >  NC_000913/328238‑328528
                                                                                                                                             |                                                                                                                                                     
ccGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCgt                                                                                                                                               <  2:899685/150‑1 (MQ=255)
 cGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGCTGATCCCGGCgct                                                                                                                                              <  1:1117930/150‑3 (MQ=255)
 cGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGCTGATCCCGGCgct                                                                                                                                              <  1:537162/150‑3 (MQ=255)
  ggggCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGCTGATCCCGGCgctg                                                                                                                                             <  2:968585/150‑4 (MQ=255)
  ggggCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCgtag                                                                                                                                             <  1:1197954/150‑1 (MQ=255)
   gggCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGGGCCGGGATCAGCCCGGCgtagt                                                                                                                                            <  2:996695/150‑1 (MQ=255)
     gCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTAc                                                                                                                                          <  1:1001896/150‑1 (MQ=255)
        gATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACGCAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTcc                                                                                                                                       >  1:1149032/1‑150 (MQ=255)
         aTTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCa                                                                                                                                      <  2:254460/150‑1 (MQ=255)
          tttCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAg                                                                                                                                     >  1:1477654/1‑150 (MQ=255)
            tCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCa                                                                                                                                   <  1:1170878/150‑1 (MQ=255)
                 aGGGCAATCTGGATCGGGTAGTGCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTcc                                                                                                                              <  2:116531/150‑1 (MQ=255)
                  gggCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCg                                                                                                                             >  2:272668/1‑150 (MQ=255)
                        tCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCgg                                                                                                                       <  1:64891/150‑1 (MQ=255)
                        tCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCgg                                                                                                                       <  1:709284/150‑1 (MQ=255)
                           ggATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTaa                                                                                                                    >  2:1313860/1‑150 (MQ=255)
                           ggATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTaa                                                                                                                    <  2:1481314/150‑1 (MQ=255)
                                 ggTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGatat                                                                                                              >  1:868504/1‑150 (MQ=255)
                                  gTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATc                                                                                                             <  2:869195/150‑1 (MQ=255)
                                          aTGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTg                                                                                                     <  2:1334153/150‑1 (MQ=255)
                                          aTGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTg                                                                                                     <  1:495494/150‑1 (MQ=255)
                                           tGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCCCGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGa                                                                                                    >  2:964158/1‑150 (MQ=255)
                                             cgcCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAgg                                                                                                  >  2:1655312/1‑150 (MQ=255)
                                               ccAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGtt                                                                                                >  1:1446043/1‑150 (MQ=255)
                                                  aTCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAATACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTcc                                                                                             >  2:275372/1‑150 (MQ=255)
                                                      cTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAat                                                                                         <  2:17129/150‑1 (MQ=255)
                                                          cACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGc                                                                                     >  1:80763/1‑150 (MQ=255)
                                                          cACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGc                                                                                     <  2:864914/150‑1 (MQ=255)
                                                                 cccAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCg                                                                              >  2:902244/1‑150 (MQ=255)
                                                                  ccAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCgg                                                                             <  1:1545839/150‑1 (MQ=255)
                                                                  ccAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCgg                                                                             <  2:444008/150‑1 (MQ=255)
                                                                  ccAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCgg                                                                             <  2:868504/150‑1 (MQ=255)
                                                                    aGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTTCCTTCCAGCGCCGGGATGAGGCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGTTTTCCGAATATGCTTTTCCGgtg                                                                           <  2:1401305/150‑1 (MQ=255)
                                                                     gCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGgtgt                                                                          <  1:1037073/150‑1 (MQ=255)
                                                                     gCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGgtgt                                                                          <  2:1000704/150‑1 (MQ=255)
                                                                           cgcgGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAggg                                                                    >  1:647086/1‑150 (MQ=255)
                                                                           cgcgGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAggg                                                                    >  2:1205591/1‑150 (MQ=255)
                                                                                  ggTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAg                                                             <  2:1106822/150‑1 (MQ=255)
                                                                                       acacAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTg                                                        >  1:283142/1‑150 (MQ=255)
                                                                                        cacaAAGGACGTTTCACGCAACGTGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAat                                                                                         <  1:801623/116‑1 (MQ=255)
                                                                                        cacaAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGc                                                       >  2:1567171/1‑150 (MQ=255)
                                                                                        cacaAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGc                                                       >  2:393678/1‑150 (MQ=255)
                                                                                        cacaAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTACCGGTGTCGAGGGTTTCCAGTTTTGc                                                       <  2:1477654/150‑1 (MQ=255)
                                                                                        cacaAAGGACGTTTCACGCAACGAGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAat                                                                                         >  2:801623/1‑116 (MQ=255)
                                                                                         acaAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCg                                                      <  1:1567171/150‑1 (MQ=255)
                                                                                         acaAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCg                                                      <  1:393678/150‑1 (MQ=255)
                                                                                         acaAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCg                                                      <  2:1315928/150‑1 (MQ=255)
                                                                                         acaAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCg                                                      <  2:257236/150‑1 (MQ=255)
                                                                                          caAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGa                                                     >  2:754822/1‑150 (MQ=255)
                                                                                           aaaGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAg                                                    <  1:503639/150‑1 (MQ=255)
                                                                                           aaaGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAg                                                    <  2:285126/150‑1 (MQ=255)
                                                                                           aaaGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAg                                                    <  2:1399679/150‑1 (MQ=255)
                                                                                             aGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGtt                                                  <  2:1446043/150‑1 (MQ=255)
                                                                                             aGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGtt                                                  <  1:871630/150‑1 (MQ=255)
                                                                                              ggACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTc                                                 <  2:1149032/150‑1 (MQ=255)
                                                                                                  gTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCa                                             >  1:394655/1‑150 (MQ=255)
                                                                                                     tCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTg                                          <  1:1205591/150‑1 (MQ=255)
                                                                                                        cGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGt                                       >  1:473789/1‑150 (MQ=255)
                                                                                                          cAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGttg                                     >  2:1620697/1‑149 (MQ=255)
                                                                                                          cAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTc                                     >  2:113068/1‑150 (MQ=255)
                                                                                                              gggATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGc                                 <  2:460126/150‑1 (MQ=255)
                                                                                                                gATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTACGGTGTCGAGGGTTTCAAGTTTTGAGAGTAGGTCATTGCGTTGACGCAg                               >  2:704362/1‑150 (MQ=255)
                                                                                                                 aTCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGa                              <  1:272668/150‑1 (MQ=255)
                                                                                                                  tCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGaa                             >  2:1027706/1‑150 (MQ=255)
                                                                                                                   cTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAat                            >  2:159173/1‑150 (MQ=255)
                                                                                                                   cTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTGGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAat                            <  2:365316/150‑1 (MQ=255)
                                                                                                                    tGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAata                           >  2:1405540/1‑150 (MQ=255)
                                                                                                                     ggCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAatat                          >  2:1271278/1‑150 (MQ=255)
                                                                                                                     ggCTGCCTTCCAGCGCCGGGATCAGCCCGGAGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAatat                          >  2:576056/1‑150 (MQ=255)
                                                                                                                      gCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATc                         >  2:407764/1‑150 (MQ=255)
                                                                                                                      gCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATc                         >  1:917384/1‑150 (MQ=255)
                                                                                                                         gCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAAc                      <  1:902244/150‑1 (MQ=255)
                                                                                                                          ccTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTGCAGTTTTGCGAGTTCGTCATTGCGTTCAAGCAGAATATCAACg                     >  1:924636/1‑150 (MQ=255)
                                                                                                                          ccTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACg                     >  2:222241/1‑150 (MQ=255)
                                                                                                                          ccTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACg                     >  1:978055/1‑150 (MQ=255)
                                                                                                                              actacgcCGGGCTGATCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGccc                 >  2:154606/5‑150 (MQ=255)
                                                                                                                               ctacgcCGGGCTGATCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCg                >  2:1246757/4‑150 (MQ=255)
                                                                                                                                 gcgcCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGCGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGAc              >  2:751168/1‑150 (MQ=255)
                                                                                                                                   gcCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGc            >  1:110895/1‑150 (MQ=255)
                                                                                                                                   gcCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGc            >  2:1253698/1‑150 (MQ=255)
                                                                                                                                     cGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAg          >  2:152206/1‑150 (MQ=255)
                                                                                                                                      gggATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGa         >  1:1043549/1‑150 (MQ=255)
                                                                                                                                      gggATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGa         >  2:43880/1‑150 (MQ=255)
                                                                                                                                       ggATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGaa        >  1:1262974/1‑150 (MQ=255)
                                                                                                                                           cAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATAcg    >  2:1346416/1‑150 (MQ=255)
                                                                                                                                            aGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATAcgc   >  2:1083716/1‑150 (MQ=255)
                                                                                                                                            aGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATAcgc   >  1:429010/1‑150 (MQ=255)
                                                                                                                                             gCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATAcgcg  >  2:421626/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
CCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCG  >  NC_000913/328238‑328528

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: ‑