Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,033,822 T→A 7.9% L30Q (CTG→CAG)  dfp → bifunctional 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,033,8220TA 7.9% 215.0 / 5.6 72L30Q (CTG→CAG) dfpbifunctional 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase
Reads supporting (aligned to +/- strand):  ref base T (30/36);  new base A (3/3);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCGCAGGCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGG  >  NC_002947/6033680‑6033956
                                                                                                                                               |                                                                                                                                       
cgccgcAGGCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCAc                                                                                                                                  >  1:169250/1‑150 (MQ=255)
      aGGCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CCCCTGGAGCACGgcgcg                                                                                                                            >  2:175092/1‑150 (MQ=255)
        gCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGTGCTCCAGGAG‑ACGgcgaat                                                                                                                           <  1:671911/149‑4 (MQ=255)
        gCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGTGCTCCAGGAG‑ACGgcgaat                                                                                                                           >  2:671911/1‑146 (MQ=255)
        gCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCa                                                                                                                          <  2:169250/150‑1 (MQ=255)
          aTTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCTAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATGCGCCGTGCTCCAGGAG‑ACGgcgaac                                                                                                                           <  2:1126264/147‑4 (MQ=255)
          aTTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATGCGCCGTGCTCCAGGAG‑ACGgcgaat                                                                                                                           >  1:1126264/1‑144 (MQ=255)
           tttCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGt                                                                                                                       >  1:642034/1‑150 (MQ=255)
           tttCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGt                                                                                                                       >  2:1080830/1‑150 (MQ=255)
               aaaCGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTgcgc                                                                                                                   >  2:660969/1‑150 (MQ=255)
               aaaCGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTAGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTgcgc                                                                                                                   >  1:627936/1‑150 (MQ=255)
                 aCGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGcgt                                                                                                                 >  1:798325/1‑150 (MQ=255)
                     cgtcgCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCgtcgtc                                                                                                             <  1:660969/150‑1 (MQ=255)
                     cgtcgCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCgtcgtc                                                                                                             <  2:798325/150‑1 (MQ=255)
                             gCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCg                                                                                                     >  2:1610784/1‑150 (MQ=255)
                               tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCa                                                                                                            >  2:1379053/1‑141 (MQ=255)
                               tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCa                                                                                                            <  1:1379053/141‑1 (MQ=255)
                               tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCggg                                                                                                   <  1:1610784/150‑1 (MQ=255)
                               tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCggg                                                                                                   <  1:175092/150‑1 (MQ=255)
                               tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCggg                                                                                                   <  2:642034/150‑1 (MQ=255)
                                 tcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGggt                                                                                                 <  1:1080830/150‑1 (MQ=255)
                                  cacggtaaaaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGggtg                                                                                                <  2:1579029/144‑1 (MQ=255)
                                     tGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGgtggtg                                                                                             >  1:500775/1‑150 (MQ=255)
                                     tGGAAAAAAGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGgtggtg                                                                                             <  1:1040201/150‑1 (MQ=255)
                                      ggAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGc                                                                                            <  1:1238729/150‑1 (MQ=255)
                                       gAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑a                                                                                          >  2:603615/1‑150 (MQ=255)
                                         aaaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCggg                                                                                                   <  2:8499/140‑1 (MQ=255)
                                         aaaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCggg                                                                                                   >  1:8499/1‑140 (MQ=255)
                                         aaaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGa                                                                                        >  2:1096291/1‑150 (MQ=255)
                                          aaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTTGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGaa                                                                                       <  2:1186726/150‑1 (MQ=255)
                                          aaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGggtg                                                                                                >  1:665586/1‑142 (MQ=255)
                                          aaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGggtg                                                                                                <  2:665586/142‑1 (MQ=255)
                                            aGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAAtt                                                                                     <  1:1068958/150‑1 (MQ=255)
                                             ggCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATtc                                                                                    <  2:1209922/150‑1 (MQ=255)
                                              gCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATtca                                                                                   >  1:805075/1‑150 (MQ=255)
                                                                 gTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCa                                                                <  2:884976/150‑1 (MQ=255)
                                                                 gTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCa                                                                <  2:1131343/150‑1 (MQ=255)
                                                                    tCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAg                                                               >  2:196885/1‑148 (MQ=255)
                                                                    tCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAg                                                               <  1:196885/148‑1 (MQ=255)
                                                                      gCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGACGCGCACCTGCGcgccgt                                                                                                              <  2:1563914/100‑1 (MQ=255)
                                                                      gCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGACGCGCACCTGCGcgccgt                                                                                                              >  1:1563914/1‑100 (MQ=255)
                                                                      gCAAGCGCATCATTCTCGCCGTGGGAGGCCGAATTGCCGCCTCCGAAAGCGCAGAGCTGATTCGCCGT‑CCCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGGATTCATCACCCCGCTGACCCTGCGGgcgc                                                           >  2:525706/1‑150 (MQ=255)
                                                                            gcATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAg                                                               >  1:1056818/1‑140 (MQ=255)
                                                                            gcATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAg                                                               <  2:1056818/140‑1 (MQ=255)
                                                                            gcATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCgg                                                     <  2:1024416/150‑1 (MQ=255)
                                                                            gcATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCgg                                                     <  1:1096291/150‑1 (MQ=255)
                                                                             cATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGAGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCggg                                                    <  2:805075/150‑1 (MQ=255)
                                                                              aTCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGc                                                   <  2:132726/150‑1 (MQ=255)
                                                                                 gTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCAc                                                >  1:340225/1‑150 (MQ=255)
                                                                                  ttCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATTACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCAcc                                               <  1:611264/150‑1 (MQ=255)
                                                                                     tcGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCg                                            >  2:1608082/1‑150 (MQ=255)
                                                                                         cGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTg                                                                                                                                        <  1:557234/55‑1 (MQ=255)
                                                                                         cGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTg                                                                                                                                        >  2:557234/1‑55 (MQ=255)
                                                                                                ggcggcATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGAc                                 <  2:340225/150‑1 (MQ=255)
                                                                                                   ggcATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGTGCCACCCCGTACACATGGACCTg                              >  1:1376674/1‑150 (MQ=255)
                                                                                                    gcATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGc                             <  1:408737/150‑1 (MQ=255)
                                                                                                         gccgccTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGATGTCGGGCCACCCCGTACACATGGACCTGCTCGAc                        >  2:877428/1‑150 (MQ=255)
                                                                                                                        gcCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCt         >  1:1169672/1‑150 (MQ=255)
                                                                                                                         ccGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTg        >  2:602510/1‑150 (MQ=255)
                                                                                                                          cGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCCGc       <  1:877428/150‑1 (MQ=255)
                                                                                                                           gAGCTGATTCTCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTa                                          >  2:1537459/1‑114 (MQ=255)
                                                                                                                           gAGCTGATTCTCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTa                                          <  1:1537459/114‑1 (MQ=255)
                                                                                                                              cTGATTCGCCGT‑CCCCTGGAGCACGGCGCGCAGGTGCGCGTCGTAATGACCC‑GGACGACGCGAGAATGCAGCTCCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATg   <  1:603615/150‑1 (MQ=255)
                                                                                                                               tGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATgg  <  2:500775/150‑1 (MQ=255)
                                                                                                                                 aTTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑a                                                                                          <  2:1554105/60‑1 (MQ=255)
                                                                                                                                 aTTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑a                                                                                          >  1:1554105/1‑60 (MQ=255)
                                                                                                                                   gcGCCGTGCTCCAGGAG‑ACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATtc                                                                                    <  2:968879/63‑1 (MQ=25)
                                                                                                                                   gcGCCGTGCTCCAGGAG‑ACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATtc                                                                                    >  1:968879/2‑64 (MQ=25)
                                                                                                                                              cTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCa     >  2:693868/1‑133 (MQ=255)
                                                                                                                                              cTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCa     >  2:693807/1‑133 (MQ=255)
                                                                                                                                              cTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCa     <  1:693868/133‑1 (MQ=255)
                                                                                                                                              cTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCa     <  1:693807/133‑1 (MQ=255)
                                                                                                                                               tGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAg                                                               >  2:1548401/1‑74 (MQ=255)
                                                                                                                                               tGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAg                                                               <  1:1548401/74‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                       
CGCCGCAGGCATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCGT‑CTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGC‑AGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGG  >  NC_002947/6033680‑6033956

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