Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 32,911 G→T 5.2% E699* (GAA→TAA) ‡ carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091332,9110GT 5.2% 462.7 / 11.1 153E699* (GAA→TAA) ‡carBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base G (74/71);  new base T (3/5);  total (77/76)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

Only 100 of 155 total aligned reads displayed.

AACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCG  >  NC_000913/32762‑33060
                                                                                                                                                     |                                                                                                                                                     
  cTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGAGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGaa                                                                                                                                                     >  2:803041/1‑150 (MQ=255)
   tGGCGCGCGCGCTGGAAGCTGCTGGCGTACCTGTTATAGGCACCATCCCGGATGATATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGaaa                                                                                                                                                    <  1:156948/150‑1 (MQ=255)
       gcgcgcgcTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGc                                                                                                                                                <  2:1511760/150‑1 (MQ=255)
        ggcgcgcGGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCg                                                                                                                                               <  1:143195/149‑1 (MQ=255)
            cgcTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTag                                                                                                                                           >  2:869375/1‑150 (MQ=255)
            cgcTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTag                                                                                                                                           <  1:1866027/150‑1 (MQ=255)
            cgcTGGAAGCTGCTGGCGTACCGGTTATCGCCACCAGCCCGGATGCTATCGACCGGGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCCATTTCAATAGCGGTaa                                                                                                                                           <  1:77628/150‑2 (MQ=255)
              gtGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTagag                                                                                                                                         <  2:543693/149‑1 (MQ=255)
              cTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTagag                                                                                                                                         >  1:1022400/1‑150 (MQ=255)
              cTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTagag                                                                                                                                         >  2:529355/1‑150 (MQ=255)
              cTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTagag                                                                                                                                         >  1:896321/1‑150 (MQ=255)
              cTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTagag                                                                                                                                         <  1:66436/150‑1 (MQ=255)
                ggAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGaa                                                                                                                                       >  1:1278199/1‑150 (MQ=255)
                ggAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCCATTTCAATAGCGGTaacggt                                                                                                                                       <  1:819562/150‑6 (MQ=255)
                  aaGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGGTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCCATTTCAATAGCGGTaacggtgg                                                                                                                                     <  1:485407/150‑8 (MQ=255)
                     cggctgGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCCGAAGACCGGGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGTAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGa                                                                                                                                  <  2:64106/148‑1 (MQ=255)
                       gctgGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGCATCAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGaaa                                                                                                                                >  1:417994/1‑150 (MQ=255)
                       gctgGCGTACCGGTAATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGaaa                                                                                                                                <  1:442584/150‑1 (MQ=255)
                        ctgGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAag                                                                                                                               <  1:1746223/150‑1 (MQ=255)
                          ggCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAagag                                                                                                                             >  1:1706309/1‑150 (MQ=255)
                            cGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGAGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGCGCGTCGGAACCTGAAACAACCGGCGCACGCCCCCGTGACCGCTATTGAAATGGCGGTAGAGACGGCGATAGAGAt                                                                                                                           >  1:298106/1‑150 (MQ=255)
                              tACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTg                                                                                                                         <  1:1306051/150‑1 (MQ=255)
                              tACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTg                                                                                                                         <  1:1529183/150‑1 (MQ=255)
                              tACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTg                                                                                                                         <  2:1305914/150‑1 (MQ=255)
                                ccGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCGTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGc                                                                                                                       <  2:171462/150‑1 (MQ=255)
                                 cGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCt                                                                                                                      <  2:903822/150‑1 (MQ=255)
                                 cGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCt                                                                                                                      <  2:413928/150‑1 (MQ=255)
                                  ggTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTTTTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTa                                                                                                                     >  1:1623689/1‑150 (MQ=255)
                                      aTCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGGACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACCACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCg                                                                                                                 >  1:1158837/1‑150 (MQ=255)
                                         ggCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCtg                                                                                                              >  2:1205512/1‑150 (MQ=255)
                                             ccAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCtggtgg                                                                                                          >  1:1077578/1‑150 (MQ=255)
                                               aGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTa                                                                                                        <  1:1413232/150‑1 (MQ=255)
                                                  ccGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTAcgt                                                                                                     <  2:478531/150‑1 (MQ=255)
                                                     gATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTAcgtccg                                                                                                  <  1:1580043/150‑1 (MQ=255)
                                                         cTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCtt                                                                                              <  2:118034/150‑1 (MQ=255)
                                                          tATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGACACACCCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTa                                                                                             >  2:513034/1‑150 (MQ=255)
                                                           aTCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTAc                                                                                            >  1:612710/1‑150 (MQ=255)
                                                            tCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTAATGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACg                                                                                           <  2:1648406/150‑1 (MQ=255)
                                                              gACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGtt                                                                                         >  2:1790483/1‑150 (MQ=255)
                                                              gACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGtt                                                                                         >  2:1571902/1‑150 (MQ=255)
                                                                  gTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcg                                                                                     >  2:1641738/1‑150 (MQ=255)
                                                                   tGCAGAAGACCGTGAACGCTTCCTGCATGCGGTTGCGCGTCTGAAACTGAAACAACCGGCGAACGCCCCCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGATACCCGCTGGGGTTACGTCCGTCTTACGTTCTcgg                                                                                    >  2:1127326/1‑150 (MQ=255)
                                                                    gCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcggc                                                                                   >  2:695462/1‑150 (MQ=255)
                                                                     cAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcggcg                                                                                  <  2:339244/150‑1 (MQ=255)
                                                                      agaagaCCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcggcgg                                                                                 <  1:1802106/150‑1 (MQ=255)
                                                                         cgaCCGGGAACGCTTCCAGCATGCGGTGGATCGTCTGAACCTGAAACACCCGGCGACCGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGTTACGTCCGTCTTACGTTCTCGGCGGTCg                                                                              <  2:1562168/149‑1 (MQ=255)
                                                                         agaCCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCg                                                                              <  2:1300713/150‑1 (MQ=255)
                                                                         agaCCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCg                                                                              <  2:610097/150‑1 (MQ=255)
                                                                         agaCCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCGATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCg                                                                              >  1:952352/1‑150 (MQ=255)
                                                                             cGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCg                                                                          <  2:1205569/150‑1 (MQ=255)
                                                                                gAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATg                                                                       <  2:1706309/150‑1 (MQ=255)
                                                                                gAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATg                                                                       <  2:896321/150‑1 (MQ=255)
                                                                                gAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATg                                                                       <  1:529355/150‑1 (MQ=255)
                                                                                       tCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCg                                                                >  1:1464028/1‑150 (MQ=255)
                                                                                         cAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTc                                                              <  2:746099/150‑1 (MQ=255)
                                                                                          aGCATGCGGGTGAGCGTCTGAAACTGAAACTACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTCGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCt                                                             <  2:266562/150‑1 (MQ=255)
                                                                                           gCATGCGGTTGAGAGTCTGAAACTGAAACAACCGGCGAACGCCACCGTAACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTAGCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTa                                                            <  1:617191/150‑1 (MQ=255)
                                                                                            cATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTAt                                                           >  2:1729873/1‑150 (MQ=255)
                                                                                            cATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTAt                                                           >  2:413668/1‑150 (MQ=255)
                                                                                               gCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGAc                                                        >  2:504012/1‑150 (MQ=255)
                                                                                                 ggTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGa                                                      >  2:203112/1‑150 (MQ=255)
                                                                                                    tGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGc                                                   <  1:1409713/150‑1 (MQ=255)
                                                                                                     gAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGCAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCt                                                  >  1:491614/1‑150 (MQ=255)
                                                                                                      aGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTg                                                 >  1:1046929/1‑150 (MQ=255)
                                                                                                            tgaaactgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGc                                           >  1:1304197/1‑150 (MQ=255)
                                                                                                             gaaactgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCg                                          >  2:1049416/1‑150 (MQ=255)
                                                                                                               atctgaaacACCCGGCGAAGGCCACCGTTACCGCTAGTGTAATGGCGGTAGTGAAGGCGTAAGTGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAATTCGTCTATGACGATGCTGACCTGCGTc                                        <  1:1247066/148‑1 (MQ=255)
                                                                                                               aactgaaacAACCGGCGATCGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTc                                        <  2:1022400/150‑1 (MQ=255)
                                                                                                               aactgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTc                                        >  1:725272/1‑150 (MQ=255)
                                                                                                                actgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCg                                       >  1:1181829/1‑150 (MQ=255)
                                                                                                                 ctgtaacAACCGGCGAACGCCTCCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGTTAGAGATTGGATACCCGCTGGTGGTACTTACGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGc                                      <  2:400120/150‑1 (MQ=255)
                                                                                                                 ctgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGc                                      >  2:702308/1‑150 (MQ=255)
                                                                                                                   gaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTa                                    >  2:124298/1‑150 (MQ=255)
                                                                                                                    aaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTAc                                   >  2:1473423/1‑150 (MQ=255)
                                                                                                                    aaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACACGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATAGTCTATGTCGAAGCTGACCTGCGTCGCTAc                                   >  2:746637/1‑150 (MQ=255)
                                                                                                                       caacCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTc                                >  2:1461282/1‑150 (MQ=255)
                                                                                                                       caacCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGGGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTc                                <  1:857936/150‑1 (MQ=255)
                                                                                                                        aacCGGCGTACGCCACCGTTACCGCTATTGAATTGGCGGTAGAGACTGCGAAAGAGATTGGCTCCCCGCTGGTGGTACGTCCGTCTTACTGTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTcc                               <  2:277171/150‑1 (MQ=255)
                                                                                                                         acCGGCGAACGCCACTGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCa                              >  2:1728177/1‑150 (MQ=255)
                                                                                                                          ccGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAg                             >  2:1592996/1‑150 (MQ=255)
                                                                                                                          ccGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGAAGAAGCAGACCTGCGTCGCTACTTCCAg                             >  2:304967/1‑150 (MQ=255)
                                                                                                                           cGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGa                            >  1:243599/1‑150 (MQ=255)
                                                                                                                              cGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGAcgg                         >  1:521232/1‑150 (MQ=255)
                                                                                                                                  cGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGt                     <  2:1464028/150‑1 (MQ=255)
                                                                                                                                    ccttctctACCGCCATTTCAATAGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCa                   >  1:1139645/8‑150 (MQ=255)
                                                                                                                                    ccACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCa                   >  1:1074112/1‑150 (MQ=255)
                                                                                                                                       ccGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGATGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCg                <  1:1834294/150‑1 (MQ=255)
                                                                                                                                       ccGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCg                <  1:610686/150‑1 (MQ=255)
                                                                                                                                        ctctACCGCCATTTCAATAGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCgt               >  2:976802/4‑150 (MQ=255)
                                                                                                                                         gTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCgtg              <  1:1049994/150‑1 (MQ=255)
                                                                                                                                          ttACCGCTATCGAAATGGCGGTAGAGAAGGCGAACGAGATTGGCGACCCGCTGGTGGTACGTCCGTCATACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCgtgt             <  1:1106664/150‑1 (MQ=255)
                                                                                                                                             ccGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTa          <  1:560097/150‑1 (MQ=255)
                                                                                                                                             ccGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTa          <  1:1620288/150‑1 (MQ=255)
                                                                                                                                             ccGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTa          >  1:1028307/1‑150 (MQ=255)
                                                                                                                                               gCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGCCCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAAc        >  2:1066673/1‑150 (MQ=255)
                                                                                                                                               gCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAAc        >  1:1126578/1‑150 (MQ=255)
                                                                                                                                               gCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGGGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAAc        <  1:291790/150‑1 (MQ=255)
                                                                                                                                                  aTTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGCCCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGAt     >  1:1638932/1‑150 (MQ=255)
                                                                                                                                                    tGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATgc   >  2:1271259/1‑150 (MQ=255)
                                                                                                                                                     gAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATgcg  >  2:1290057/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCG  >  NC_000913/32762‑33060

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