Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,321,427 A→T 5.9% L507Q (CTG→CAG)  trpE ← anthranilate synthase subunit TrpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,321,4270AT 5.9% 270.0 / 5.3 84L507Q (CTG→CAG) trpEanthranilate synthase subunit TrpE
Reads supporting (aligned to +/- strand):  ref base A (39/40);  new base T (3/2);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

TGGTTGCGGTAAATCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGC  >  NC_000913/1321284‑1321573
                                                                                                                                               |                                                                                                                                                  
tGGTTGCGGTAAATCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAg                                                                                                                                              >  1:1599920/1‑150 (MQ=255)
      cGGTAAATCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCggg                                                                                                                                        <  1:1324012/150‑1 (MQ=255)
         tAAATCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCTGTACTGcgcgctatt                                                                                                                                     <  1:498727/150‑6 (MQ=255)
         tAAATCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCtt                                                                                                                                     <  1:1001964/150‑1 (MQ=255)
            aTCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGt                                                                                                                                  <  1:1080370/150‑1 (MQ=255)
                ccacTTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTAct                                                                                                                              <  2:1532615/150‑2 (MQ=255)
                ccacGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACg                                                                                                                              <  1:793749/150‑1 (MQ=255)
                    gTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGtt                                                                                                                          <  1:1295813/150‑1 (MQ=255)
                         gCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTtcgtc                                                                                                                     >  2:637484/1‑150 (MQ=255)
                          cccATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTtcgtcg                                                                                                                    >  1:6769/1‑150 (MQ=255)
                          cccATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTtcgtcg                                                                                                                    >  1:801202/1‑150 (MQ=255)
                          cccATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTtcgtcg                                                                                                                    <  1:1482575/150‑1 (MQ=255)
                          cccATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCGGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTtcggcg                                                                                                                    >  2:671751/1‑150 (MQ=255)
                                   gcgcAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGAc                                                                                                           <  1:151722/150‑1 (MQ=255)
                                      cAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGc                                                                                                        >  2:1082007/1‑150 (MQ=255)
                                      cAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCAATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGc                                                                                                        >  1:427340/1‑150 (MQ=255)
                                        aCTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCgg                                                                                                      <  2:500480/150‑1 (MQ=255)
                                        aCTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCgg                                                                                                      <  1:27232/150‑1 (MQ=255)
                                        aCTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTGACGGGTTTCGTCGGCTTCCGACTGCgg                                                                                                      >  2:1433282/1‑150 (MQ=255)
                                         cTGATCTGCCAGGTTGTACGCAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGa                                                                                                     <  2:435265/150‑1 (MQ=255)
                                          tGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGaa                                                                                                    >  1:1301466/1‑150 (MQ=255)
                                               tGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGaa                                                                                               >  1:1323899/1‑150 (MQ=255)
                                               tGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGaa                                                                                               >  2:502332/1‑150 (MQ=255)
                                                 ccAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATc                                                                                             <  1:1082007/150‑1 (MQ=255)
                                                     gTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAgg                                                                                         <  1:193751/150‑1 (MQ=255)
                                                      ttGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGa                                                                                        <  2:1461446/150‑1 (MQ=255)
                                                       tGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGAc                                                                                       <  1:960726/150‑1 (MQ=255)
                                                            gTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTacac                                                                                  >  2:966016/1‑150 (MQ=255)
                                                             tAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACAcc                                                                                 >  2:1559592/1‑150 (MQ=255)
                                                             tAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACAcc                                                                                 >  2:1453755/1‑150 (MQ=255)
                                                               aaaGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAg                                                                               <  2:1352901/150‑1 (MQ=255)
                                                                 agagTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCa                                                                             <  1:1623151/150‑1 (MQ=255)
                                                                  gagTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCAc                                                                            >  1:1090661/1‑150 (MQ=255)
                                                                  gagTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGGGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCAc                                                                            <  2:1237525/150‑1 (MQ=255)
                                                                       gATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCt                                                                       <  2:1451896/150‑1 (MQ=255)
                                                                         tattatCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTg                                                                     <  1:160081/150‑1 (MQ=255)
                                                                                  gcagcaGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGc                                                            <  2:801202/150‑1 (MQ=255)
                                                                                  gcagcaGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGc                                                            <  1:637484/150‑1 (MQ=255)
                                                                                  gcagcaGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGc                                                            <  1:1312830/150‑1 (MQ=255)
                                                                                   cagcagAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCg                                                           <  2:467003/150‑1 (MQ=255)
                                                                                       agAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATAc                                                       <  1:1045142/150‑1 (MQ=255)
                                                                                       agAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTGCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATAc                                                       <  1:1039625/150‑1 (MQ=255)
                                                                                        gAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATAcc                                                      <  1:523477/150‑1 (MQ=255)
                                                                                         aaTGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCTCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCATTGACTACACCAGCACCCGCTTGCACGGTTGCGATAGCg                                                     <  2:741686/150‑1 (MQ=255)
                                                                                            gTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGttt                                                  >  1:1094586/1‑150 (MQ=255)
                                                                                              cAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTc                                                >  2:1338256/1‑150 (MQ=255)
                                                                                              cAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTc                                                <  2:160740/150‑1 (MQ=255)
                                                                                                gCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGTACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTcca                                              <  2:622753/150‑1 (MQ=255)
                                                                                                 ccATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTccac                                             <  1:502332/150‑1 (MQ=255)
                                                                                                   aTCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTccacca                                           >  2:692076/1‑150 (MQ=255)
                                                                                                     cAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAgc                                         >  1:253421/1‑150 (MQ=255)
                                                                                                     cAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAgc                                         >  2:444547/1‑150 (MQ=255)
                                                                                                     cAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAgc                                         >  2:160343/1‑150 (MQ=255)
                                                                                                      aGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAgcg                                        >  2:948943/1‑150 (MQ=255)
                                                                                                      aGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAgcg                                        >  2:335080/1‑150 (MQ=255)
                                                                                                       gAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAgcgc                                       >  2:1009434/1‑150 (MQ=255)
                                                                                                       gAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAgcgc                                       >  2:734618/1‑150 (MQ=255)
                                                                                                          aGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAgg                                                                                         >  2:1070682/1‑97 (MQ=255)
                                                                                                          aGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAgg                                                                                         <  1:1070682/97‑1 (MQ=255)
                                                                                                            tctcCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGa                                                          <  1:1189027/126‑1 (MQ=255)
                                                                                                            tctcCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGa                                                          >  2:1189027/1‑126 (MQ=255)
                                                                                                             ctcCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCg                                 >  1:696586/1‑150 (MQ=255)
                                                                                                              tcCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCgg                                <  1:160343/150‑1 (MQ=255)
                                                                                                                    gCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATcaca                          <  2:1323899/150‑1 (MQ=255)
                                                                                                                     cATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACaa                         >  1:1471856/1‑150 (MQ=255)
                                                                                                                      atgatgCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAAt                        >  2:72145/1‑150 (MQ=255)
                                                                                                                          tgCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAg                    <  1:805433/150‑1 (MQ=255)
                                                                                                                          tgCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAg                    >  2:882087/1‑150 (MQ=255)
                                                                                                                             gcgGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGgtg                 <  2:6769/150‑1 (MQ=255)
                                                                                                                                 tGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCga             <  1:72145/150‑1 (MQ=255)
                                                                                                                                  ggCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCgag            >  2:941132/1‑150 (MQ=255)
                                                                                                                                  ggCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCgag            <  2:1329079/150‑1 (MQ=255)
                                                                                                                                    cAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGAt          >  1:322126/1‑150 (MQ=255)
                                                                                                                                    cAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGAt          >  1:563352/1‑150 (MQ=255)
                                                                                                                                     aaTAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATc         >  2:1345266/1‑150 (MQ=255)
                                                                                                                                       tAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGc       >  2:1031911/1‑150 (MQ=255)
                                                                                                                                       gcccgcgcTGTACTGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCgag            <  1:1514555/142‑1 (MQ=255)
                                                                                                                                       gcccgcgcTGTACTGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCgag            >  2:1514555/4‑145 (MQ=255)
                                                                                                                                        cccgcgcTGTACTGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGcc      >  1:721857/3‑150 (MQ=255)
                                                                                                                                        cccgcgcTGTACTGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGcc      >  1:634529/3‑150 (MQ=255)
                                                                                                                                         gcgcgcAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCa     >  2:764541/1‑150 (MQ=255)
                                                                                                                                         gcgcgcAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCa     >  1:224778/1‑150 (MQ=255)
                                                                                                                                           gcgcAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATg   >  1:1441906/1‑150 (MQ=255)
                                                                                                                                            cgcAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATgc  <  2:1599920/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TGGTTGCGGTAAATCACCACGTTATGCCCATTGCTGCGCAACTGATCTGCCAGGTTGTACGTAAAAGAGTCGATATTATCGAGCAGCAGAATGTCAGCCATCAGAAAGTCTCCTGTGCATGATGCGCGGTGGCAATAGCGCGCAGTACAGCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGC  >  NC_000913/1321284‑1321573

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: