Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,692,978 A→C 6.4% L429R (CTC→CGC)  uidB ← glucuronide transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,692,9780AC 6.4% 286.5 / 4.8 95L429R (CTC→CGC) uidBglucuronide transporter
Reads supporting (aligned to +/- strand):  ref base A (42/46);  new base C (2/5);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCATGGCCACTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAG  >  NC_000913/1692829‑1693119
                                                                                                                                                      |                                                                                                                                             
gCCATGGCCACTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGa                                                                                                                                               <  2:37496/150‑1 (MQ=255)
    tGGCCACTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCgg                                                                                                                                           >  2:1814/1‑150 (MQ=255)
    tGGCCACTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACAACGATTTCTTTGAATTTTTTATCCGTGAGCgg                                                                                                                                           >  2:264961/1‑150 (MQ=255)
      tacACTATTTTTCTAATAGTTGCACCTTTGATGTTCTTATTATTTTTATTGAATATTAATTTGTTATATCGCTGATTAATTTCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGAt                                                                                                                                         <  1:599622/148‑1 (MQ=255)
       ccACTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATa                                                                                                                                        <  2:928614/150‑1 (MQ=255)
          cTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATaaaa                                                                                                                                     >  1:1714498/1‑150 (MQ=255)
           tATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAAc                                                                                                                                    <  2:595759/150‑1 (MQ=255)
           tATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGCTCACGGATaaaaa                                                                                                                                    <  1:1645170/150‑2 (MQ=255)
           tATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACAATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAAc                                                                                                                                    <  1:1343500/150‑1 (MQ=255)
             tttttCTCATAGTTTCTCCTTTTATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTTATTGCTGCTGCACTTTTTTACGATTTTCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGCTCACGGATaaaaaat                                                                                                                                  <  1:753242/150‑4 (MQ=255)
             tttttCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCa                                                                                                                                  >  1:160321/1‑150 (MQ=255)
             tttttCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCa                                                                                                                                  <  1:503060/150‑1 (MQ=255)
              ttttCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGCTCACGGATaaaaaatt                                                                                                                                 <  2:694781/150‑5 (MQ=255)
                ttCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGat                                                                                                                               <  1:820062/150‑1 (MQ=255)
                 tctcATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGata                                                                                                                              <  2:1074807/150‑1 (MQ=255)
                 tctcATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGCTCACGGATaaaaaattcaa                                                                                                                              <  1:438108/150‑8 (MQ=255)
                 tctcATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGCTCACGGATaaaaaattcaa                                                                                                                              <  1:415331/150‑8 (MQ=255)
                   tcATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCAGTGAGCGGATAAAACCAGataat                                                                                                                            >  1:245126/1‑150 (MQ=255)
                   tcATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCAGTGAGCGGATAAAACCAGataat                                                                                                                            >  1:444102/1‑150 (MQ=255)
                    cATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGataata                                                                                                                           >  2:616867/1‑150 (MQ=255)
                    cATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGATGAACTTTTTTACGATTATCAATTTCAACAACGATTTCTTTGAATTTTTTATCCTTGAGCTGATAAAAACAGataata                                                                                                                           >  1:1492670/1‑150 (MQ=255)
                      tAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAAc                                                                                                                         <  1:924732/150‑1 (MQ=255)
                       aGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACg                                                                                                                        <  2:888129/150‑1 (MQ=255)
                          tGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAAt                                                                                                                     <  1:453856/150‑1 (MQ=255)
                            cACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGc                                                                                                                   <  1:724924/150‑1 (MQ=255)
                             aCCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGcc                                                                                                                  <  1:1048068/150‑1 (MQ=255)
                             aCCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGcc                                                                                                                  <  2:1476646/150‑1 (MQ=255)
                              ccTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGcc                                                                                                                  >  1:360100/1‑150 (MQ=255)
                              ccTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCa                                                                                                                 >  2:1486827/1‑150 (MQ=255)
                              ccTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCAGTGAGCGGATAAAACCAGATAATAACGAATGCCa                                                                                                                 >  1:1444390/1‑150 (MQ=255)
                                 ttGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTa                                                                                                              <  1:178506/150‑1 (MQ=255)
                                  tGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAg                                                                                                             <  2:524375/150‑1 (MQ=255)
                                    aTGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCa                                                                                                           <  1:1486827/150‑1 (MQ=255)
                                      gTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATa                                                                                                         >  1:1350629/1‑150 (MQ=255)
                                       ttCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATaa                                                                                                        <  1:264961/150‑1 (MQ=255)
                                           gATTATTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCg                                                                                                   >  2:504569/1‑150 (MQ=255)
                                                       tGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTaa                                                                                        >  2:342044/1‑150 (MQ=255)
                                                       tGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACAAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTaa                                                                                        >  2:1007762/1‑150 (MQ=255)
                                                        gAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCTCAAGGTACTAAg                                                                                       >  2:880229/1‑150 (MQ=255)
                                                         aaTATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAgg                                                                                      <  1:1664019/150‑1 (MQ=255)
                                                          atatTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGc                                                                                     <  2:1231294/150‑1 (MQ=255)
                                                          atatTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGc                                                                                     <  2:1714498/150‑1 (MQ=255)
                                                           tatTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCa                                                                                    <  2:160186/150‑1 (MQ=255)
                                                              taattaGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAAtt                                                                                 >  2:892122/1‑150 (MQ=255)
                                                                   attGATATCGGAGATTAATTGATTCTTCACTTTTTTACTATTATCAATTTAAACAACGATTTCTTTGAATTTTTTATACGTTAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGAtgt                                                                            <  2:1492670/148‑1 (MQ=255)
                                                                    gTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGAtgtg                                                                           <  2:160321/150‑1 (MQ=255)
                                                                     tGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGc                                                                          <  1:669022/150‑1 (MQ=255)
                                                                      tatatCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCg                                                                         <  1:508938/149‑1 (MQ=255)
                                                                       atatCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCgg                                                                        <  2:1196428/150‑1 (MQ=255)
                                                                           cGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGc                                                                    >  2:66554/1‑150 (MQ=255)
                                                                           cGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGc                                                                    <  1:1814/150‑1 (MQ=255)
                                                                           cGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGc                                                                    >  1:876926/1‑150 (MQ=255)
                                                                           cGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTACGGATGc                                                                    >  1:1521568/1‑150 (MQ=255)
                                                                           cGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGACAGTAGAATAAATCCGCAAGGTAATAAGGCAATTGATGTGCGGATGc                                                                    >  1:62596/1‑150 (MQ=255)
                                                                            gCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGcc                                                                   >  2:1155821/1‑150 (MQ=255)
                                                                                  tAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCataat                                                             >  1:468714/1‑150 (MQ=255)
                                                                                      tgctgctgcACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATACAACCAGATAATAACGAATGCCAGTAGCATAAATCCACAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACt                                                         >  1:97133/1‑150 (MQ=255)
                                                                                       gctgctgcACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACtt                                                        >  1:992857/1‑150 (MQ=255)
                                                                                        ctgctgcACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGTTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTc                                                       >  1:595367/1‑150 (MQ=255)
                                                                                         tgctgcACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGCGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGGGCGGAGGCCCATAATAACTTcc                                                      >  2:1545970/1‑150 (MQ=255)
                                                                                         tgctgcACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTcc                                                      >  1:666670/1‑150 (MQ=255)
                                                                                          gctgcACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCg                                                     <  1:832509/150‑1 (MQ=255)
                                                                                          gctgcACTTTTTTACGATTAACAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAAGCAGATAATAACTAATGCCAGTATCATAAATCCTAAAGGTACTAAGCCAATTGATGAGAGTATGCCCATAATAACTTCCg                                                     >  2:1242801/1‑150 (MQ=255)
                                                                                           ctgcACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCgg                                                    <  1:1187443/150‑1 (MQ=255)
                                                                                             tcACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTTAATTTTTTATCCTTGAGCGGATAAAACCAGATAATAACGAATGCCATAAGCATAAATCCGCAATGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCg                                                  <  1:1235387/149‑1 (MQ=255)
                                                                                             gcACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCg                                                  >  2:1524447/1‑150 (MQ=255)
                                                                                             gcACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCg                                                  <  2:1195862/150‑1 (MQ=255)
                                                                                                cTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGttt                                               >  2:1579721/1‑150 (MQ=255)
                                                                                                cTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAAACAGAAAATAACGAATGCCACTAGCATAAATCCGCAAGGTACTAAGGCAATTCATGTGCGGATGCCCATAATAACTTCAggtcttt                                               >  2:414873/1‑145 (MQ=255)
                                                                                                 ttttttACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTg                                              <  1:29259/150‑1 (MQ=255)
                                                                                                  tttttACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGc                                             >  2:809038/1‑150 (MQ=255)
                                                                                                     ttACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCAcc                                          >  2:601389/1‑150 (MQ=255)
                                                                                                        cGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGa                                       >  2:1498808/1‑150 (MQ=255)
                                                                                                        cGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGa                                       >  1:978869/1‑150 (MQ=255)
                                                                                                             aTCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGc                                  >  1:1365849/1‑150 (MQ=255)
                                                                                                             aTCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGc                                  <  2:1546195/150‑1 (MQ=255)
                                                                                                              tCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCg                                 <  1:1071547/150‑1 (MQ=255)
                                                                                                                      aaCCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATcc                         <  2:1350629/150‑1 (MQ=255)
                                                                                                                              tttCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAAcc                 <  1:342044/150‑1 (MQ=255)
                                                                                                                                 cTTTGAATTTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAAcccc               <  2:360100/150‑1 (MQ=255)
                                                                                                                                  tttGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCaa             >  2:1349576/1‑150 (MQ=255)
                                                                                                                                     gAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTACGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGGTATATCCGCTTAACCCCAAAAt          >  1:1665133/1‑150 (MQ=255)
                                                                                                                                     gAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCATAAt          >  1:999671/1‑150 (MQ=255)
                                                                                                                                     gAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAAt          <  1:834999/150‑1 (MQ=255)
                                                                                                                                      aaTTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATa         <  1:1308555/150‑1 (MQ=255)
                                                                                                                                         ttttttATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATaaa       <  1:681564/150‑1 (MQ=255)
                                                                                                                                         ttttttATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATaaa       <  1:187391/150‑1 (MQ=255)
                                                                                                                                         ggttttATCCGCTCACGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATaaa       >  1:1644108/3‑150 (MQ=255)
                                                                                                                                          tttttATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAg      >  2:756624/1‑150 (MQ=255)
                                                                                                                                          tttttATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAg      <  1:1657948/150‑1 (MQ=255)
                                                                                                                                          tttttATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAg      <  1:1155821/150‑1 (MQ=255)
                                                                                                                                            tttATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGc    <  2:245126/150‑1 (MQ=255)
                                                                                                                                            tttATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGc    <  2:444102/150‑1 (MQ=255)
                                                                                                                                             ttATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCa   <  2:978869/150‑1 (MQ=255)
                                                                                                                                             ttATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCACATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCa   >  2:1206486/1‑150 (MQ=255)
                                                                                                                                              tATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAg  >  2:1099273/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                             
GCCATGGCCACTATTTTTCTCATAGTTGCACCTTTGATGTTCTGATTATTTTTATTGAATATTAATTAGTGATATCGCTGATTAATTGCTGCTGCACTTTTTTACGATTATCAATTTCAACCACGATTTCTTTGAATTTTTTATCCGTGAGCGGATAAAACCAGATAATAACGAATGCCAGTAGCATAAATCCGCAAGGTACTAAGGCAATTGATGTGCGGATGCCCATAATAACTTCCGGCGTTTGCACCTGATTGGCGATATATCCGCTTAACCCCAAAATAAAGGCAG  >  NC_000913/1692829‑1693119

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: