Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 770,810 T→C 7.4% intergenic (+199/‑648) mngB → / → cydA alpha‑mannosidase/cytochrome d terminal oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913770,8100TC 7.4% 448.3 / 16.7 143intergenic (+199/‑648)mngB/cydAalpha‑mannosidase/cytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base T (65/67);  new base C (5/6);  total (70/73)
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.28e-01

Only 100 of 151 total aligned reads displayed.

GCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCAAA  >  NC_000913/770662‑770955
                                                                                                                                                    |                                                                                                                                                 
gcgcTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTatat                                                                                                                                                   <  1:1633655/149‑1 (MQ=255)
 cgcTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATAtg                                                                                                                                                  <  1:480373/149‑1 (MQ=255)
  gcTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATAtgt                                                                                                                                                 <  1:598841/149‑1 (MQ=255)
    tttGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTgg                                                                                                                                               <  4:847809/149‑1 (MQ=255)
    tttGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTgg                                                                                                                                               >  1:467521/1‑149 (MQ=255)
       gTTCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCtt                                                                                                                                            >  3:435369/1‑149 (MQ=255)
        ttCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCttt                                                                                                                                           >  2:467827/1‑149 (MQ=255)
        ttCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCttt                                                                                                                                           >  3:556718/1‑149 (MQ=255)
        ttCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCttt                                                                                                                                           >  3:1667920/1‑149 (MQ=255)
         tCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTa                                                                                                                                          >  2:1410566/1‑149 (MQ=255)
         tCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACATATACGTGGCttga                                                                                                                                          <  4:1253440/149‑3 (MQ=255)
           aTGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACATATACGTGGCttgaga                                                                                                                                        <  1:490081/149‑5 (MQ=255)
              ccGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAttttt                                                                                                                                     <  1:563722/149‑1 (MQ=255)
                ggATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTaa                                                                                                                                   <  2:467521/149‑1 (MQ=255)
                  aTACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACa                                                                                                                                 >  2:1514044/1‑149 (MQ=255)
                      ggCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAaat                                                                                                                             >  3:94524/1‑149 (MQ=255)
                       gCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAaata                                                                                                                            >  2:1311226/1‑149 (MQ=255)
                          tGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAataata                                                                                                                         >  4:1634311/1‑149 (MQ=255)
                           gAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAAtaataa                                                                                                                        >  3:530835/1‑149 (MQ=255)
                            aaCGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAAc                                                                                                                       <  3:927863/149‑1 (MQ=255)
                             aCGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAAcc                                                                                                                      >  3:1639741/1‑149 (MQ=255)
                               gCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCt                                                                                                                    >  3:19915/1‑149 (MQ=255)
                               gCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCt                                                                                                                    <  1:116776/149‑1 (MQ=255)
                                 tttATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTgg                                                                                                                  >  4:208247/1‑149 (MQ=255)
                                 tttATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTgg                                                                                                                  <  1:390323/149‑1 (MQ=255)
                                 tttATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTgg                                                                                                                  <  4:421196/149‑1 (MQ=255)
                                  ttATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTggg                                                                                                                 >  2:1248985/1‑149 (MQ=255)
                                    aTTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAttt                                                                                                                                       <  4:1655470/125‑1 (MQ=255)
                                    aTTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAttt                                                                                                                                       >  3:1655470/1‑125 (MQ=255)
                                     ttCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGa                                                                                                              >  2:428169/1‑149 (MQ=255)
                                     ttCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACATATACGTGGCttga                                                                                                                                          <  2:889113/121‑3 (MQ=255)
                                     ttCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACATATACGTGGCttga                                                                                                                                          >  1:889113/1‑119 (MQ=255)
                                       cGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGt                                                                                                            <  4:676372/149‑1 (MQ=255)
                                           cTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAAt                                                                                                        <  2:1530357/149‑1 (MQ=255)
                                           cTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAAt                                                                                                        <  2:117090/149‑1 (MQ=255)
                                             aCAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCATTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATta                                                                                                      <  4:716677/149‑1 (MQ=255)
                                              cAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGa                                                                                                              <  1:436199/140‑1 (MQ=255)
                                               aaaaGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATtata                                                                                                    >  2:1074047/1‑149 (MQ=255)
                                               aaaaGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATtata                                                                                                    >  3:26726/1‑149 (MQ=255)
                                               aaaaGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACATATACGTGGCttga                                                                                                                                          >  2:743838/1‑109 (MQ=255)
                                               aaaaGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACATATACGTGGCttga                                                                                                                                          <  1:743838/111‑3 (MQ=255)
                                                aaaGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATtataa                                                                                                   <  1:747304/149‑1 (MQ=255)
                                                           aTTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATaa                                                                                        <  4:435369/149‑1 (MQ=255)
                                                           aTTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATaa                                                                                        <  1:1452349/149‑1 (MQ=255)
                                                            ttCAATATATTGCAGAGATGGTGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAg                                                                                       <  2:797256/149‑1 (MQ=255)
                                                            ttCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAg                                                                                       >  3:1526220/1‑149 (MQ=255)
                                                                atatatTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTaa                                                                                   <  3:339117/149‑1 (MQ=255)
                                                                   tatTGCAGAGATGATGTAGACACTGGAAAGCGTAGCGCATCAGGTAATTTTGGGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATATTTATAAGTTAACTa                                                                                >  4:372317/1‑149 (MQ=255)
                                                                      tGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAAt                                                                             <  4:1392737/149‑1 (MQ=255)
                                                                      tGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTATAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAAt                                                                             <  1:467827/149‑1 (MQ=255)
                                                                       gCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATg                                                                            >  4:1635457/1‑149 (MQ=255)
                                                                         agagATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGtt                                                                          >  4:1465656/1‑149 (MQ=255)
                                                                            gatgatGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAat                                                                       >  1:1431014/1‑149 (MQ=255)
                                                                               gatgTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGt                                                                           <  2:1167017/142‑1 (MQ=255)
                                                                               gatgTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGt                                                                           >  1:1167017/1‑142 (MQ=255)
                                                                                atgTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAAc                                                                                  >  3:1366821/1‑134 (MQ=255)
                                                                                atgTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAAc                                                                                  <  4:1366821/134‑1 (MQ=255)
                                                                                atgTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTg                                                                   <  1:1057193/149‑1 (MQ=255)
                                                                                atgTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTg                                                                   <  4:1526220/149‑1 (MQ=255)
                                                                                    aGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTgg                                                                                                                  <  4:1654846/98‑1 (MQ=255)
                                                                                       cacTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGt                                                            >  3:1469963/1‑149 (MQ=255)
                                                                                       cacTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGt                                                            <  4:684329/149‑1 (MQ=255)
                                                                                        acTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTaa                                                                                                         >  2:778917/1‑103 (MQ=255)
                                                                                           ggCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGAt                                                        <  4:1311266/149‑1 (MQ=255)
                                                                                              aaGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATta                                                                                           >  1:1395756/1‑111 (MQ=255)
                                                                                              aaGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATta                                                                                           <  2:1395756/111‑1 (MQ=255)
                                                                                                    aGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCtt                                               <  3:739959/149‑1 (MQ=255)
                                                                                                     gcgcATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCttt                                              <  1:151847/149‑1 (MQ=255)
                                                                                                         aTCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTttatta                                          <  4:1469963/149‑1 (MQ=255)
                                                                                                         aTCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTttatta                                          <  2:86855/149‑1 (MQ=255)
                                                                                                              gTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAAt                                     <  2:1514408/149‑1 (MQ=255)
                                                                                                               tAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATc                                    <  3:208247/149‑1 (MQ=255)
                                                                                                                        gTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAActc                           >  1:1119634/1‑149 (MQ=255)
                                                                                                                        gTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAActc                           >  3:709691/1‑149 (MQ=255)
                                                                                                                         tttATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAActct                          >  2:715525/1‑149 (MQ=255)
                                                                                                                          ttATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTTTAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTg                         <  2:240991/149‑1 (MQ=255)
                                                                                                                           tATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGc                        <  4:197329/149‑1 (MQ=255)
                                                                                                                             tCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCggg                                                              <  1:1239825/109‑1 (MQ=255)
                                                                                                                             tCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCATGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCggg                                                              >  2:1239825/1‑109 (MQ=255)
                                                                                                                              cTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGcgtc                     <  2:1124951/149‑1 (MQ=255)
                                                                                                                               ttCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAAc                                                                                  <  1:1341945/87‑1 (MQ=255)
                                                                                                                               ttCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAAt                                                                             <  3:1405833/92‑1 (MQ=255)
                                                                                                                               ttCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGcgtcg                    >  4:359479/1‑149 (MQ=255)
                                                                                                                               ttCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGcgtcg                    >  1:1527407/1‑149 (MQ=255)
                                                                                                                                tCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGcgtcgt                   >  4:612843/1‑149 (MQ=255)
                                                                                                                                  aCTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATtataa                                                                                                   >  1:1627470/1‑67 (MQ=255)
                                                                                                                                  aCTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTTTGCATTTTTTAGTAATCCTGAAACTCTCCGTCGTAt                 >  3:1052729/1‑149 (MQ=255)
                                                                                                                                    tctcAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAAc                                                                                  <  4:782066/82‑1 (MQ=255)
                                                                                                                                    aataaaGCCACATATACGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTa               >  3:1645261/5‑149 (MQ=255)
                                                                                                                                      tcAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCt                                  <  3:1084838/128‑1 (MQ=255)
                                                                                                                                      tcAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCt                                  >  4:1084838/1‑128 (MQ=255)
                                                                                                                                      tcAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGc             <  1:1248985/149‑1 (MQ=255)
                                                                                                                                      tcAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGATGAGTTAATTATTATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGc             >  4:1239980/1‑149 (MQ=255)
                                                                                                                                        aaaccACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCa           >  4:1528742/4‑149 (MQ=255)
                                                                                                                                          gCCACATATACGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATtataatataa                                                                                              >  3:1538410/1‑64 (MQ=39)
                                                                                                                                           ccACCTATATGTGGCTTTATTTTTAAAAAAATAATAACCTTGGGTGACGTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTAGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTg        <  4:1052729/149‑1 (MQ=255)
                                                                                                                                            cACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGa       <  3:372317/149‑1 (MQ=255)
                                                                                                                                            cACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGa       <  4:94524/149‑1 (MQ=255)
                                                                                                                                                 atatGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCaaa  >  3:881839/1‑149 (MQ=255)
                                                                                                                                                  tatGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGt                                                            <  2:158004/90‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCAAA  >  NC_000913/770662‑770955

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: