Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 847,215 A→T 8.4% intergenic (‑96/+43) glnP ← / ← glnH glutamine transporter subunit/glutamine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913847,2150AT 8.4% 214.3 / 9.6 84intergenic (‑96/+43)glnP/glnHglutamine transporter subunit/glutamine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (40/37);  new base T (2/5);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-02

ACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCG  >  NC_000913/847066‑847360
                                                                                                                                                     |                                                                                                                                                 
aCCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATa                                                                                                                                                   <  1:1182318/150‑1 (MQ=255)
         cAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGGAGAgg                                                                                                                                          <  2:1281576/150‑1 (MQ=255)
         cAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGGAGAgg                                                                                                                                          <  1:1395748/150‑1 (MQ=255)
          aGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAggg                                                                                                                                         >  1:187185/1‑150 (MQ=255)
           taaGCTGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATTATCAAATTCAGGGTAATTAATGGCGCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGa                                                                                                                                        <  1:1587500/149‑1 (MQ=255)
           gAAGCTGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATTGGGAGAggat                                                                                                                                        <  2:684105/150‑3 (MQ=255)
           gAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGGAGAggat                                                                                                                                        <  1:747964/150‑3 (MQ=255)
           gAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGa                                                                                                                                        >  2:571399/1‑150 (MQ=255)
           gAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGa                                                                                                                                        <  1:1117574/150‑1 (MQ=255)
             aGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAAc                                                                                                                                      <  1:984676/150‑1 (MQ=255)
               cGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATAAAGGGCCCCGGAGCGGGGTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTg                                                                                                                                    <  2:1044985/150‑1 (MQ=255)
                ggAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGGAGAggattaggg                                                                                                                                   <  2:1405252/150‑8 (MQ=255)
                   aTGCCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAg                                                                                                                                <  2:69700/150‑1 (MQ=255)
                        aGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGt                                                                                                                           >  1:687226/1‑150 (MQ=255)
                         ggCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCCAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGACCTGCCAGGGCggt                                                                                                                          >  1:1396249/1‑148 (MQ=255)
                          gCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGCGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGttt                                                                                                                         >  1:150647/1‑150 (MQ=255)
                           ccAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGtttt                                                                                                                        <  1:667236/150‑1 (MQ=255)
                             aGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTAc                                                                                                                      <  1:500935/150‑1 (MQ=255)
                              gATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACa                                                                                                                     >  2:224121/1‑150 (MQ=255)
                               ttGGCACTCCAGTCAAACTGCATAAGTTGTACCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAg                                                                                                                    <  1:628544/149‑1 (MQ=255)
                                tGGCTCTCCAGTCAAACTGCATATGTTGTTCGTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAgg                                                                                                                   <  2:794469/150‑1 (MQ=255)
                                 ggCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGgt                                                                                                                  <  2:342994/150‑1 (MQ=255)
                                   cACTCCAGTCAAACTGCATATGTTGTTCATGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGCGAGAGGGAACTGCCAGTGCGATTTACAGgtgt                                                                                                                >  1:366052/1‑150 (MQ=255)
                                    aCTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTa                                                                                                               <  1:100310/150‑1 (MQ=255)
                                     cTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAg                                                                                                              >  1:1401631/1‑150 (MQ=255)
                                        cAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGt                                                                                                           <  2:1243610/150‑1 (MQ=255)
                                         aGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGtt                                                                                                          >  2:1093952/1‑150 (MQ=255)
                                         aGTCAAACTGCATATGTTGTTCCTGTGACCGTGGTGTAAATAATCAAAGTCAGGGTAAGTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCAAATCCTCTCACCATAGGGGAGAGGGAACTGCCAGGGCGGTGTACAGGTGTAGCGtt                                                                                                          >  2:142140/1‑150 (MQ=255)
                                          gTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGtta                                                                                                         <  1:888016/150‑1 (MQ=255)
                                             aaaCTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGttatta                                                                                                      <  2:526572/150‑1 (MQ=255)
                                             aaaCTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGttatta                                                                                                      >  2:908612/1‑150 (MQ=255)
                                                cTGCATATGTTGTTCCTGTTACCGGGGTGTAAATAATCCAATTCATGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTACTCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTAttt                                                                                                   <  1:293034/150‑1 (MQ=255)
                                                 tGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTc                                                                                                  <  2:30481/150‑1 (MQ=255)
                                                    atatGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGt                                                                                               <  1:1003666/150‑1 (MQ=255)
                                                      ttgttgtTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTc                                                                                             <  2:796506/149‑1 (MQ=255)
                                                      atGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTc                                                                                             >  2:1185423/1‑150 (MQ=255)
                                                       tgttgtTCCTGATACCGTGGTGTAAATAATCAAATTCCGGGTAATTTATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCACTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACATGTGTAGCGTTATTATTTCGGTTCa                                                                                            <  1:1119031/150‑1 (MQ=255)
                                                           gttCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTAc                                                                                        >  1:467256/1‑150 (MQ=255)
                                                            ttCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTAcc                                                                                       <  2:1246311/150‑1 (MQ=255)
                                                             tCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCg                                                                                      >  1:578914/1‑150 (MQ=255)
                                                             tCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCg                                                                                      <  1:353672/150‑1 (MQ=255)
                                                               cTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGaa                                                                                    >  1:822569/1‑150 (MQ=255)
                                                                tGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAAc                                                                                   >  2:1451118/1‑150 (MQ=255)
                                                                  ttACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCa                                                                                 <  1:224121/150‑1 (MQ=255)
                                                                      cGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAtttt                                                                             >  1:32468/1‑150 (MQ=255)
                                                                           tgtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTa                                                                        <  2:657000/150‑1 (MQ=255)
                                                                            gtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAg                                                                       >  1:943917/1‑150 (MQ=255)
                                                                                  aatCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCg                                                                 <  2:254178/150‑1 (MQ=255)
                                                                                   atCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                >  1:917563/1‑150 (MQ=255)
                                                                                     cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTg                                                              <  2:822569/150‑1 (MQ=255)
                                                                                       aaTTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTa                                                            <  2:187185/150‑1 (MQ=255)
                                                                                          tCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGt                                                         <  2:1427434/150‑1 (MQ=255)
                                                                                           cAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGtt                                                        >  1:821923/1‑150 (MQ=255)
                                                                                               gTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGt                                                    <  1:562273/150‑1 (MQ=255)
                                                                                                        ggCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCa                                           >  2:796166/1‑150 (MQ=255)
                                                                                                         gCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTTCCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAg                                          >  2:1312538/1‑150 (MQ=255)
                                                                                                          ccccGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAgg                                         >  2:19590/1‑150 (MQ=255)
                                                                                                           cccGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAggg                                        >  1:1558042/1‑150 (MQ=255)
                                                                                                               gAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAGCGAGTTTTTGGAGATGTCGTTGTAAGTTCCGTTCTCGCGCAGGGtttt                                    >  1:286272/1‑150 (MQ=255)
                                                                                                                   ggggTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAAc                                <  2:145218/150‑1 (MQ=255)
                                                                                                                    ggTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAAcg                               <  2:943917/150‑1 (MQ=255)
                                                                                                                    ggTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAGTTTTTGTATATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAAcg                               >  1:865413/1‑150 (MQ=255)
                                                                                                                      tttGCGCCCCTCACCCTAATCCGCTCCCCCTAGGGGAGAGGGAACTGCCACTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAAcgcg                             <  1:69091/150‑1 (MQ=255)
                                                                                                                            ccccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTg                       <  2:185943/150‑1 (MQ=255)
                                                                                                                            ccccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTg                       >  1:430045/1‑150 (MQ=255)
                                                                                                                              ccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTATCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGAc                     <  2:768517/150‑1 (MQ=255)
                                                                                                                               cTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACt                    <  2:824921/150‑1 (MQ=255)
                                                                                                                                 cACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACttt                  >  2:313109/1‑150 (MQ=255)
                                                                                                                                  aCCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTg                 >  1:1062148/1‑150 (MQ=255)
                                                                                                                                  aCCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTg                 >  2:54572/1‑150 (MQ=255)
                                                                                                                                  aCCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGGGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTg                 >  1:1067648/1‑150 (MQ=255)
                                                                                                                                   cccTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGt                >  1:901959/1‑150 (MQ=255)
                                                                                                                                    ccTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTc               >  2:782601/1‑150 (MQ=255)
                                                                                                                                     cTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCa              >  2:486663/1‑150 (MQ=255)
                                                                                                                                      gttcccTCTCCCCTATGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCAc             >  2:1233893/5‑150 (MQ=255)
                                                                                                                                       aaTCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACg            >  2:563839/1‑150 (MQ=255)
                                                                                                                                          ccTCTCCCCTATGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAg         >  2:1051925/1‑150 (MQ=255)
                                                                                                                                           ctctCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGACCCATTTTTTGTATATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCCCGCAGc        >  2:199513/1‑150 (MQ=255)
                                                                                                                                           ctctCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGc        >  1:1076658/1‑150 (MQ=255)
                                                                                                                                            tctcCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCt       >  2:1726390/1‑150 (MQ=255)
                                                                                                                                             ctcCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtc      <  2:629000/150‑1 (MQ=255)
                                                                                                                                               ccccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgt    <  2:150647/150‑1 (MQ=255)
                                                                                                                                                cccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGAGTTATTATTTCGGTTCAGTACCGAACCAGGGTGTGTAGAATTCGTTGTAAGTAACGTTCTCGAGCAGGGTTTTCAACGAGCCGATGACATTGTCACGCCGCtcgtc   >  2:1066330/1‑150 (MQ=255)
                                                                                                                                                 ccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtcg  >  2:211216/1‑150 (MQ=255)
                                                                                                                                                 ccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtcg  <  1:19590/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCG  >  NC_000913/847066‑847360

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