Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,472,437 T→C 7.4% N33D (AAT→GAT)  pyrL ← pyrBI operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,472,4370TC 7.4% 185.9 / 5.2 80N33D (AAT→GAT) pyrLpyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (37/33);  new base C (4/2);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-01

TTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCA  >  NC_000913/4472295‑4472574
                                                                                                                                                |                                                                                                                                         
tttGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGAtt                                                                                                                                           >  2:243307/1‑143 (MQ=255)
     tttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTGTCTCCTGACGCCTGCGCAAAAAAAAGCCCCTCGATTGAggg                                                                                                                                      >  1:639299/1‑143 (MQ=255)
     tttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTCTCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGAGTGAgcg                                                                                                                                      >  1:695636/1‑141 (MQ=255)
       tgaGTTTAGCCGCTGGCGCCCGCACACGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTATGATATAGCGGATTAGGCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGc                                                                                                                                    <  1:587241/141‑1 (MQ=255)
       tCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAGGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCGAAAAAAAGCCCCTCGATTGAGGGGc                                                                                                                                    >  2:651893/1‑143 (MQ=255)
       tCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCGCTCGATTGAGGGGc                                                                                                                                    >  1:145688/1‑143 (MQ=255)
        caggtgcgccgcTGACGCCAGCACCAGAGTAAGGTCATCGCGACTAAGGGCGAATATGGAAATGATATGTTTCTGATATAGCAGATTAGCCATCTTTTATCTCCTGAAAAAAAAAAAAAAAAAAGCCCCTCAATCGAGGGGCt                                                                                                                                   <  2:395018/137‑1 (MQ=255)
         atgttCGCCGCTGTCGCCAGCCCCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCCGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTg                                                                                                                                  <  1:52945/140‑1 (MQ=255)
          gTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGAGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTgg                                                                                                                                 >  1:370038/1‑143 (MQ=255)
           tttCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGGCTAAGGACGTTTATGGAAATGATATGTTTCTGAAAGAGCGGAGTAGCCAGCGATTATCTCCTAACGCGTGGGCAGAAAAAAGCGCCTCGGTTGAgggggtagg                                                                                                                                >  2:209095/1‑138 (MQ=255)
                ccgcTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGACTgg                                                                                                                           >  1:459257/1‑143 (MQ=255)
                ccgcTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGAAATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCGAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATgg                                                                                                                           >  1:517127/1‑143 (MQ=255)
                  gcTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGt                                                                                                                         >  1:541460/1‑143 (MQ=255)
                    tgtggacAGCACTAGGATAAGGTTATCTGTGCTAAGTTAGTATGTGGAAATGATATGTTTCTGATATAGGGGATTAGTCATCTTTTATTTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGa                                                                                                                       <  2:800912/137‑1 (MQ=255)
                             cgccAGACTAAAGTCAGCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCATGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGga                                                                                                              <  2:368361/141‑1 (MQ=255)
                              aCCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTGCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaa                                                                                                             <  2:270095/143‑1 (MQ=255)
                                cAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaaga                                                                                                           >  1:298309/1‑143 (MQ=255)
                                    ttAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaaaa                                                                                                       >  2:679100/1‑143 (MQ=255)
                                    ttAAGGTCATCGCGACTAAGGTCGTTTATGGAAAGGATATGTGTCTGATATAGCGGATTAGCGCATATAGAAGTTCTTACGGCTGGGCCAATAAAATCCCCTCTATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaaaa                                                                                                       <  1:381443/143‑1 (MQ=255)
                                            aTCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCTATCCTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGc                                                                                               <  1:735352/143‑1 (MQ=255)
                                               gcgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGAGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGGGATCAACGAGAGGAACAACAGCAGGCCAg                                                                                            >  2:480657/1‑143 (MQ=255)
                                                ctacgatggTCGCTTTCGGAAATGATATGTTTCGGATACAGCGGATTCGCCGTCTTTTATCTCCTGACGACCGGGCACAACAAGGCCGCTCGATTGTGGGGCTCGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCAAGc                                                                                           <  2:209783/136‑1 (MQ=255)
                                                 gACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCg                                                                                          >  1:800192/1‑143 (MQ=255)
                                                  aCTAAGGTCGGTGATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAAGGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGt                                                                                         <  1:43774/143‑1 (MQ=255)
                                                    tAAGGTCGTTTATGGAAATGATATGTTTCTGATACAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCGAAAAAAAGCCCCTCGATTGCGGGGCTTGGTATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCt                                                                                       >  1:305135/1‑143 (MQ=255)
                                                      atgtagTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCTTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCttt                                                                                     <  2:157647/138‑1 (MQ=255)
                                                        ttagTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGACAAAGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCttttt                                                                                   <  1:782458/140‑1 (MQ=255)
                                                        gTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGACAAAAAAAAGCCCCTCGATTGAGGGTCCGGGAATGGGTGGTCGACGGGAAGAAAACCGGCAGGCCCGCGTCttttt                                                                                   >  1:755384/1‑143 (MQ=255)
                                                         tCGTGTATGCAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtttttt                                                                                  <  2:517127/143‑1 (MQ=255)
                                                          cGTTTATGGACATGATATGTTGCTGATATAGCGGATTAGGCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTc                                                                                 <  2:262757/143‑1 (MQ=255)
                                                             ttATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGa                                                                              <  1:408238/143‑1 (MQ=255)
                                                                  gAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATGTCTTGACGACTGAGCAAAAAAACACCGCATTCAGCGTGGGCTGGGAATGGGTGATCCACGGGAAGACAAACGGCAGGCCAGCGTCTTTTTTCCGACGCgg                                                                         >  1:190695/1‑143 (MQ=255)
                                                                  gAAATGATATGTTTCTGATATAGCGCGTTAGCCATCTTGTATCTCCTGAAGCCTGGGCAAAAAAAAAACCGTGTGATGAGTTGCTGGGAATGGGGGAGCAACGGGAAGAAAAACGGCAGGCGAGCGGCGTTTTTCAGACGCgg                                                                         >  2:584725/1‑143 (MQ=255)
                                                                        atatGTTTCTGATATCGCGGATTAGCCATCTTCCATCGCCGGACGCCTGGGCAAAAAAAAGCCCCTCGATGGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGAc                                                                   <  1:163508/143‑1 (MQ=255)
                                                                          atGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGACTGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGAGGTAAGCCaa                                                                 >  2:158989/1‑143 (MQ=255)
                                                                            gTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTGTTTTCAGCCGCGGTAAGCCaaaa                                                               >  2:5332/1‑143 (MQ=255)
                                                                            gTTGCTGATGTAGGTGATTGGCCATCTTGTATCGCCGGACGCCTGGGCAAAAAATAGCCCCTCGATGGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaaaa                                                               <  2:421696/143‑1 (MQ=255)
                                                                                 tGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCg                                                          <  1:115578/143‑1 (MQ=255)
                                                                                     ataGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTCTCAGACGCGGTAAGACAAAATGTGGAacg                                                      >  2:234050/1‑142 (MQ=255)
                                                                                     ataGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGCCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAaca                                                      <  1:367308/143‑1 (MQ=255)
                                                                                       aGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGTGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacaca                                                    <  1:407604/143‑1 (MQ=255)
                                                                                       aGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacaca                                                    <  2:541460/143‑1 (MQ=255)
                                                                                        gCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGGATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACACAATGTCGAacacac                                                   >  2:682482/1‑143 (MQ=255)
                                                                                        gCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCCAacacac                                                   >  1:266194/1‑143 (MQ=255)
                                                                                         cGGCTTAGCCATCGTTGATCGCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACt                                                  <  1:510683/143‑1 (MQ=255)
                                                                                           gATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACACAATGTCGAAAACACTGa                                                >  2:33968/1‑143 (MQ=255)
                                                                                              tAGCCATCTTTTGTCTCCGGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAAcc                                             <  1:106101/143‑1 (MQ=255)
                                                                                                gCCATCTTGTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGCGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCAt                                           <  2:267477/143‑1 (MQ=255)
                                                                                                 ccATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATa                                          <  1:325305/143‑1 (MQ=255)
                                                                                                 ccATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGAGTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATa                                          <  1:385398/143‑1 (MQ=255)
                                                                                                 ccATCTTTTATCTCCTGAAGCCTGGGCAAAAAAAAGCGCTTCGATTGAGGCACTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAAGCGGTAAGACAAAAATTCGAACACACAGAACCATa                                          >  2:175588/1‑143 (MQ=255)
                                                                                                 ccATCTTCTCTCCCCTGACGCCTGGGCAAAAAAAAGCCCCTCGAGTGGGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATa                                          <  1:41588/143‑1 (MQ=255)
                                                                                                 ccATCCTTTATCTCGTGACGCCTGGGCAAAAAAACGCCCGTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATa                                          <  2:248699/143‑1 (MQ=255)
                                                                                                   cATCTTTTATCTCCTGACGCCTGGGGAAGAAAAGGCGGCTTGAGTGAGGGGGTGGGAATGGGTGATCAAGGGGGAGAGAAACGGCAGGCGAGCGTCTTTTTTCAGGCGCGGTAAGAGAGAATGTCGAGCGCACTGAACCATAc                                         >  1:806431/1‑143 (MQ=255)
                                                                                                    aTCTTTTATCTCCTGACGCCTGGGCAAAAACAATACCTGAAGTTGAGGGGCTGGGGATGGAAGAAGAACAGAAATTAATACGGCATTCCATCGTCAACATTCAGATGCGGAAAGACAAAAGGTCGAACACACTTAACGATACa                                        >  2:290676/1‑143 (MQ=255)
                                                                                                    aTCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAGATGTCGAACCCACTGAACCATACa                                        >  1:502590/1‑143 (MQ=255)
                                                                                                    aTCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAAGCATACa                                        >  2:1454/1‑143 (MQ=255)
                                                                                                     tCTTTTATCTCCTGACGCCTCGGCAAAAAAAGCCCCCCGATTGCGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtc                                      <  1:295333/143‑1 (MQ=255)
                                                                                                       ttttATCTCCTGACGCCTGGGCAAAACAAAGCCCCTCGATTGCGGGGCTGGGAATGGGTGATCAACGGGAATAAAAACGACAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcc                                     >  2:786234/1‑143 (MQ=255)
                                                                                                          tATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcctcc                                  <  1:554873/143‑1 (MQ=255)
                                                                                                           atgcccggACGCGTGGGCAAAAAAAAGCCCCACGATGGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTccc                                 <  1:238327/136‑1 (MQ=255)
                                                                                                              tcCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACGCACTGAACCATACATCCTCCCGGc                              >  1:200420/1‑143 (MQ=255)
                                                                                                                   aCGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAAtt                         >  2:603151/1‑143 (MQ=255)
                                                                                                                   aCGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAAtt                         <  1:225396/143‑1 (MQ=255)
                                                                                                                   aCGCCCGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAAtt                         <  1:243307/143‑1 (MQ=255)
                                                                                                                     gCCTGGGCAAAAAAAATCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCCGCGTCTTTTTTCAGACGCGGTAAGACAAACTTTCGAACACACTGAACCATACATCCTCCCGGCAAATTGt                       <  2:806564/143‑1 (MQ=255)
                                                                                                                     gCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGt                       >  2:160147/1‑143 (MQ=255)
                                                                                                                        tgtgCAGAAAAGAGCCCCTGGGTGGAGGGCCTGGGAATGGGTGGGCAAGGGGAAGAAAAACGGCAGGCCAGCGTCTTTGTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCg                    <  2:462536/140‑1 (MQ=255)
                                                                                                                          ggCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAAGGGGTGATCAAGGGGAAGAAAAACGGCAGTCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGc                  >  1:367313/1‑143 (MQ=255)
                                                                                                                           gCAAAAAAGCGCCCCTCGAGTGAGGGGCAGACCGTGATATACCAACGGTAAGAAAAGCGGCGGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCa                 <  2:169385/143‑1 (MQ=255)
                                                                                                                           gCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGACCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCAGCAAATTATCCGGCa                 >  1:585185/1‑143 (MQ=255)
                                                                                                                             caaaaaaaGCCCCTCGATTGAGGGGCTGGGCATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTCCTCCTTTTGTTTTGTGCAGAAATGTCGAGCACACTGAACCAAACATCCACCCCGCAACTTGTACGGCCtt               >  2:202516/2‑143 (MQ=255)
                                                                                                                              aaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCCttg              >  1:180998/1‑142 (MQ=255)
                                                                                                                               aaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATtat             >  2:417181/1‑143 (MQ=255)
                                                                                                                                aaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAAGAGAATAGACTTAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATtata            <  1:420704/143‑1 (MQ=255)
                                                                                                                                   ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACCCTGAACCATACCTCCTCCCGGCAACTCGTCCGGCATTATACt          >  2:182967/3‑143 (MQ=255)
                                                                                                                                    caGCCCCTCAATCGAGGGGCTGGGAAGGGGTGAGCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTc         >  1:176398/2‑143 (MQ=255)
                                                                                                                                       ccccTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGACCCATACATCCTCCAGGCAGATTGGCCGGCATGATACTCATc      >  2:692562/1‑143 (MQ=255)
                                                                                                                                          cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGCCAAACGCCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACGAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTc   <  2:800192/143‑1 (MQ=255)
                                                                                                                                          cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGAGTCTTTTTTTAGACGCGGTAAGACAAAATGTAGAAGACACTGAACCATAATTCCTCCCGGTAAATTGTCGGGCATTATACTCATCGTc   >  2:286113/1‑143 (MQ=255)
                                                                                                                                           tCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCCGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCa  <  2:535739/143‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                         
TTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCA  >  NC_000913/4472295‑4472574

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: