Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 928,452 T→A 5.0% Q248H (CAA→CAT cydC ← glutathione/cysteine ABC transporter export permease/ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913928,4520TA 5.0% 342.2 / 4.8 114Q248H (CAA→CATcydCglutathione/cysteine ABC transporter export permease/ATPase
Reads supporting (aligned to +/- strand):  ref base T (46/61);  new base A (3/3);  total (50/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 116 total aligned reads displayed.

GTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAGTTGTTGGCGAT  >  NC_000913/928303‑928600
                                                                                                                                                     |                                                                                                                                                    
gTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACGGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCt                                                                                                                                                      >  2:785504/1‑150 (MQ=255)
 ttACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCtt                                                                                                                                                     >  1:666095/1‑150 (MQ=255)
    cTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCg                                                                                                                                                  <  2:1881120/150‑1 (MQ=255)
         gCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCTCCAATGAGCAGCATTATCGCTTGCGACAAt                                                                                                                                             >  2:1170680/1‑150 (MQ=255)
         gCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAAt                                                                                                                                             >  2:597217/1‑150 (MQ=255)
         gCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAAt                                                                                                                                             <  1:995989/150‑1 (MQ=255)
          ccAGTGATTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGTGTTGAGCATTGCCGCCAACGACGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGAGAAAATg                                                                                                                                            >  1:588772/1‑150 (MQ=255)
           cAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACTGGTCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGc                                                                                                                                           <  1:1234606/150‑1 (MQ=255)
           cAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGc                                                                                                                                           <  1:1365274/150‑1 (MQ=255)
           cAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGc                                                                                                                                           <  2:141111/150‑1 (MQ=255)
           cAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGc                                                                                                                                           >  2:1467898/1‑150 (MQ=255)
             gTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTGTCGCAAGCGATAATGctg                                                                                                                                         <  1:1705039/150‑3 (MQ=255)
              tGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGt                                                                                                                                        <  1:427673/150‑1 (MQ=255)
               gCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTGTCGCAAGCGATAATGctgct                                                                                                                                       <  2:1028080/150‑5 (MQ=255)
               gCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTc                                                                                                                                       >  2:875338/1‑150 (MQ=255)
               gCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTc                                                                                                                                       <  1:1790352/150‑1 (MQ=255)
                cTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGATCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGCGCAGCATTATCGCTTGCGACAATGCGGGCa                                                                                                                                      >  2:167430/1‑150 (MQ=255)
                 ttCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAAGGCGGTCAg                                                                                                                                     >  1:937008/1‑150 (MQ=255)
                  tCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCCGCATTATCGCTTGCGACAATGCGGTCAGt                                                                                                                                    >  1:1392103/1‑150 (MQ=255)
                  tCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGt                                                                                                                                    >  1:455185/1‑150 (MQ=255)
                   cAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCTGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGtt                                                                                                                                   <  1:511296/150‑1 (MQ=255)
                    aaaCGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTc                                                                                                                                  >  2:184542/1‑150 (MQ=255)
                        gcgGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCAAAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCACTTCAGAt                                                                                                                              >  1:222174/1‑150 (MQ=255)
                          ggCTAACGCGCAGAAGACAAACAGGGCAAATAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTg                                                                                                                            <  2:1703043/150‑1 (MQ=255)
                            cTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCCCCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGCCAATGCGGTCAGTTCAGATTGAc                                                                                                                          >  1:897578/1‑150 (MQ=255)
                            cTAACGCGCAGAAGACAAACAGCGCAATAAACGCGCCGGGGAGAGCATTGCCGCCAACGCCGCCAGACGGCATCCACAGCATCAGGATCTCCGCTACCGCGGCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGAc                                                                                                                          <  2:366397/150‑1 (MQ=255)
                                cgcgCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCg                                                                                                                      <  2:1333197/150‑1 (MQ=255)
                                   gcAGAAGACAAACAGGGCAATAAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTg                                                                                                                   <  1:796651/150‑1 (MQ=255)
                                    cAGAAGACTACCAGGGGAATTAAAGCGCCGGGTTGAGCATTGCCGCCTACGCAGCCAGACGACATCCACAGCATCAGGAGCACCGCTAACGCGCCTATGATCAGCTTAATCGCTTGCGACAATGAGGTCAGTTCAGATTGACGGCGTTgc                                                                                                                  <  1:1241510/150‑1 (MQ=255)
                                    cAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTgc                                                                                                                  <  2:258765/150‑1 (MQ=255)
                                    cAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTgc                                                                                                                  <  1:1744581/150‑1 (MQ=255)
                                         gaCAAACAGGTCAATAAAAGCGACGTGTTGATCATTGCCGCCAAAGCCTCCAGACTCCATCCACAGCATCATGATCACCGCTAACGCTCCAATGAGCAGCATTATCGCTTGCTACAATGCGGTCAGATCAGATTTACGGCGTTGCGCTTc                                                                                                             <  1:1106862/150‑1 (MQ=255)
                                            aaaCAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAg                                                                                                          <  2:494568/150‑1 (MQ=255)
                                              aCAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCTTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGcc                                                                                                        <  1:1350569/150‑1 (MQ=255)
                                              aCAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGcc                                                                                                        <  2:1404487/150‑1 (MQ=255)
                                                cgggCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCAt                                                                                                      <  2:543733/149‑1 (MQ=255)
                                                aGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCAt                                                                                                      <  1:1068939/150‑1 (MQ=255)
                                                aGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCAt                                                                                                      <  1:528869/150‑1 (MQ=255)
                                                 gggCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCAtt                                                                                                     <  1:1820563/150‑1 (MQ=255)
                                                   gCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATCAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCGCATTGCCGGCGTTGCGCTTCCGGCCATTGa                                                                                                   >  1:136646/1‑150 (MQ=255)
                                                      aTTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAAtt                                                                                                >  1:1855653/1‑150 (MQ=255)
                                                        tAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTc                                                                                              <  1:1747116/150‑1 (MQ=255)
                                                          aCGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTg                                                                                            <  1:1890052/150‑1 (MQ=255)
                                                            gcgcCGGGTTGTGAAGTGCCGCCCACGCCGCCAGACGCCATCCCCAGCATCAGGCTCACCGCTTACGCGCCACTGAGCAGCATTATCGCTTGGGACTATGCGGGCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTa                                                                                          <  2:282600/150‑1 (MQ=255)
                                                             cgcCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCTCCGCTAACGCGCCAATGAGCAGGATTATCGCTTGCGCCAATGCGGTCTGTTCAGATTGACGGCGTTGCGGGTCGCGCCATTGAGTGTCTGTAt                                                                                         >  1:302712/1‑150 (MQ=255)
                                                             cgcCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTAt                                                                                         <  2:1635308/150‑1 (MQ=255)
                                                             cgcCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTAt                                                                                         <  1:245997/150‑1 (MQ=255)
                                                               ccGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATtc                                                                                       >  2:1595585/1‑150 (MQ=255)
                                                                cGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTGCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATtct                                                                                      <  2:952772/150‑1 (MQ=255)
                                                                        gCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTTACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTtgc                                                                              >  1:980365/1‑150 (MQ=255)
                                                                         cATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTtgcg                                                                             <  2:1791506/150‑1 (MQ=255)
                                                                          aTTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTtgcgt                                                                            <  2:1671395/150‑1 (MQ=255)
                                                                          aTTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCGTGCGACAATGCGGTCAGGTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTtgcgt                                                                            >  1:1106700/1‑150 (MQ=255)
                                                                           ttGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTtgcgtg                                                                           <  2:21201/150‑1 (MQ=255)
                                                                            tGCCGCCAACTCCGCCAGACGCCATCCACAGCATCAGGATCACCGATTACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGTGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTtgcgtgc                                                                          <  1:101622/150‑1 (MQ=255)
                                                                             gccgccAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCCGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCCGCCATTGAATTTCTGTATTCTCTAGTtgcgtgcg                                                                         >  1:551017/1‑150 (MQ=255)
                                                                              ccgccAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGa                                                                        <  2:834565/150‑1 (MQ=255)
                                                                               cgccAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGAt                                                                       <  1:207034/150‑1 (MQ=255)
                                                                               cgccAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGAGAt                                                                       >  2:1782772/1‑150 (MQ=255)
                                                                                  cAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAAc                                                                    <  2:1591835/150‑1 (MQ=255)
                                                                                   aaCGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACg                                                                   <  1:1496648/150‑1 (MQ=255)
                                                                                    aCGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGa                                                                  >  2:1667805/1‑150 (MQ=255)
                                                                                       ccgccAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCCTTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGAtag                                                               >  1:577143/1‑148 (MQ=255)
                                                                                       ccgccAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCg                                                               >  2:1873586/1‑150 (MQ=255)
                                                                                       ccgccAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCg                                                               >  1:773922/1‑150 (MQ=255)
                                                                                           cAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTgg                                                           <  1:653597/150‑1 (MQ=255)
                                                                                                    tCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCaaaaa                                                  >  1:829984/1‑150 (MQ=255)
                                                                                                     ccACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAAt                                                 >  2:205061/1‑150 (MQ=255)
                                                                                                           cATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtcag                                           >  2:92123/1‑150 (MQ=255)
                                                                                                              cAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtcagctc                                        <  1:1196059/150‑1 (MQ=255)
                                                                                                              cAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtcagctc                                        <  1:1706923/150‑1 (MQ=255)
                                                                                                               aGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtcagctca                                       >  2:478320/1‑150 (MQ=255)
                                                                                                               aGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtcagctca                                       >  1:725873/1‑150 (MQ=255)
                                                                                                                ggATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATTGCTGGCACCAAAAATGGtcagctcag                                      >  1:1757827/1‑150 (MQ=255)
                                                                                                                ggATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtcagctcag                                      <  2:774160/150‑1 (MQ=255)
                                                                                                                ggATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtcagctcag                                      >  1:1503357/1‑150 (MQ=255)
                                                                                                                 gATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtcagctcagc                                     >  1:535852/1‑150 (MQ=255)
                                                                                                                         cTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTTCGTGCGATTACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCtt                             >  2:339883/1‑150 (MQ=255)
                                                                                                                         cTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCtt                             <  2:1825235/150‑1 (MQ=255)
                                                                                                                           aaCGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGc                           <  1:581490/150‑1 (MQ=255)
                                                                                                                           aaCGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGc                           <  2:545523/150‑1 (MQ=255)
                                                                                                                           aaCGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATAGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGc                           <  2:90142/150‑1 (MQ=255)
                                                                                                                            aCGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAATAATGGTCAGCTCAGCTTGCCCTTGCa                          <  2:1036053/150‑1 (MQ=255)
                                                                                                                                    actgaccgcaTTGTCGCAAGCGATAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGcc                  >  2:1279736/8‑150 (MQ=255)
                                                                                                                                       gaccgcaTTGTCGCAAGCGATTATGCGGTCAGTTCAGATTGCCGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAACAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTc               >  1:468941/5‑150 (MQ=255)
                                                                                                                                       atcagcaTTATCGCTTGCGACAATGCGGGCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTACCCTTGCAGCCAGGCCGTc               <  2:63855/148‑1 (MQ=255)
                                                                                                                                       agcagcATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGta               >  1:728751/1‑149 (MQ=255)
                                                                                                                                       agcagcATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGACCGTc               >  1:1813757/1‑150 (MQ=255)
                                                                                                                                        gcagcaTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCa              >  2:1798252/1‑150 (MQ=255)
                                                                                                                                         ccgcaTTGTCGCAAGCGATAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAg             >  2:1312128/3‑150 (MQ=255)
                                                                                                                                         cagaaTTATAGCTTTCGACAATGCGGTCAGTTCAGATTGACGGCGGTGCGCTTCCAGACATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAg             <  2:80418/150‑1 (MQ=255)
                                                                                                                                          agcaTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAgt            >  1:1776073/1‑150 (MQ=255)
                                                                                                                                          agcaTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCATTGCAGCCAGGCCGTCAgt            >  1:1430413/1‑150 (MQ=255)
                                                                                                                                            caTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCAAGCAATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAgttg          >  1:1792341/1‑150 (MQ=255)
                                                                                                                                                aTCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGATGCGCTTCCCGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAGTTGTTgg      <  1:758026/150‑1 (MQ=255)
                                                                                                                                                 tCGCATGCGACACTGGGGTCAGGTCCGATTGACGGCGTTGCGCTTCCAGCCATTGACTTTCTGTATTCTCAAGATGCGTGCGATTACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTTCCCTTGCAGCCAGGCCGTCAGTTGTTGGc     <  2:302712/150‑1 (MQ=255)
                                                                                                                                                   gCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAGTTGTTGGCGa   >  1:1892728/1‑150 (MQ=255)
                                                                                                                                                   gCTTGCGACAATGAGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCACTTGTTGGCGa   >  1:1648173/1‑150 (MQ=255)
                                                                                                                                                    cTTGCGACAATGCGGTCAGTTCAGATTGCCGGCGGGGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGAGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAGTTGTTGGCGAt  <  2:937008/150‑1 (MQ=255)
                                                                                                                                                    cTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAGTTGTTGGCGAt  <  2:680458/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAGTTGTTGGCGAT  >  NC_000913/928303‑928600

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: